INTERNATIONAL 
JOURNAL OF ABSTRACTS 


Si wip iGA lb THEORY 
AND METHOD 


VOLUME 2 + NUMBER 3 


1961 


PUBLISHED FOR 
THE INTERNATIONAL STATISTICAL INSTITUTE 
OLIVER AND BOYD 


VOLUME 2 + NUMBER 3 


COPYRIGHT © 1961 INTERNATIONAL STATISTICAL INSTITUTE 


Abbreviation (World List Style) : 
Int. F. Abstr., Statist. Theory 


Short Title: 
Statistical Theory Abstracts 


PRINTED IN GREAT BRITAIN BY 
OLIVER AND BOYD LTD., EDINBURGH 


INTERNATIONAL JOURNAL OF ABSTRACTS 
STATISTICAL THEORY AND METHOD: 


COVERAGE OF JOURNAL 
THE aim of this journal of abstracts is to give complete coverage of papers in the field of statistical 
theory and new contributions to statistical method. Papers which report only applications or examples 
_ of existing statistical theory and method will not be included. ‘There are approximately two hundred 
and fifty journals published in various parts of the world which are wholly or partly devoted to 
the field of statistical theory and method and which will be brought within the scope of this journal of 
abstracts. A complete list of journals covered is printed in the annual Index Supplement. In the case of the 
following journals, however, being those which are wholly devoted to statistical theory—all contributions, 
whether a paper, note or miscellanea, will be abstracted : 

Annals of Mathematical Statistics 

Biometrika 

Journal, Royal Statistical Society (Series B) 

Bulletin of Mathematical Statistics 

Annals, Institute of Statistical Mathematics 
Within the larger group of journals, which are not wholly devoted to statistical theory and method, there 
are some journals which have the vast majority of their contributions in this field. These journals, there- 
fore, will be abstracted on a virtually complete basis : 

Biometrics 

Metrika 

Metron 

Review, International Statistical Institute 

Technometrics 

Sankhya 
After experience in the publication of this journal of abstracts it may be found desirable to add further 
journals to these lists. In any case readers may be assured that all papers properly to be included in this 
journal of abstracts will be included irrespective of such notification. If any reader of this journal discovers 
a paper which happens to have been overlooked, the General Editor will be pleased to be informed so 
that the appropriate abstract can be made: always provided that the date of publication is after 1st 
October 1958, when the abstracting for this journal commenced. 

In addition to the ordinary journals, there are two kinds of publication which fall within the scope 
of this journal of abstracts. They are the experiment | and other research station reports—which occur 
particularly in the North American region—atid “the | Teports of conferences, symposia and seminars. 
Whilst these latter may be included in the book review sections of journals it is unusual for any individual 
contribution to be noted at any length. These publications are, in effect, special collections of papers 
and for this reason the appropriate arrangements will be made for them to be included in this journal. 
By the same token, abstracts of papers given at conferences and reproduced in an appropriate journal will 
be disregarded until the definitive publication is available. 


FORMAT OF JOURNAL ; “i 
The abstracts will be up to 400 words long—the recommendation of UNESCO for the “ long 


abstract service: they will be in the English language although the language of the original paper will 
be indicated on the abstract together with the name of the abstractor. In addition, the address of the 
author(s) will be given in sufficient detail to facilitate contact in order to obtain further detail or request 


an offprint. Suitable indexes will be provided annually in a supplement. 
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The scheme of classification has been developed upon lines that will facilitate the transfer to punched 
cards of the code numbers allocated to each abstract: to allow for future development it is suggested that use 
is made of 4-column fields. Each abstract will have two classification numbers: the primary number in 
heavy type to indicate the basic topic of the paper and the secondary number in brackets to take account of 
the most important cross-reference. The sheets of the journal are colour-coded according to the twelve main 
sections of the classification and it should be noted that it is the main section number of the primary 
classification which determines the colour code for each abstract. It is believed that this method of 
colour-coding the pages will provide a distinctive visual aid to the identification of abstracts both when 
the journal is in bound form or dismantled and filed on cards or in binders. The format and simple 
binding allows of the following alternative treatments by users of the journal : 

(a) Leave intact as a shelf-periodical. 

(6) Split and filed in page form according to the main sections of the classification. 

(c) Split and guillotine (single cut) each page ready for : 

(i) pasting onto standard index cards. 
(ii) filing in loose-leaf or other binders—for which the appropriate holes are punched. 


It will also be possible for those users who need a completely referenced file to insert skeleton cards or 
sheets according to the secondary classification number provided on each abstract. 
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Volume 2. CORRECTIONS 


Author’s name should read, Emannelli, F. 
Language of original paper was Polish. 


AKAIKE, H. & SAIGUSA, Y. (Institute of Statistical Mathematics, Tokyo) 


0.6 (0.1) 


On a min-max theorem and some of its applications—Jn English 
Ann. Inst. Statist. Math., Tokyo (1960) 12, 1-5 (1 reference) 


In this paper the authors show that for a matrix with 
elements arranged in decreasing order in each row, we 
have only to select successively the rows with the smallest 
diagonal elements in order to rearrange the rows of the 
matrix, so that we obtain the minimum of the maximum 
diagonal elements of the rearranged matrices; that is, 
we have only to take as the new jth row the one with 
the smallest jth element in the remaining rows. 

Two examples of the application of this result are 
given: 

(i) application to the analysis of rational economic 
behaviour of a consumer who tries to keep 
minimum the maximum expenses in buying 
successively a set of commodities with gradually 
decreasing prices 


preparation to assure a good convergence rate in 
applying Seidel’s process of successive approxi- 
mation for the solution of simultaneous linear 
equations. 


(ii) 


(H. Akaike) 
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DANZIG, G. B. (RAND Corp., Santa Monica, Cal.) 


0.8 (9.1) 


On the significance of solving linear programming problems 


with some integer variables—In English 
Econometrica (1960) 28, 30-42 (10 references) 


This paper deals with the consideration that recent 
proposals by Gomory [“ Essentials of an algorithm for 
integer solutions to linear programmes,” Bull. Amer. 
Math. Soc. (1958) 64] and by Beale [“ A method of 
solving linear programming problems when some but 
not all of the variables must take integer values,” Statist. 
Tech. Res. Group (1958) Princeton] for solving linear 
programmes involving integer-valued variables appear 
sufficiently promising that it is worthwhile to systemati- 
cally review and classify problems that can be reduced 
to this class and thereby solved. 

The author commences by discussing the “ cutting 
plane approach ”’ which was first proposed and demon- 
strated by Fulkerson, Johnson & Danzig, G. B. [‘* Solu- 
tion of a large scale travelling salesman problem,” J. 
Operat. Res. Soc. (1954) 4, 393-410] and further explored 
by Manne & Markowitz [“ On the solution of discrete 
programming problems,” Econometrica (1957) 25, 19}. 
Gomory developed a theory of automatic generating 
“cutting planes” and the theory was generalised to 
the case where some variables are constrained and some 
continuous by Beale. 
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The second section deais with the fixed charge problem 
and is followed by the application of the principles to 
the orthogonal latin-square problem: in this context 
the author remarks on the fact that Bose & Shrikhande 
recently proved that Euler’s famous conjecture about 
the non-existence of orthogonal latin-squares of certain 
even orders was false [Notices of the American Math. 
Soc. (1959); abstract No. 558-27]. 

The fifth section discusses the classic problem of 
colouring a map using at most four colours so that no 
two regions having a boundary in common have the 
same colour, but the author states that ““ The two ways 
that are given to colour constructively a particular map 
do not contribute anything to a proof or the truth or 
falsity of the conjecture except that an efficient way 
for solving particular problems on an electronic com- 
puter may provide a counter example ”’. 

Historically, non-linear, non-convex and combinatorial 
problems are areas where classical mathematics almost 
always fail and it is therefore significant that a reduction 
to class of linear programming involving integer-valued 


variables can be made. 
(W. R. Buckland) 
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FERRERI, C. (School of Statistics, Palermo) 


0.2 (—.-) 


On a mathematical theory of the survival and the fit of some new functions 


to the Sicilian population—ZJn Italian 


Ann. Fac. Econ. Com. Palermo (1958) 12, 1-70 (18 references, 13 tables, 5 figures) 


Quiquet’s theory on survival functions is generalised. 
The author shows the practical impossibility of choosing 
rationally a function representing the survivals for all 
ages of a population. He later supposes a subdivision 
of the survivor’s distribution into the following three 
periods: (a) childhood, (6) middle age, (c) old age. 
He assumes, as discriminant values, the years corres- 
ponding to the minimum and to maximum of death. 

The author shows that the Lazarus function may 
describe the middle age better than many other functions; 
for childhood and old age he proposes, respectively, the 
new functions 


I(x) = As*g*’b** and I(x) = A(w—x)*b* 


where the parameters are posed with regard to bio- 
logical, sociological and demographic knowledge of 
human mortality. 
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FORTET, R. (University of Paris) 


After some critical remarks on Vincent’s criterion 
determining the oldest age, the author using the 
mean-square method fits his functions to the survival 
distribution of the Sicilian population. He gets a new 
death table in general agreement with the empirical date. 


(F. Galantino) 


0.0 (0.8) 


The algebra of Boole and its applications in operational research—IJn French 


Trab. Estadist. (1960) 11, 111-118 (9 references) 


Boolean algebra is defined as having two elements, say 
0 and 1. Boolean addition and multiplication are 
defined. Boolean functions are defined as functions of 
Boolean variables, that is, variables which can only 
take values 0 and 1. Boolean equations are defined 
as f(x, .--, X,) = 0, where f is a Boolean function and 
x, a Boolean variable. A set of equations of n Boolean 
variables is obtained as a solution for a condition 
imposed on the n Boolean variables. A method to 
solve systematically a condition C expressed by an 
equation of the type f(x, --., x,) = 0 is given. 

A function f(x,, --- X,) of 2 Boolean variables is said 
to be a whole algebraic function if it can be expressed 
as a polynomial in the sense of the ordinary algebra 
on the x; and x, A whole algebraic function can 
always be transformed in a linear function in the sense 
of the ordinary algebra if convenient supplementary 
variables and conditions with whole coefficients in the 
initial and supplementary variables are introduced. 

An algebraic Boolean programme is defined as 
follows: if x1, -»-) X, are Boolean variables for which 
a set of conditions R holds and f(x,, .--, X,) is a whole 
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algebraic function of the x;, we have to find the elements 
of the set Y,) in which the x; are defined that maximise 
f and satisfy R. Boolean programmes are very frequent 
in nature but they are very difficult to solve. A 
systematic method of proving all the 2” elements of V(,) 
is given in order to solve a Boolean programme. An 
example is studied. Equalities and inequalities for whole 
algebraic functions are studied. For these we have the 
following theorem: every algebraic Boolean programme 
where the condition R is expressed by algebraic equa- 
tions or inequations can be carried to a system of whole 
algebraic equations or inequations. It is possible to 
solve a whole algebraic equation or inequation by 
carrying it to an algebraic Boolean equation, at least 
when the numerical coefficients are whole numbers; 
this is enough for practical applications. Two examples 
are studied. 


(M. Melendez) 
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GREVILLE, T. N. E. (Office of the Quartermaster General, Washington) 0.6 (6.1) 
Some applications of the pseudo-inverse of a matrix—In English 
SIAM Rey. (1960) 2, 15-22 (13 references) 


An mxn matrix A of rank r>0 can be expressed as a 
product A = BC where B is an mxr matrix and C is 
an rxn matrix, and both are of rank r.. Then the 
pseudo-inverse of A is given by A” =C1(CC’)"! 
(B’B)~'B"™ where the superscript T denotes transpose. 
The author describes certain advantages in using the 
pseudo-inverse for obtaining orthogonal polynomials 
over a discrete domain. He also derives a recursive 
algorithm for the pseudo-inverse which is useful in 
fitting polynomials of successively higher degree. Atten- 
tion is drawn to the fact that ‘‘ pseudo-inverse ” is the 
same as the “conditional inverse’’ as used by Bose 
(Department of Statistics, University of North Carolina). 
This idea has widespread application in the general 
theory of linear estimation. 


(W. T. Wells) 
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HANSON, M. A. (University of New South Wales, Sydney) 0.8 (7.5) 
Errors and stochastic variations in linear programming—ZJn English 
Aust. J. Statist. (1960) 2, 41-46 (2 references) 


This paper points out that small errors in the data of 
a linear programme may cause a large variation in the 
desired maximum. Such errors should therefore be 
closely investigated, as their practical effects may out- 
weigh the advantages of linear programming. 

Problems of linear programming amount to determin- 
ing an n-dimensional vector x (with non-negative 
elements) subject to the condition Ax = b where A is 
an mxn matrix and b an m-dimensional vector, in 
such a way that the scalar product ex is maximised. 

The effects of fluctuation in the coefficients of A, 
b, c, and of the wrong choice of basic variables on the 
maximum of cx are examined. 

A simple example is given in which the coefficients 
of A alone are subject to error; it is shown that a co- 
efficient of variation of two per cent. in these causes an 
increase of 50 per cent. in the maximum. The author 
concludes “‘ that one should not accept the results of 
linear programming without an analysis of the errors 
involved ”’. 


(J. Gani) 
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JIRINA, M. (Math. Inst., Czechoslovak Academy of Sciences, Prague) 


0.1 (10.6) 


Ordinary differential equations with random coefficients and random 


right-hand side—IJn English 


Second Hungarian Mathematical Congress (1961) 


The author is concerned with the differential equation 


r 
arc.) ty. 


m=0 

We shall suppose that c,, are complex-valued random 
variables. The smallest o-algebra generated by them 
will be denoted by t. A wide-sense stationary random 
distribution will be called t-stationary, after Ito, if the 
wide-sense stationarity condition is required for con- 
ditional expectations with respect to t instead of ordinary 
expectations. Under the condition that y is a t- 
stationary random distribution, the necessary and 
sufficient condition for the existence of a t-stationary 
random distribution x which satisfies the equation is 
presented. 

The general form of x is described by means of its 
spectral decomposition. This makes possible the descrip- 
tion of the general form of the corresponding correlation 
distribution and spectral measure. 


(M. Jitina) 
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KNODEL, W. (Technical High School, Vienna) 
A generalised transportation problem—Jn German 


0.8 (-.-) 


Unternehmensforschung (1960) 4, 1-12 (5 references) 


A generalised transportation problem is treated in which 
each supplier has the possibility of producing a number 
of different products. An optimum solution minimising 
total cost must give the optimum shipment plan and 
state at the same time how the production capacity of 
each supplier is used for the different products. This 
problem can be treated as a linear programming problem 
and solved by the simplex technique. To avoid the 
computational burden, it is shown how this problem 
can be reduced to the ordinary transportation problem 
which can be solved by the stepping-stone method. 
This can be done by the two methods of introducing 
fictional transportation costs and localities. 

In the first case, the production of different com- 
modities by a single factory is treated as different supply 
points. Similarly the points of destination are dealt 
with. No delivery x,, from supply point P, to a 
destination B, is allowed, thus x,, = 0, ifu + v. If the 
stepping-stone method is applied, it seems rather 
difficult to take account of this condition in the actual 
calculations. However, if the corresponding unit trans- 
portation cost a,, is selected as a,, = 00 or, to allow 
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numerical calculations, is given a sufficiently high finite 
value, this route will not appear in the optimal solution 
thus fulfilling the above condition. In the second case, 
if total production capacity surmounts total demand, 
inequalities appear. These can be converted into 
equalities by introducing slack variables. This can be 
interpreted as adding a fictive point of destination By to 
which the amounts of goods not produced are being 
shipped at zero transportation costs. The total amount 
By equals total production capacity minus total demand. 

The method is applied to two special types of problems: 
in the first one each supplier can produce any combina- 
tion of different products within his capacity limits. In 
the second one, it is assumed that each supplier can 
produce any one of the commodities, but only one at 
a given time. Finally the more general case is treated, 
in this case a producer is only willing to devote part of 
his capacity to products of comparatively low profit 
margins. A formal proof for this method is given, 
this can also be applied to similar linear programming 
problems of a more general type. 


(H. Giilicher) 
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PATTERSON, H. D. (Rothamsted Experimental Station) 


A further note on a simple method for fitting an exponential curve—Jn English 


Biometrika (1960) 47, 177-180 (5 references) 


The author is concerned with the fitting of the exponential 
regression curve «—fp*(0<p<1) which he has discussed 
previously [“‘ A simple method for fitting an asymptotic 
regression curve’? Biometrics (1956) 12, 323-329 and 
later Patterson & Lipton “ An investigation of Hartley’s 
method for fitting an exponential curve”? Biometrika 
(1959) 46, 281-292; abstracted in this journal No. 1/602, 
4.3]. 

In this paper he presents a modification of his approxi- 
mate method which leads to increased efficiency without 
adding considerably to the arithmetic. Up to twelve 
equally spaced values of x are considered, and the 
estimation of e involves simply the root of a quadratic 
equation whose coefficients are given by linear functions 
of the ordinates. It is suggested that unless p is very 
small and the number of values of x very large this 
method will give estimates of the parameters which are 
of high efficiency and which are close to the least- 
squares estimates. Numerical examples are given. 


(Florence N. David) 
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PIEHLER, J. (Leuna-Werke, Merseburg, Germany) 
A note on sequencing problems—J/n German 


0.8 (-.-) 


Unternehmensforschung (1960) 4, 138-142 (3 references) 


Under certain restrictions the sequencing problem for 
m jobs and n machines is reduced to the travelling- 
salesman problem. These problems are concerned with 
situations where each of m jobs requires time on some 
or all of n different machines. The times required by 
the m jobs differ from machine to machine. 

The problem is to determine a sequence of the m 
jobs; a permutation of the integers 1, 2, ..., m, which 
minimises the total elapsed time. In the following case 
it is supposed that all jobs require time on each of the 
machines and that the sequence of machines must be 
the same for all jobs. Later it is supposed that there 
are no possibilities of intermediate storing. 

Let a,, be the time required by job i on machine k. 
Then the total elapsed time is given by 
m m-1 k el 
es ial ee » max pa in my ois 
k=1 Baa Sk Sty — el y=1 ¢ 
It is to be noted that Z is to be evaluated for any 
sequence of jobs. In order to minimise Z, or to find 
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the sequence minimising Z, one must determine the 
matrix 


k k-1 
max <>, 4,— >) a,,- forj#? 

A;; = 1Sk<m([vy=1 van 
M forj=i 


and solve m travelling-salesman problems, given by the 
matrix A;;, so that in each of these travelling-salesman 
problems one of the m jobs is treated as the last one in 
the sequence. 


A numerical example illustrates the method presented. 


(W. Dinkelbach) 
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VERHAGEN, A. M. W. 


Aust. J. Statist. (1960) 2, 122-127 (10 references) 


The paper starts with an a priori justification and 
introduction of Gompertz and logistic curves in dealing 
with growth data, and makes a comparison of the 
shapes of these curves. A close approximation may be 
obtained by matching the branches of the Gompertz 
and logistic curves on either side of the inflection point 
by exponential functions. Functional transformations 
can be made on the Gompertz and logistic functions 
and the three parameters estimated. The success of 
the Gompertz and logistic curves as graduation formulae, 
despite operational difficulties, appears due to the fact 
that, irrespective of the parameter values, they have 
the sigmoid shape characteristic of most growth processes 
in addition to having three adjustable parameters 
available for the matching of the observations. The 
Gompertz curve and the logistic curve are discussed in 
relation to the allometric and linear transformations. 
The rest of the paper is devoted to the investigation 
of the general functional form rH(p+qt) in relation to 
these transformations. It is established that the class 
of exponential functions £ exp (yf) is the only one with 
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WOLFE, P. 


(Commonw. Sci. Industr. Res. Organisation, Melbourne) 
Growth curves and their functional form—Jn English 


0.2 (10.9) 


the general form rH(p+qt) such that any two members 
are related by an allometric transformation. The 
exponential function does not lose its form under any 
allometric transformation. The Gompertz does not 
lose its form if the constant term y, is zero. 


(A. M. W. Verhagen) 


0.8 (-.-) 


The simplex method for quadratic programming—ZIn English 
Econometrica (1959) 27, 382-398 (9 references, 2 figures) 


After defining “‘ quadratic programming” the author 
states four problems which are of some interest: regres- 
sion, efficient production, the “ portfolio”? problem, 
and convex programming. The main difference between 
the method proposed by the author and other methods 
for the computational solution of quadratic program- 
ming problems lies in his use of the simplex method for 
linear programming. This is developed in two forms: 
the ‘‘ short ” and the ‘‘ long ” method, the latter arising 
on solving the following quadratic problem for all A420: 
minimise f(/, x) = Apx+4x'Cx 

where x20 and Ax = b: C is an n by n symmetric 
matrix, and p a one by n matrix. 

After discussing some preliminary points the author 
outlines the computation in section three; followed by 
an example in both the forms mentioned above. Section 
five contains the relevant proofs for the reversions given 
in section three and the final section modifies the com- 
putations for other types of constraints. 


(W. R. Buckland) 
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CLUNIES-ROSS, C. W. & RIFFENBURGH, R. H. (Virginia Poly., Blacksburg, 1.4 (6.4) 


and Univ. Hawaii) 
Geometry and linear discrimination—Jn English 


Biometrika (1960) 47, 185-189 (1 reference) 


This paper is concerned with the problem of classifying 
an individual into one of two multivariate populations. 
The discriminator is minimax and is assumed to be 
linear in the measurements; that is to say it is assumed 
to be a hyperplane in the multidimensional space. 

It is shown that provided an equation (authors’ 4.1) 
has a negative root, the optimum discriminant hyper- 
plane passes between the centres of the two populations, 
that is to say that the optimum discriminator dis- 
criminates between the means of the two populations. 

When the equation (4.1) does not have a negative 
root there is no clear-cut result. The best discriminating 
hyperplane will not pass between the population centres 
and the minimax procedure favours the population 
with the greatest weight factor. 


(Florence N. David) 
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FADDEYEV, D. K. & VOROBJEV, N. N. (Leningrad) 1.1 (-.-) 
Continualisation of conditional probabilities—Jn Russian 
Teor. Veroyat. Primen. (1961) 6, 116-118 


The probability measure p on finite set R is called 
interior if u(a)>0 for any ae R. The set of all interior 
measures on R is denoted by W(R). The following 
theorem is proved: 

there exists a mapping ¢ of W(R) into Euclidean 
space E of suitable dimension with two properties: 


(i) All conditional probabilities 
Rom) 6 2 : 
p(a | A) = —— (ae AcR) 
H(A) 


are uniformly continuous functions of ¢u on the whole 
set @W(R) in the sense of the metric on E 


(ii) The closure of ¢W(R) in E is homeomorphic to 
the closed simplex of suitable dimension. 


(D. K. Faddeyev) 
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HAJEK, J. (Math. Inst., Czechoslovak Academy of Sciences, Prague) 


1.5 (4.9) 


Limiting distributions in simple random sampling from a finite population—Jn English 
Publ. Math. Inst. Hung. Acad. Sci. (1960) 5, 361-374 (5 references) 


The author considers an infinite sequence of simple 
random sampling experiments without replacement, the 
vth of which has the size n, and refers to a population 
of size N, with values y,,,..., yyy, Let €, be the sum of 
yy; 8 in the sample. Necessary and sufficient conditions 
concerning {y,, m,, N,} are given under which the 
distribution function of &, converges to a limiting 
distribution law: n,>0o and N,—n,>0. This result 
contains as a special case that of Erdés & Rényi [Publi. 
Math. Inst. Hung. Acad. Sci. (1959) 4, 49-61; abstracted 
in this journal No. 1/171, 1.5] for the convergence to 
the normal distribution law. In addition it is shown 
that the condition given by them is not only sufficient 
but also necessary. 


If n,/N, does not converge to 0, normal limiting 


distribution and limiting distribution formed as a 
convolution of a normal distribution and some two- 
points distributions are possible. If n,/N, converges to 
0, the variance preserving limiting distribution may be 
only infinitely divisible. The conditions for this are 
the same as if the sampling were carried out with 
replacement. 
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HAJEK, J. 


(Math. Inst., Czechoslovak Academy of Sciences, Prague) 


The proof is based on a model in which an experiment 
producing simultaneously a simple random sample and 
a so-called “‘ Poisson sample ”’ is defined. In “‘ Poisson 
sampling ’’, each unit of the population has indepen- 
dently the probability n,/N, of being selected; the sample 
size is a random variable with expectation n,. The 
experiment is defined in the manner that the simple 
random sample should contain the Poisson sample or 


Ny 
vice versa. Let Y,= )) y,; and 4 and n* be the sums 
i=1 


of the differences y,;— Y, for the simple random sampling 
elements and for the Poisson sampling elements, 
respectively. Then the following inequality holds: 


%*\2 
E(n—n*) < Lest) 
Dn* n, N,-n, 


As n* can be regarded as the sum of independent terms 
the main theorems can be proved by application of the 
well-known theory of independent random variables. 


(K. Sarkadi) 


1.5 (4.9) 


Necessary and sufficient conditions for asymptotic normality of a scalar product 


of two randomly permutated vectors—IJn English 


Second Hungarian Mathematical Congress (1961) 


In this paper the author presumes a sequence of couples 
of vectors with real components: 


Ree) rie, Bun) | 


Let (X,,, ..., Xyy,) be a random vector that takes on 
the N,! permutations of (4,1, ..., dyy,) with equal 
probabilities. Consider the statistic 


19.0) 
vy 


=1 


Ny 
fT; = pak Dy: X yi 


and search for conditions of asymptotic (v>oo) nor- 
mality of T, with respective parameters &(7,) and D(7)). 

This problem was first formulated by Wald & 
Wolfowitz in 1944 who also derived a sufficient condition 
which was strengthened by Noether in 1949. The most 
far-reaching result was achieved by Motoff in 1957 who 
proved that the Lindeberg type is sufficient in fact, 
not only in the problem under consideration, but also 
in its generalisation by Hoeffding. The Lindeberg type 
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condition is also necessary, however, as is shown in 
the following theorem: 


provided that the a’s and the b’s satisfy a particular 

limiting condition, the statistic 7, is asymptotically 

normal with parameters &(T,) and D(7,) if, and only if, 
yy \ (iz d,)*(b,; —b,)? 


lim (i, J) € Svc 
Ny N 


Ney (Gyranhy 6) 
1 


tS ge 


=0 


for every t = 0, where 
Sye a } (i, j) : Gi G,)*(b,,—b,)* 


Inks Ny 
=a aliNe yy (a,;—a,)* » (ba bts 
i=1 i=1 


(J. Hajek) 
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HASMINSKI, R. Z. (Moscow) 1.5 (10.1) 
On limit distributions for sums of conditionally independent random variables—Jn Russian 
Teor. Veroyat. Primen. (1961) 6, 119-125 (13 references) 


This paper deals with limit distributions for sums », 
which become independent when a certain path x,(n = 0, 
1, 2, ...) of a Markoff chain is defined. The dependence 
between {y,} and {X,} is expressed more exactly by 


Pr{n, <Z,, N2<Z2, sieey Nn<Zn | Xo, X,, skola Ags 
= Prin, <Z, | Xo, X1}Prin2<x2| Xo, X1, X2}x 
ee eS eee ILA os A aeissh al a} 


Bsr 


Let X, be the path of a continuous Markoff process. 
It is shown that the study of the limit distributions of 


of limit distributions of sums y,. This reduction is 
illustrated for the case where X, is a single-dimensional 
diffusion process. The limit distribution for ¢(t) coin- 
cides with distributions stated by Feller [Trans. Amer. 
Math. Soc. (1949) 67, 98-119]. 


| 
t 

f= | S(X,)ds as t>oo can be reduced to the study 
19) 


(R. Z. Hasminski) 
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LINNIK, Y. V. (Steklov Institute of Mathematics, Leningrad) 1.5.2:3) 
Some new limit theorems for the sums of independent random variables—Jn English 
Second Hungarian Mathematical Congress (1961) (2 references) 


This is a somewhat extended version of the report sent 
to the fourth Berkeley Symposium, and briefly reported 
elsewhere [ Dokl. Acad. Nauk, SSSR (1960) 133]. The 
zones of integral and normal convergence for big 
deviations were considered and Petroff [Dokl. Acad. 
Nauk, SSSR (1960) 134], using a suitable extension of 
the author’s method, succeeded in obtaining “ satis- 
factory’’ necessary and sufficient conditions for the 
convergence to the set of “ Cramérian tails”. The 
above-mentioned results on normal convergence follow 
from the last results as a particular case. 

The problem of finding the conditions for convergence 
for big deviations to the set of laws other than stable 
laws and “‘ Cramérian tails ”’ still remains open. 


(Y. V. Linnik) 
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MOGYORODI, J. (Eétvés University, Budapest) 


CCE ee 


1.5 (1.6) 


On limiting distributions for sums of a random number of independent 


random variables—In English 


Second Hungarian Mathematical Congress (1961) 


Let us consider a sequence {¢,' of independent random 
variables and put (,= €,+,+...+¢, Let {B,} bea 
sequence of positive numbers tending to infinity such 
that Pr(¢,/B,<x) tends to a proper distribution function 
F(x) as n>+oo. Further, let {v,} be a sequence of 
positive integer-valued random variables. We mention 
that nothing is supposed about the dependence of v, 
on the random variables €,. 
The following two theorems hold: 


Theorem 1. Let us suppose that v,/n converges in 
probability to one as 2+ 00. In order that 


[C,,/Bn)<X] 
should converge to the law F(x), it is necessary and 
sufficient that the following condition 


fim lim sup B,/Braci +0} 

67-0 no 
be satisfied. This theorem is in a certain sense a 
generalisation of a theorem by Anscombe [Proc. Camb. 
Phil. Soc. (14952) 48, 600-607}. 
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NAGAEYV, S. V. (Tashkent) 


More exact statements of limit theorems for homogeneous Markov chains—Jn Russian 
Teor. Veroyat. Primen. (1961) 6, 67-86 (11 references) 


This paper contains several theorems defining more 
precisely the convergence of homogeneous Markoff 
chains with an abstract state space to the Gaussian 
distribution. Moreover, limit theorems for large 
deviations are proved. The proofs are analogous to 
the ones in an article of Nagaev [Teor. Veroyat. Primen. 
(1957) 2, 389-416]. 


(S. V. Nagaev) 
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Theorem 2. Let ¢,=€,+¢,+...+&, where the 
random variables €, are independent and identically 
distributed with mean value zero, variance one and let 
&(€{)<A: A is independent of k. Let us suppose that 
v,/ converges in probability to a positive random 
variable. Then 


: ce 1 ee nes 
lim Pr{ 2” <x} = Ee.” ay, 
n>+c ea ) \/(2n) ie Z 


This theorem is a generalisation of a theorem by Rényi 
[Acta Math. Acad. Sci. Hung. (1960) 11, 97-102; 
abstracted in this journal, No. 1/557, 1.5]. 


: 
(J. Mogyorddi) 
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PALASTI, I. (Math. Inst., Hungarian Academy of Sciences, Budapest) 


On some random space filling problems—ZJn English 


Second Hungarian Mathematical Congress (1961) 


Place as many as possible disjoint domains at random, 
and with uniform distribution, into a rectangle with 
sides x and y: these disjoint domains being congruent 
and parallel to a given convex domain D. The question 
arises, how a large part of the rectangle will be filled on 
average with these domains? 

Let us denote by M(x, y) the mean value of the 
number of domains placed in the above-mentioned 
way. The plausible hypothesis is set up that there 
exists a constant B>0 such that 

| M(x, +x2, y)— M(x1, y)— M(x2, y) | S By 
and | M(x, VitV2)—- M(x, y)—-M(, Ya) | S Bx. 

A theorem of Hyers concerning functional inequalities 
is generalised for two dimensions as follows: let f(x, y) 
be a Borel-measurable function of two variables, satis- 
fying the following conditions: 
| Aor +, YS, Y)—fl%2, y) |S By for x, 20, x2 20, 
y20and | f(x, ¥1+¥2)—f% 1) -f% Y2) | S$ Bx for x20, 
yi, 290, y. 20. 

B>0 is a constant, not depending on the variables 
xX; ¥;, (i = 1, 2) and there exists a constant p>O such 
that f(x, y) = 0 if 0<x<p or 0<y<p then 

| fx, y)—axy | $BOx+y) 
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PREDETTI, A. (Catholic University of Milan) 


EE EE 


1.4 (1.5) 


where a is constant. Thus 
lim f(x, y)/xy = a 


x,yro 
with the aid of which the limit of the above hypothesis 
is shown to be 

lim M(x, y)/xy = a(D) 


x, yo 
If Dequals the unit square S,, then «(.S,) = c? where the 
constant c is the constant obtained in the analogous prob- 
lem for one dimension which has been solved by Rényi 
[Mag. Tud. Akad. Kut. Mat. Intézet. Kézlemenyei (1958) 
109-127: abstracted in this journal No. 1/18, 1.4], 


ea t4_,-4 
c=2| (1-e~‘) exp(-2 | mane iu) dt ~0:748 
0 ty) u 


and thus the value of c? is approximately equal to 0-56. 
This means that about 56 per cent. of the area of a 
large rectangle will be covered by randomly placed unit 
squares, with sides parallel to those of the rectangle. 

Analogously, it isconjectured that, for then-dimensional 
problem, on average the c"th per cent. of the volume 
of a large n-dimensional parallelotope will be covered by 
n-dimensional unit cubes, with sides parallel to those 
of the parallelotope, and placed at random. 

(I. Palasti) 


12000 


Characteristic features of a general scheme of repeated trials with 


dependent probabilities—Jn Italian 


Riv. Ital. Econ. Demogr. Statist. (1960) 14, 85-93 (2 references) 


The author considers a particular scheme by Woodbury 
according to which the law of dependency of the 
probabilities is linear. 

After having proved that this scheme can be con- 
sidered as resulting from the application of a Polya 
scheme to a Poisson scheme, he determines the mean 
and variance of the number of successes. Jn particular, 
it is shown that the said mean value has an expression 
formally analogous to the compounded amount, principal 
value and interest, of a rent. 

The author then gives a procedure by which the 
parameters included in the formulae obtaining for the 
mean and variance can be determined approximately. 


(V. Levis) 
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Some results associated with Bochner’s theorem—Jn Russian 
Teor. Veroyat. Primen. (1961) 6, 87-93 (9 references) 


In this paper the authors discuss some results associated 
with Bochner’s theorem. Let {P,} be a family of 
probability distributions in a separable Hilbert space; 
or more generally, in a space ¥ = Y* conjugate to 
countably-Hilbert space Y: let {y,} be the family 
of corresponding characteristic functionals. They in- 
vestigate whether or not there exists a locally convex 
topology T having the property that the relative com- 
pactness of {P,} is equivalent to uniform, with respect 
to «, continuity of {y,}. 

They prove that a topology of this kind does not exist 
except for the case of countably-Hilbert nuclear space Y. 


(Y. V. Prohoroff) 


2/449 


RENYI, A. (Math. Inst., Hungarian Academy of Sciences, Budapest) 
Limit theorems concerning random walk problems—J/n Hungarian 
Magy. Tud. Akad. III Oszt. K6zl. (1960) 10, 149-169 (20 references, 1 figure) 


In this paper the author summarises the most important 
results in connection with the random walk; first of all 
on the real line. He gives new proofs of some known 
results. 

For example, he applies the generalised Borel-Cantelli 
lemma [Erdés & Rényi, “On Cantor’s series with 
convergent £1/q,,"’, Ann. Univ. Sci. Budapest, Sect. Math. 
(1959) 2, 93-100; abstracted in this journal No. 1/543, 
1.5] to prove the recurrence theorem of Polya. A simple 
proof of the arc sine law is obtained by showing that the 
characteristic function of the number of steps among 
the first n steps, after which the random walk point 
is to the right of the origin, is essentially equal to a 
Legendre polynomial. 


(P. Révész) 
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REIMANN, J. (Math. Inst., Hungarian Academy of Sciences, Budapest) 


1.3 (0.6) 


Matrix-theoretical treatment of random walk problems—Jn English 


Second Hungarian Mathematical Congress (1961) 


Random walk on a lattice-point region of an n-dimen- 
sional-space with absorbing and reflecting barriers, 
leads to hyper-matrices, which can be expressed by the 
direct products of the matrices obtained in the one- 
dimensional case. Thus, the results obtained for one 
dimension can be applied in the investigations of random 
walks of higher dimensions. 

With the matrix method outlined above, the following 
problems are solved: 


(i) in the case of random walk in a plane and in space, 
respectively, on a lattice-rectangle with absorbing 
barriers, the probability is calculated that a particle 
starting from a given interior lattice-point will arrive 
in steps to a predescribed interior or absorbing lattice- 
point, 

(ii) in the case of absorbing barriers on a n-dimen- 
sional rectangle with given size and given starting point, 
a formula is given for the probability of the particle 
being “‘ alive ” after N steps, 

(iii) in the case of symmetrical random walk, the 
expected value of the number of steps before absorption 
is calculated, 


2/451 


REIMANN, J. (Math. Inst., Hungarian Academy of Sciences, Budapest) 
Matrix-theoretical treatment of random walk problems—ZIn English 
Second Hungarian Mathematical Congress (1961) 


(vi) An iterative matrix-theoretical method is given 
for the solution of the problem mentioned in (iv) above, 
in the case of a lattice-point region of arbitrary form. 

(vii) In case of random walk on a lattice-rectangle 
with reflecting barriers, explicit formulae are given for 
the solution of the problem stated in (i). 

(viii) A special case of random walk in a plane is 
investigated; namely, when the probability of moving 
in different directions is not constant, as in above cases, 
but is a function of the coordinates. In connection 
with this problem, the Ehrenfest-model of diffusion is 


generalised. 


(J. Reimann) 
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(iv) for an asymmetrical random walk in two and 
three dimensions with absorbing barriers, the probability 
that a particle starting in an interior point will be 
absorbed after an undetermined number of steps in a 
predescribed absorbing point is calculated. That is, in 
case of the random walk in a plane or on a lattice 
rectangle consisting of m+2 rows and n+2 columns, 
this probability is 


JE (ip \Eae 
Ae ee. (LNG oe ls 
; m+1\p 5 / (rs) 


. oe iln oa pln sh(n+1—j)9,, 
1=1 m+] m+1  sh(n+1)9, : 


2ch9,= : ao (ee cost 
iGs) rs n+1 


This probability has also been determined for the 
symmetrical random walk by McCrea and Whipple. 

(v) The problems mentioned before are also answered 
in case of partly absorbing and partly reflecting barriers. 


where 
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REVESZ, P. (Eétvés University, Budapest) 


1.6 (10.1) 


A remark on the random ergodic theorem—ZJn English 


Second Hungarian Mathematical Congress (1961) 


In this paper the author poses a problem in connection 
with a theorem stated by Ulam & Neumann and later 
discussed by Rényi; the original theorem is as follows: 

let 7 be a measure space of the measure preserving 
transformations defined on {X, ¥, } with a probability 
measure P. Let us consider the product space 


T* = TXT Xo, 


where 7; = 7(i = 1, 2, ...) and the product measure 
p= => BP OSP Xie. 
Let f(x)e Lyx, ,) then 


Pr \tin Y) f(T... Tx) >f *(x) almost everywhere = | 
k=1 


where (7, T>, ...) is a point of Z* and f*(x)eLix. ys 
Rényi raised another problem: if there is a sequence 
T,, T>, ... of measure preserving transformations, under 


what conditions does 1/n > f(T,...T,x) converge 
k=1 


almost everywhere or in mean? An interesting problem 
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SARKADI, K. 
On random walk problems—In English 


(Math. Inst., Hungarian Academy of Sciences, Budapest) 


in connection with this question is the following: if 
from the set 7 a sequence 7, 7, ... is chosen at random, 
not necessarily with the same distributions; for instance 
with the distributions P,, P,, ..., but independently, 


under what conditions does 1/n ) f(T,...T,x)>f*(x) 
k=1 


hold almost everywhere or in mean with probability 
one? Some conditions are given: for instance, it is 
proved that if 


[| semper, 


| TeZ 


Sm || fitz) msi-1i-* 


ge B) 

0<e<l 
for every f(x) for which fe L?, | f= 0 then the result 
holds. It is used to prove a strong law of large numbers 
in the theory of Markoff chains. 


(P. Révész) 


1.3 (5.6) 


Second Hungarian Mathematical Congress (1961) (2 references) 


In this paper the author gives a new proof for a known 
theorem concerning the so-called Galton-test [see Feller, 
An introduction to probability theory and its applications 
(1957) New York: Wiley, 68-73] and some connected 
problems are considered. For the above-mentioned 
theorem a combinatorial proof [Hodges, “* Galton’s 
rank-order test’? Biometrika (1955) 42, 261-262] is 
already known; the proof given here is simpler. 

The combinatorial formulation of the mentioned 
theorems are as follows: 

let us consider the random sequence S of n (+1)’s 
and n(—1)’s. Let s, be the sum of the first k members 
of the sequence and let us denote by v = v(S) the 
smallest, and by w = w(S) a randomly chosen number 
from the values for which the relationship 


S, = Sy = Max S;, 
k 


holds. Further, let w = u(S) be the number of positive 


numbers in the sequence S;, 53, -++, San-1- 
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The theorem concerning Galton’s test states that 
Rr) ee (0, Ls 2). 
A theorem of Sparre Andersen states that 
PrGr= 127) p= er= ee 
The following theorem is also proved: 
Prwi= 7) al 2n. 


The above theorems are proved by giving one to one 
transformations mapping the sequences S into sequences 
S’ and S” for which u(S’) = u(S)—1, v(S”) = v(S)-1 
and w(S”) = w(S)—1 respectively. 

The third theorem can be extended for a sequence 
of n (1/n)’s and m (1/m)’s as well. In statistical termin- 
ology: the theorem can be extended for samples of 
unequal size. 


(K. Sarkadi) 
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TINTNER, G. (Iowa State University, Ames) 


1.8 (0.7) 


The application of decision theory of probability to a simple inventory problem—ZJn English 


Trab. Estadist. (1959) 10, 239-247 (12 references) 


The author uses an elementary form of the inventory 
problem in order to illustrate some ideas in the field 
of decision theory and probability. 

In order to illustrate this he supposes a commodity 
that is bought for $b, and sold for $b, per indivisible 
unit. Assuming some hypotheses about sales it is 
necessary to determine the optimum number of items 
to be bought. Three cases are discussed in this paper: 


(i) A probability distribution of sales exists and is 
known. The objective might be either the 
maximisation of the mathematical expectation of 
profit, or the minimisation of variance with equal 
mathematical expectation. 

The probability distribution is not known. The 
problem might be considered as a game between 
the individual and nature. The objective might 
be reached with either the minimax strategy or 
the minimum regret; regret equivalent to loss 
which is represented by the difference between 
the profits actually made and the maximum profit 


(ii) 
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VIETORIS, L. (Math. Inst., Universitat Innsbruck) 


which might have been realised by the given 
strategy. Alternatively one might maximise a 
weighted average of the maximum and minimum 
profit which can be realised. 

The probability distribution of the state of nature 
is not known and the optimum is obtained by 
use of an a priori probability distribution: either 
an equal probability is assumed for possible 
demand values, or an application of Carnap’s 
theory of probability is utilised for hypotheses 
about probability of possible demand values. 


(iii) 


(J. M. Garcia) 


1.0 (4.0) 


Remarks and estimates concerning induction—Jn German 
Monatsh. Math. (1960) 64, 233-250 (14 references) 


Technically, this is in essence a sharpening along lines 
suggested by J. V. Uspensky [/ntroduction to Mathe- 
matical Probability (1937) New York: McGraw-Hill], of 
Bernoulli’s (Ars conjectandi) estimation of probabilities 
through observed relative frequencies. 

More informally, the topic is discussed in context 
with a degree-of-confirmation concept of probability; 
no mention is, however, made of the standard literature, 
e.g., Carnap [The Continuum of Inductive Methods (1952) 
Chicago: University of Chicago Press]. While the 
author explains a degree of confirmation as “a relation 
between the accredited hypotheses and the accredited 
affirmation ” and is therefore in agreement with Carnap’s 
theory of c-functions, it appears plausible that quite 
different points of departure, like that of de Finetti or 
L. J. Savage would do equally well for Vietoris’ specific 
purpose. 

In particular, the author takes issue with Bertrand 
Russell’s position [Das menschliche Wissen, Zurich: 
Holle] regarding the justification of inductive inferences 
in principle. According to Russell some extra-logical 
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properties of the real world alone, not mathematical 
probability theory, must provide that justification. It 
is the author’s contention that extra-logical properties, 
called “primary” by him do indeed provide that 
justification, but via the relative frequencies which depend 
on “ physical constants”? and which, as “ degrees of 
possibility,” form the conventional basis of probability 
theory and statistics. 

Of more technical content is the author’s rejection 
of a probabilistic analogue of statement-logical contra- 
position (i.e., “if A, then B; then if non-B, non-A’), 
by means of an inverted use of the Clopper-Pearson 
formula [Biometrika (1934) 26, 404] for confidence 
limits for which tables have been provided by Dudin- 
Barkovski & N. V. Smirnoff [Teoria veroyatnostei i 
matematicheskaya statistika v tekhnike (1955) Moscow: 
Gittl. 


(E. M. Fels) 
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With unknown means p and g respectively. The problem 
is to maximise the expected sum of n independent 
sions, each being from X or Y, where after cach 
_ observation the experimenter may select the variate to 
- besampled. The author restricts himself to the following 
Situation: observe alternate pairs, a decision being 
made ai cach stage by watching the random walk 
—«O = X, — Y,+---+ X,— ¥, with absorbing barriers at 
_ ta. If —a<U,<2a obtains, observe another pair: 

f U,2zxs—-2) obtaims, stop observing pairs and 
' observe only X’s (¥’s) for the remainder of the n steps. 
Lat K be the step at which either absorption occurs or 
else K = p(n = 2p). 

The best possible outcome in expectation is no where 
o=max(p, g). The loss function is defined as 
L, =no—W, where W is the expected sum of all n 
steps. It is proved that L, = (o—1)[n+(—1)EK)/ 
(+1), where « = p+g—o and u = A1—72)/x1—«). 
A strategy consists simply of a choice of a positive 
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VOGEL, W. (Univ. Tibingen, Germany and Univ. Chicago) 
A sequential design for the Two-Armed Bandit—JIn English 


EL, W. (Univ. Tubingen, Germany and Univ. Chicago) 
_ A sequential design for the Two-Armed Bandit—In English 
Amn. Math. Statist. (1960) 31. 430-443 (& references) ar 


| Consider two independent Bernoulli variates X¥ and Y 


1.8 (9.7) 


integer asp. A strategy minimising L, cannot be 
found without some a priori knowledge of p and q. 
This is the approach of Bradt, Johnson & Karlin [Amnn. 
Math. Statist. (1956) 27, 1060-1074]. 

Using €Uy = (p—Q@)@K, a slight generalisation is 
proved in this paper and another formula, derived in 
terms of absorbing probabilities at each barrier and the 
probability of stopping at yp, the author obtains the 
inequality €Kso{u’—1)(o—sXu"+1) with the left- 
hand side approaching the right as no. Let M, be 
the approximate loss-function obtained by substituting 
the upper bound for @K into the expression for L,. 
The author then uses M, to choose an optimum strategy 
in the minimax sense. 

Now M, is a function of 2, c, « and by regarding a as 
a continuous variable on 0<a<y, the author obtains 
the theorem for n=4 that the minimum over =< of the 
maximum over (c, t) of M, equals the quantity attained 
by interchanging these extrema. Calling the saddle 
point coordinates a(n), o{n), <(n) the author gives their 


continued 
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Aun. Math. Statist. (1960) 31. 430-443 (8 references) 


_ asymptotic behaviour as n> <0 as well as o{n)—7(n). 
_.and a{n)logu{n). This allows the author to propose 
the simple formula a(n) = (0-292,/n) for n> 100 as the 
optimum minimax strategy. Some considerations are 
given as to when the approximation M, is good and 
“the result applicable. but no error estimates are given. 
Another approximation for L,, for n large. is obtained 
heuristically by considering a sequence of random walks 
of the type discussed, approaching a Weiner-Process, 
and by choosing the appropriate absorbing boundaries 
and using well-known formulae for the expected time 
to reach these boundaries. On this basis the author 
conjectures a limiting formula for L,/, in with 2, G,—T, 
‘being O(,/n) and oc, 1,~pAO. 1)- “The author then 
 adduces some evidence for his suspicion that this 
"limiting formula has a saddle point. In the last section 
some conditions are discussed which allow this procedure 
~ to be generalised to other than Bernoulli variates. 


(S. C. Saunders) 
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AMATO, V. (Faculty of Economics, University of Catania) 
A new theoretical scheme of a factorial type and the distribution of 


plural births—Jn Italian 


Riv. Ital. Econ. Demogr. Statist. (1960) 14, 33-64 (9 references, 14 tables) 


The author, noting the systematic difference between 
the empirical frequencies of some distributions of plural 
births and the theoretical frequencies derived from the 
geometric law, suggests a less simple model which 
appears to give a better “‘ fit ” to the actual distributions. 

Such a model belongs to the class of factorial models 
with dependent probabilities, and its choice is justified 
by the ample heterogeneity of the partial distributions 
entering into the total distribution of the observed 
phenomenon and by the dependency between the relevant 
events. 

The author first shows the better “fit” obtaining 
with the theoretical distribution he suggests and then 
proves how this distribution converges to the geometric 
one. A simple procedure is also given for calculating 
the moments. 


(V. Levis) 
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BARDWELL, G. E. (University of Denver, Colorado) 


On certain characteristics of discrete distributions—In English 


Biometrika (1960) 47, 473-476 (1 reference) 


The author starts from the remark, “‘ It has been stated 
that for the hypergeometric, binomial, Poisson and 
geometric distributions the relationship holds that the 
mean deviation is equal to twice the variance multiplied 
by the frequency function at [], where [j/] is the greatest 
integer not greater than the mean’”’. 

He is concerned to point out that this relationship 
is true for the binomial and Poisson, true for the geo- 
metric distribution for integral values of y, but does not 
hold for the hypergeometric distributions. Further 
properties of the discontinuous frequency functions 
mentioned above are derived. 


(Florence N. David) 
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2.5 (0.2) 


2.1 (2.4) 


BOROVKOFF, A. A. 


7a 


Local theorems and moments for maximums of sums of lattice restricted 


components—I/n Russian 


Teor. Veroyat. Primen. (1961) 6, 108-110 (1 reference) 


Let ¢,, €,, ..., be independent latticed random variables; 
eel Meso imax sre, 4.) +€,;). 


1Svsn 
_ Both the formulae for Pr(S, = x) and for the first 
moments of s, are obtained by the author in this paper. 


(A. A. Borovkoff) 
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BUEHLER, R. J. 


(Statist. Lab., lowa State Univ., Ames, Iowa) 


2.2 (6.8) 


An application of regression to frequency graduation—In English 
Biometrics (1960) 16, 659-670 (11 references, 3 tables) 


Given observed values (x) grouped in a frequency 
distribution, it is desired to derive a function (x) 
which gives a smooth representation in terms of some 
standard distribution g(x). Hammersley & Morton 
[Biometrika (1954) 41, 296-301] utilised a method of 
weighted least-squares to determine the parameters « 
and f in the transformation ¢(x) = «+x, where g(x) 
~ is arbitrary. 

In this paper the author generalises this transformation 
to 


Wx) = Y Bible 

where the #, are any number m of arbitrary specified 
functions. The f; are asymptotically efficient estimates 
obtained by the method of weighted least-squares pro- 
posed by Aitken [Proc. Roy. Soc. Edinb. (1935) 55, 
42-48]. Simultaneous asymptotic minimisation of chi- 
square is also noted. A major advantage of the method 
is that estimates of the percentile points y; of g(x) are 
given by an explicit formula rather than the implicit 
formula used by Edgeworth [J. R. Statist. Soc. (1898) 
61, 671-699]. 
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The author proposes an F test which is helpful in 
determining the correct number of parameters in the 
transformation. However, the use of a multi-parameter 
expression of the form shown above has the inherent 
difficulty that @ may have maxima or minima which 
leads to a multi-valued inverse function 1/¢. The 
problem is somewhat resolved if it can be shown that 
¢@ is monotone in the range of interest. 

As an example, the method was applied to the 
breadth-of-bean data of Johannsen reported for example 
by M. G. Kendall & Stuart [The Advanced Theory of 
Statistics (1958) London: Griffin] to effect a cubic 
transformation (m = 4) to normality. Details of the 
calculations are shown. Graduated percentile points 
obtained by this method are compared with fitting by 
three other types of four-parameter equations. Both 
the Pearson type IV curve obtained by Pretorius [Bio- 
metrika (1930) 22, 109-223] and that obtained by 
Johnson [Biometrika (1949) 36, 149-176] using a type 
Sy translation give superior fits using y? as a rough 
guide. For the purpose of transforming to normality 
the method proposed in this paper is superior to that 
of unweighted least-squares. 

(J. E. Dunn) 


COHEN, A. C. (University of Georgia, Athens, U.S.A.) 


EE <<< ee 


2.8 (4.3) 


An extension of a truncated Poisson distribution—Jn English 
Biometrics (1960) 16, 446-450 (6 references, | table) 


This paper is concerned with the estimation of para- 
meters for a truncated Poisson distribution with missing 
zero class. Although the distribution to be considered 
might be studied from a strictly abstract point of view, 
‘the author points out that it is desirable to keep in 
mind the practical example of the distribution of 
organisms among colony sites under circumstances 
such that no migration occurs between colonies. 

The maximum-likelihood estimate is obiained for one 
parameter. The maximum likelihood estimating equa- 
tion is obtained and easily solved for the second para- 
meter by interpolating linearly in Table I of an earlier 
article by Cohen [Biometrics (1960) 16, 203-211; 
abstracted in this journal No. 2/46, 4.3]. 

The asymptotic variance-covariance matrix of the two 
parameters is obtained and the parameters are shown 
to be asymptotically independent. 
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CRAWFORD, J. R. & WALSH, J. E. (System Development Corp., Cal.) 


The results of this paper are illustrated by using an 
example from Beall & Rescia [Biometrics (1953) 9, 354- 
386]. A tabulation of the observed data is presented 
together with expected frequencies using results of this 
paper, and for comparison expected frequencies based 
on a complete Poisson distribution and on five con- 
tagious distributions considered by Beall & Rescia. 
As modified here, the superiority of the truncated Poisson 
distribution is clearly demonstrated in fitting the observed 
distribution of this example. 


(L. J. Walker) 


2.7 (2.1) 


Empirical examination of Edgeworth series—Jn English 
Ann. Inst. Statist. Math., Tokyo (1960) 12, 13-26 (4 references, 3 tables) 


In this paper the authors describe a numerical in- 
vestigation carried out for five types of probability 
approximations by using the first seven terms of the 
Edgeworth series expansion for the distribution of a 
continuous random variable 7. These approximations 
are applied to three types of probability expressions 
Wj i = 1), FO)—h ep) = Pr(—fSTSt), and 
. F(t)—F(—t+1) = Pr(—t+1ST<z2), where T is taken 
to have zero mean and unit variance. 
Calculations are performed for each combination of 
certain values of the third, fourth, and the fifth moment 
about the mean. For those ¢ values in the set —4-00 
(0-25) 4:00, the authors obtain significant limits 7,, t, 
(t;<t,) which indicate that a probability expression 
satisfies the required monotonicity conditions and neither 
exceeds unity nor falls below zero for computed ¢ such 
that t;<t<t,. In addition to the significant limits, the 
values of F(O) and F(1:75)—F(—1-75) are given for 
each case investigated: the tables of these results are 
also given. 


(M. Huzii) 
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GAEDE, K. W. (Technische Hochschule, Miinchen) 


2.0 (2.2) 


On the distribution of the roots of random algebraic equations—In German 
Math. Nachr. (1960) 21, 81-107 (6 references, 10 figures) 


In this paper the author discusses a polynomial of 
degree n in one variable x with the leading coefficient 1: 
the other coefficients A,, A>, ..., A, are of the form 
A, = S,+iT,, the S, and 7, being real random 
variables. It is supposed that the joint distribution 
function of these 2n random variables is known and 
that the density function exists. The roots are suitably 
arranged in order to get a one-to-one correspondence 
between the coefficients and the roots of the polynomial. 
The functional determinant of the related transformation 
can be explicitly computed. With this, the joint 
frequency function of the real and imaginary parts of 
the roots can be calculated. 

In the case where all random coefficients are always 
real, they cannot be submitted to the above-mentioned 
research because of the supposed existence of the density 
function. It is assumed that the joint density function 
of the 1 real random variables S,, ..., S,, exists and is 
known. The frequency function of the roots can 
also be computed; but some cases have to be distin- 
guished according to the number of pairs of conjugate 
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GRIMM, H. 
Transformation of variates—I/n German 


(Inst. fiir Mikrobiologie und exp Therapie, Jena) 


complex roots. An additional and necessary condition 
that all roots are real for all realisations of the coefficients 
is stated. 

For the treatment of round-off errors of the coefficients 
it is possible to assume that the density function of the 
coefficients is constant in a certain bounded domain 
and is otherwise equal to zero. The density function 
of the roots does not then assume its maximum in the 
interior of the domain in which the density function is 
not identically equal to zero. 

The second part of the present paper is exclusively 
devoted to the cubic equation with real random co- 
efficients and more especially to the cubic equation in 
normal form; without a quadratic term. The in- 
vestigations mentioned above are carried out in full 
detail: examples are given assuming the distribution of 
the coefficients to be normal, rectangular and negative 
exponential. 


(W. Uhlmann) 


2.2 (6.8) 


Biom. Zeit. (1960) 2, 164-182 (159 references, 3 figures) 


This paper is concerned with the transformation of 
discrete random variates for use in analysis of variance, 
regression analysis and other parametric methods based 
on normal distribution theory. 

In the introductory paragraphs the assumptions under- 
lying parametric methods are discussed and tests are 
described for the validity of these assumptions. An 
extensive catalogue of usual transformations is given 
along with methods to show when and how to use them. 
The substantial bibliography is arranged by subjects 
and author’s names. 


(R. Wette) 
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GUMBEL, E. J. (Columbia University, New York) 
Bivariate exponential distributions—Jn English 


ae ee a Ey a Se Rn SR eR ge 


2.5 (2.1) 


J. Amer. Statist. Ass. (1960) 55, 698-707 (7 references, 3 figures) 


A bivariate distribution is not determined by its margins. 
Fréchet in fact has shown that for any given margins 
there exist infinitely many bivariate distributions. The 
author defines three bivariate distributions having 
exponential margins. Their cumulative distribution 
functions are, for x=0, y=0, 


very i= le (e+e *7 95 e 0<d6<1 
ow) Cer le axe *' 7) sas! 
F(x, y) = 1-e *—e ?+exp{-—(x"+ y”)"/"} m=1 
The first two are investigated in detail. None of the 


well-known properties of the bivariate normal distribu- 
tion—such as ellipticity of contours, linearity of regres- 
sion, or homoscedasticity are found to hold. Thus, one 
is cautioned against supposing that the normal case is 
necessarily typical. 

For F, and F, the author obtains the marginal and 
conditional densities, conditional means and variances, 
the correlation coefficient, and the correlation ratio. 


2/467 


HIGUTI, Z. (Institute of Statistical Mathematics, Tokyo) 


The conditional means and standard deviations are 
plotted for a few special parameter values. For F, 
the conditional mean of either variable decreases to 
zero as the other variable increases to infinity; for 
F, it may increase or decrease depending on the sign 
of the correlation, which in turn depends upon the 
sign of a. The correlation coefficient lies in the range 
—0-40<pS0 for F,, and —0:25<ps0-25 for F,. 
Note: there is a typographical slip in (3.6), which should 
read: &(x|y) = 1+40-—ae ”. 


(R. J. Buehler) 


2.9 (-.-) 


Remark on the weak distribution in a locally convex space. I.—Jn French 


Ann. Inst. Statist. Math., Tokyo (1960) 12, 63-67 


In this paper the author gives a necessary and sufficient 
condition for a weak distribution of X to be realisable 
in a space X, if X is a countably normed space which 
is complete and perfect. This condition is analogous 
to the one in the case of the Banach space given by 
Getoor [“‘ On characteristic functions of Banach space 
valued random variables ’’, Pacific J. Math. (1957) 7, 
885-896]. An application to the study of the random 
Schwartz distribution is discussed. 


(M. Huzii) 
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ISHII, G. & HAYAKAWA, R. (Atomic Bomb Casualty Comm., Heccohtmna) 
On the compound binomial distribution—Jn English 


Ne Ra 


2.4 (2.9) 


Ann. Inst. Statist. Math., Tokyo (1960) 12, 69-80 (8 references, 4 tables, 1 figure) 


This paper deals with compound binomial and multi- 
nomial distributions and some of their applications. 
Let a random variable XY have the positive binomial 


distribution (") p\(i—p)""* with parameter p and 
suppose p varies with probability density 
1 
B(a, b) 
(Beta distribution). The overall distribution is called 


the binomial-beta distribution. Further, the bivariate 
binomial beta distribution is defined as follows: 


Pr(x, y) 


B 1 Ny m/s a ema Ny i ee ee f 
[, () er P) ("?) eu py"? °f(p)dp, 


where f(p) is the density function of a Beta distribution. 

In this bivariate case, the authors obtain the con- 
ditional distribution when one of variates is fixed; also 
the distribution of the sum of variates. They show 
particularly that the regression of Y on Y is linear. 


{a= De =p). 3 Ops 
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In the same way, they discuss the k-variate binomial 
beta distribution and the bivariate multinomial-beta 
distribution. 

As applications they refer to comparisons of sex 
ratios in human families, analyses of absenteeism and 
the scoring of a game. 


(Y. Suzuki) 


2301) 


The characterisation of the normai distribution—Jn English 
J. Aust. Math. Soc. (1960) 1, 368-383 (18 references) 


Two lemmas are proved: (i) Let x and y be vectors 
whose elements are sets of m and n mutually independent 
random variables related by y = Ax. Then, if each 
column of A contains two non-zero elements and if 
2<m<n, moments of all orders are finite. If the 
variables x; are identically distributed, the condition 
-may be weakened so that one column of A contains 
- two non-zero terms. 

(ii) If x1, Xz, ---,X, are identically and independently 
distributed and x’Ax is independent of x’Bx, where 
A and B are non-negative definite real symmetric 
matrices and a,,b,, is not zero for at least one j, then all 
moments of x; are finite. 

These lemmas are then used to prove the well-known 
characterisations of the normal by Herschel and by 
Geary and the more recent one by Skitovitch, with 
the aid of cumulant theory. It is shown that Herschel’s 
(spherical symmetry) characterisation requires neither 
the Cramér nor the Marcinkiewicz theorem. The 
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lemmas can also be used to complete other proofs which 
assume the existence, for example, of finite second 
moments. It is then proved that, in the class of sym- 
metrical distributions, the Cochran or the Craig condi- 
tions for independence characterise the normal 
distribution; the method of proof is to show that a 
recursive relation holds between the cumulants of even 
order and that this relation must be the same as holds 
for the normal distribution. 


(H. O. Lancaster) 


MEDGYESSY, P. (Math. Inst., Hungarian Academy of Sciences, Budapest) 
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2.8 (2.5) 


On some unsolved problems in the theory of the decomposition of 


superpositions of distribution functions—Jn English 


Second Hungarian Mathematical Congress (1961) (1 reference) 


We are given the superposition of continuous distribu- 
tion functions 


N 
G&) = is PrP(x—o%3 B,) 


bgt Oe )s i, Ki 1, «-,N) 
where a section of the graph of G(x) and the type of 
F(x—a,; B,) are known; p,>0, «,, B, are, however, 
unknown. 

The problem of determining these unknown para- 
meters, based on the graph section of G(x), arises; 
it is called the “decomposition”? of G(x). In order 
to solve, at least approximately, this problem and also 
more complicated ones, the author has given several 
methods, both in papers, and in his book [*‘ On the 
decomposition of superpositions of distribution functions ”’ 
(1961) Budapest: Akadémiai Kiad6]. 

Some of these methods are based on the finding of 
a special parameter of distribution functions, the so- 
called “formant”, whose definition is as follows: if 
g(x; c) is a distribution function depending on the 
parameter c, then c is called a “‘ formant” (a, m, co) 
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RAMASUBBAN, T. A. (London School of Economics) 


ie On cy = C0) Of @<6=.65,.Cp 1S given 


Le SSS 
if 


0 xm 
a, m are constants, and ¢(x; c) is non-degenerate if 
CoS 6 <0 Or amcse,. 

There are, however, types of distribution functions 
F(x—o,; B,) such that for a superposition formed with 
them, all the mentioned methods fail. An interesting 
type is the following: 

X—- ak = 2 
Flea; By) = | Pash) dy. 
Te Buy 
It is known that the relevant density functions occur 
in physical optics yielding the intensity “‘ distribution” ina 
screen of the light radiating from a lighted gap which is 
parallel to the screen. This superposition could not, how- 
ever, be decomposed by any of the mentioned methods. 

Another superposition which is so far undecomposable 

is that of I distribution functions of order n where 


lim; o(%3 ¢) = U(x; m) = 


c7a 


= oO 


y 
SOc tiate sag 
F(x—o,3 By) = | Z - 


o ~ BL) 
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2.1 (2.4) 


Some distributions arising in the study of generalised mean differences—Jn English 


Biometrika (1960) 47, 469-473 (4 references) 


The author has previously defined the coefficients of 
mean differences of the rth order as the rth power of 
the modulus of the difference weighted by the probability 
distribution functions of the two variable and integrated 
(or summed when they are discrete) over the whole 
permissible ranges. He has previously discussed the 
first order (A,) for the positive and negative binomial, 
Poisson and geometric distributions: [Biometrika (1958) 
45, 549-556; abstracted in this journal No. 1/31, 2.1 
and Biometrika (1959) 46, 223-229; abstracted in this 
journal No. 1/195, 2.1]. He has also expressed A, in 
reasonably compact form in a previous paper. 

In this present paper he examines from a different 
point of view the general problem of expressing A, in 
explicit form. Three distributions, the binomial, Poisson 
and the geometric are treated. 


(Florence N. David) 
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ROGOZIN, B. A. (Moscow) 


On the increase of dispersion of sums of independent random variables—In Russian 


Teor. Veroyat. Primen. (1961) 6, 106-108 (6 references) 
Let €,, €, .-., €, be independent random variables, 
Or!) = sup PrixsGsx tl} 
| OL) = supPrixsei+...+¢,5x+L}, 
s= Y) 1- Quilt 


In this paper it is proved by the author that, where 
C is an absolute constant, and if L>,,,/5, then 


g(L)s 
v5 


(B. A. Rogozin) 
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ROGOZIN, B. A. (Moscow) _ 
On evaluating the concentration function—Jn Russian 
Teor. Veroyat. Primen. (1961) 6, 103-105 (3 references) 


Let €,, ..., €, be independent random variables, 
Q,{1} = sup Pr{xS¢,Sx+], 
OVE sup Prixsé,+...+6,5x+L}, 
a ay, [1-0,(1)]. 
If L=/, then 
cL 
l/s’ 


where C is an absolute constant. 


O(L)s 


(B. A. Rogozin) 
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2.1 (-.-) 


2.1 (-.-) 


GJEDDEBAEK, N. F. (A/S Ferrosan, Copenhagen) 
Contributions to the study of grouped observations. 


IV. Some comments on simple estimates—Jn English 


3.6 (4.3) 


Biometrics (1959) 15, 433-439 (10 references, 1 table) 


This paper is a continuation of a series by the author: 
previous publications were in Skandinavisk Aktuarie- 
tidskrift (1949) 32, 135; (1956) 39, 154 and (1957) 40, 
20. This present paper investigates the means and 
variances of simple estimators of the parameters of a 
normal distribution with grouped observations and where 
a priori 

(i) one parameter is known 

Gi) both parameters are unknown. 


In the first case where o is known, one estimator of 
ft is compared with another and the usual interpretation 
given to the variance of either estimator. Because the 
simple estimators of y and o” are inconsistent, the usual 
concept of efficiency cannot apply. The author proposes 
a derivation of this concept for his particular purpose 
and shows that the inconsistent estimators of w and 
o” have zero asymptotic efficiency. Analogous reasoning 
applies to the situation where yu is known. 
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Where both parameters are unknown it is shown 
that where the grouping is not equidistant the procedure to 
use is that of maximum likelihood: otherwise Sheppard’s 
correction may be used with a simple estimator based 
upon N—1 degrees of freedom. It is concluded that, 
for equidistant grouping with intervals less than twice 
the standard deviation, the simple estimator of the 
variance completed with Sheppard’s correction is as 
efficient practically as the maximum likelihood estimator 
at least for a number of observations less than one 
hundred. 


(W. R. Buckland) 


3.6 (4.3) 


Contribution to the study of grouped observations. V. Three-class 


grouping of normal observations—In English 


Skand. Aktuartidksr. (1959) 42, 194-208'(10 references, 1 table, 21 figures) 


In three previous articles [Skand. Aktuartidskr. (1949) 
32, 135; (1956) 39, 154; (1957) 40, 20; and Biometrics 
(1959) 15, 433-439; abstracted in this present journal 
No. 2/475, 3.6] the author has dealt with the efficiency 
of estimates of the mean and standard deviation obtained 
from grouped normal observations. The present paper 
is mainly confined to a detailed numerical investigation 
of a particular case. 

A sample of N is drawn from a normal population 
with mean € and standard deviation o. The sample is 
grouped into three classes with the fixed boundaries 
—0o<D,<D,<o. Thecorresponding class frequencies 
are n,, 7, andn; = N—n,—n. The normal equivalent 
deviates of n,/N and (a,+n,)/N are Y, and Y, respec- 
tively. It has been shown that, when only the three 
class-frequencies are given, the maximum likelihood 
estimates. of o and € are 


S= (D,—D,)(%—- Y,) 
M = D,- Y,S = D,— Y,S. 
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Let m = (M—6)/o, and ty =./N(M-—6/S. As N 
tends to infinity, m,/N and ty tend to be normally 
distributed with the same standard deviation, t (say). 
A table of t is given for different values of x, and x, 
where x; = (D;—6)/o. 

The exact distributions of m and trl. JN have been 
computed in the case N = 40 for 21 different combina- 
tions of x, and x,. The results are shown graphically. 
It is found that ty/,/N is much closer to normality 
than is m. The agreement between the distribution of 
trl s/ N and the limiting form is judged to be satisfactory 
when each n; has an expectation not less than five. 

Rules for a confidence interval for € are given. They 
involve an iterative procedure for finding the values of 
t to be used in the limits 


M-tt,S/,/N,  M+tt,S/,/N 


(t is the normal deviate corresponding to the desired 
confidence level). 


(B. Matérn) 


MALLOWS, C. L. (University College, London; now Columbia University, N.Y.) 


On the moments of small samples—In English 
Trab. Estadist. (1960) 11, 119-121 (1 reference) 


In this paper the author states that too small a sample 
cannot give close estimates of all the lower moments 
of the population simultaneously. Suppose we have, 
for instance, a sample of size n = 4 from a normal 
population: then the first four raw-moment estimators 
m,, have expectations 0, 1, 0, 3 respectively, but there 
do not exist four real numbers x; such that m/(x) = p/(x) 
forr = 1,2,3,4. The lowest n for which these estimators 
can take these values is n = 6 and there is just one 
possible sample which can have zero probability. For 
n=7 there is a positive probability density at the 
expectation point. 

If we are using k-statistics to estimate the cumulants 
with a normal population, if n=4, then k,, ky, ..., k4 
can simultaneously take their expected values; if n=>6, 
then this applies to k,, ..., ke and if n=9 it applies to 
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ROSSBERG, H. J. (Berlin) 
On the distributions of the differences and ratios of 
Math. Nachr. (1960) 21, 37-79 (3 references) 


Let x,, X2, -.-, X, be independent, one-dimensional 
random variables with the same distribution function 
Prax — f(x). Fhe) related’ order ‘statistics are 
designated by €,, €,, ..., €,. The present paper studies 
the distributions of the differences d,, = €,—€, and 
the ratios g,, = €,/¢, and is especially devoted to limit 
theorems (no) and to questions of statistical depend- 
~ ence or independence. The research is based on results 
of Gnedenko and Smirnoff on the limit distribution of 
Pri{é,<a,x+5b,}. 

It is proved that three types of distributions can occur 
as the limit distribution of Pr{d,,,<a,x} with suitable con- 
stants a,>0. Sufficient conditions are given for the 
occurrence of these three types. 

Analogous limit theorems are proved for the joint 
distribution Pr{d,,24,x, Ank24,y}. The following five 
cases are treated separately: 

Gi) h<k =l<m 
(ii) h<l<k<m 
(ili) h = k<l<m 
(iv) h<k<l=m 
(v) h<k<l<m. 
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k,, .... kg. For the central-moment estimators, the 
figures are »=5 for four moments and n=13 for six 
moments. 

These results are obtained from the polynomial Q(x) 
of degree n whose zeros are the observations x, ..., X,3 
it has coefficients which are elementary symmetric 
functions of the x;. 


(R. Pro) 


3.8 (1.5) 
order statistics—Jn German 


An investigation to discover in which cases the differ- 
ences d,, and d,,, (a<kS/l<m) are independent or 
dependent in the limit is carried out. 

The sufficient conditions mentioned above concerning 
Pr{d,,<a,x} turn out to be sufficient and necessary for 
all linear combinations c,é,+...+¢€, (¢,+0) to have 
non-degenerate limit distributions for fixed k if suitably 
normalised. All possible types of limit distributions 
for Pr{qy,<y,X+6,}, with suitable constants y,>0 and 
6, and continuous F(x), are derived. Again, sufficient 
conditions are given. The author specifies which of the 
conditions is fulfilled if the original distribution F(x) 
is either the normal distribution, the f-distribution or 
any other of the well-known distributions. 

In a final chapter the statistical dependence of the 
order statistics €,, the differences d,,,, the ratios q,, and 
certain functions of these variables are explored. For 
example mtiie Vectors (cj. ems G74) ANd (C1445. +». Sp) 
are independent on the condition that ¢; = a(1</<n). 
For the independence of d,, and €, and likwise of g,, 
and & necessary and sufficient conditions are set up. 
In addition sufficient conditions for the independence 
of all differences and likewise of all ratios are given. 

(W. Uhlmann) 
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ROSSBERG, H. J. (Inst. Appl. Math. & Mech., German Academy of Sciences, Berlin) 


3.8 (1.5) 


Some limiting theorems for order statistics—In German 


Second Hungarian Mathematical Congress (1961) 


Let x, ..., x, be m independent random variables with 
common distribution function Pr{x;<x} = F(x). Let 
us denote the corresponding order statistics by ¢,, ..., 
€,- Under these conditions the investigations concerning 
the limiting distribution of the order statistics €, for 
noo are essentially completed by the papers of 
Gnedenko (1943) and Smirnoff (1949), and this is done 
for “ external numbers’; k = const. or n—k = const., 
and also for “central members” (k = k,— 00, k/n>A, 
0<A<1). In particular the extremal members turned 
out to have the class of possible limiting distributions 
consisting of three proper types of distribution which 
are still dependent on parameters. 

On the basis of the results of the authors noted above 
the limiting distributions of the differences €,—€, (k>h) 
may be investigated. It turns out for k = constant and 
h = constant, if the limiting distribution of &, exists, 
then the same holds for €,—€, and this can be given 
explicitly. Concerning the question as to whether the 
class of the possible limiting distributions of €,—€, is 
exhausted, some statements can be made which, how- 
ever, do not completely solve the problem. Under the 
same condition also the limiting distributions of the 
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vector (¢,—¢,, €,—¢,) (k>h, s>1) exist. In some cases 
limiting independence between &,—€, and €,—€, appears 
if k St in the sense that the limiting distribution consists 
of the product of two distribution functions, for which 
there are given necessary and sufficient conditions. 
Even under the supposition of h = constant, k/n—A 
(0<A<1) and that of the existence of the limiting 
distributions of &, and é,, lim Pr(é,—&,<a,x+b,) may 
be calculated for appropriate sequences {a,} and {b,}. 
Concerning this question the conditions assuring the 
limiting independence of the order statistics are of 
interest. This is known to be the case for the extremal 
members €, and €, given 4 = constant and n—k = 
constant; according to Geffroy [Pub/. Inst. Statist. 
Paris (1959) 8, 36-121, 124-185; abstracted in this 
journal 2/253, 2.6] and a subsequent paper in 1960, 
even in a very strong sense. The limiting independence, 
however, holds between extremal and central members 
too; if A is constant then k>0o. 

Further investigations concern the necessary and 


1 
sufficient conditions for the stability of }\ ¢,¢;; where / 
1 


is a constant and c, different from zero. 
(H. J. Rossberg) 


3.5 (2.7) 


The distribution of regression coefficients in samples from bivariate non-normal 


populations. I. Theoretical investigation—Jn English 


Biometrika (1960) 47, 61-68 (8 references) 


The author assumes a bivariate population which may 
be represented by some terms of a double Edgeworth 
series and considers the distribution of the regression 
coefficient in samples of N drawn from such a popu- 
lation. 

The technique used is that put forward by Gayen. 
The distribution of the regression coefficient is derived 
and the statement is made that this distribution holds 
when the population correlation coefficient is not 
necessarily small; and for any population provided 
the sample size N is large. It will also hold when the 
sample size N is small for populations for which the 
fifth and higher order cumulants are negligible. 

The distribution of the f-statistic used for testing the 
significance of the regression coefficient is also derived 
under the Edgeworth population assumption both for 
the case of known and of unknown marginal variances. 
The author states that he intends to carry out a numerical 


investigation. 


(Florence N. David) 
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3.8 (5.3) 


On the deviation of two-variate empiric distribution functions—In Hungarian 
Magy. Tud. Akad. IIT Oszt. K6zl. (1960) 10, 361-372 (9 references) 


After summarising the most important one-sample and 
two-sample tests which are based on the deviations of 
the distribution functions and mentioning the result of 
Kiefer & Wolfowitz [“‘ On the deviations of the empiric 
distribution function of vector chance variables”’, Trans. 
Amer. Math. Soc. (1958) 87, 173-186] the author gives 
theorems concerning the deviations of the empiric 
distributions of two two-dimensional variates. 

et (6; 4,), (= 1, 2, ..., N be. continuous two- 
dimensional variates independently and _ identically 
distributed with distribution function F(x, y). Similarly 
the independent and identically distributed two-dimen- 
sional variates (€;, 7;), i = 1,2, ..., N have the distribu- 
tion function G(x, y). Let Fy(x, y) and Gy(x, y) denote 
the corresponding empirical distribution functions. 
Let 


DyyQ) ai ei [ Fylx, y)—Gy(x, y)] 


Dyy(y) = cee | Fy(x, y)— Gy(x, y) | 
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Let us choose the value of 7 = y randomly according 
to the probability law F(oo, y). Let P(® and Pl! denote 
the probabilities that Dyy(y) and Dyy(y) do not exceed 
k/N, respectively. The author’s theorems give the exact 


probabilities Py} and Pj\! and also their limiting values. 


(K. Sarkadi) 


3.2 (3.6) 


Note on the non-central t-test based on range—ZJn English 


Second Hungarian Mathematical Congress (1961) 


The non-central ¢ statistic based on a root-mean-square 
estimate of the population standard deviation o was 
dealt with in the original paper of Johnson & Welch 
_ [Biometrika (1940) 31, 362-389]. For purposes of 
practical application, especially in sampling inspection 
by variables, it is convenient to define the non-central 
statistic g using the standard deviation estimate based 
on mean sample range. Tables for the percentage 
points g, computed for the exact sampling distribution 
of mean range have been published recently by the 
' present author [Acta Tech., Czech. Acad. Sciences (1960) ]. 
Since it is known that the distribution of mean range 
in m independent samples each of size n can be approxi- 
mated very well by the distribution cy,/v, where c 
and v have been tabulated as functions of m and n by 
Patnaik [Biometrika (1950) 37, 78-87], it is possible 
to make several comparisons of the directly computed 
values g, with values obtained using the above approxi- 
mation. Three such comparisons have been made: 


(i) based on the original tables of Johnson and 
Welch, 
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(ii) based on the Pearson Type IV curve approxima- 
tion to the non-central f-distribution suggested by 
Merrington & E. S. Pearson [Biometrika (1958) 
45, 484-491; abstracted in this journal 1/39, 3.2], 

(iii) based on the tables of the non-central t-distribution 
published by Resnikoff & G. J. Lieberman (1957) 
Stanford. 


Agreement of the results is discussed. 


(Agnes H. Zaludova) 


CHAPMAN, D. G. & ROBSON, D. S. (Univ. Washington, Seattle and Cornell Univ., N.Y.) 


The analysis of a catch curve—Jn English 


4.3 (2.5) 


Biometrics (1960) 16, 354-368 (10 references, 1 table) 


The age distribution of a random sample, or the “ catch 
curve” from a stationary animal or fish population 
provides information on the annual survival rate s of 
the population. It is commonly assumed that there is 
some age Xo, such that for all ages x =x, the probability 
of selection is the same and the annual survival rate 
is the same. If ages are relabelled so that x) = 0 
the age distribution of the population is geometric, 
that is to say, 
x eS Oe 1D ee 


where a = 1—s is the annual mortality rate. 

The best (uniformly minimum variance unbiased) 
estimate of s is § = x/(1+x-—1/n), X being the mean 
age of the sample. The best estimate of Var § is 
S[S—(t—1)/(n+t—2)] with t = nx. 

While the model giving rise to the geometric age 
distribution may seem artificial, its plausibility is greater 
than might appear since, as the authors show, it can 
be regarded as the result of a random birth and death 
process when suitable assumptions are made. 
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The estimate § is compared to the so-called ‘‘ Jackson ” 
estimate and some modifications thereof. Assumptions 
basic to regression estimates of the survival rate are 
considered. 

Sometimes the sampling gear is selective against older 
ages; also, the assumption of a constant survival rate 
may be valid for only a limited number of age classes. 
In these cases the age distribution may be taken to be 
a truncated geometric over a reduced number of age 
classes. Tables facilitating the estimation of s by 
maximum likelihood are provided. 

A test of the basic model is given for the particular 
case where it is suspected that there may be selection 
against younger ages. Finally, the authors prove that 
for the geometric model there is no unbiased estimate 
of the instantaneous mortality rate i= —logs, but 
that a nearly unbiased estimate is given by 


i* = log (1+x-—1/n)—In x 
—(n—1)(n—2)/n(t+ 1)(n+t-1). 


(H. A. David) 


4,3 (2.3) 


Simplified estimators for the normal distribution when samples are 


singly censored or truncated—ZJn English 


Bull. Int. Statist. Inst. (1960) 37, III 251-269 (11 references, 2 tables) 


In life-testing, dosage response studies, target analyses, 
biological assays and in other related investigations the 
selection and observation of variables is often restricted 
over some portion of the range of possible population 
values; that is to say for some values of the random 
variable involved. Depending on the nature of the 
restriction, samples obtained in this way are designated 
as truncated or censored. Samples in which certain 
population values are entirely excluded from observa- 
tion are described as truncated. Those in which sample 
specimens with measurements falling in certain restricted 
intervals of the random variable may be identified and 
thus counted but not otherwise observed, while the 
remaining sample specimens may be observed without 
restriction, are described as censored. The present 
paper deals with maximum-likelihood estimation of the 
mean and variance of a normally distributed population 
from samples of these types. 
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The estimators obtained in this paper are simpler 
than those previously given for the cases considered here, 
in that the mean and variance of the population are 
estimated by the mean and variance of the sample plus 
simple corrections. These corrections involve only one 
auxiliary function and the practical application of these 
estimators consequently necessitates interpolation in 
only one table, and not in two or more, as was formerly 
necessary. 

For the case of singly-truncated or singly-censored 
samples the auxiliary functions are approximately linear 
if the intervals are not too wide, so that exact inter- 
polation between the number in the table is relatively 
easy in both cases. 


(A. C. Cohen, Jr.) 
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EPSTEIN, B. (Wayne State Univ., Detroit and Stanford Univ., Cal.) 4.3 (2.5) 


Estimation from life test data—Jn English 
Technometrics (1960) 2, 447-454 


In this paper four estimation procedures are discussed 
which can be used to make point or confidence estimates 
from life test data. The first two procedures can be 
applied when the underlying density of life is exponential: 
the third and fourth are non-parametric. 

If life testing is discontinued after a fixed number, ;, 
of items have failed, where items on test may or may 
not be replaced and where the number of items initially 
on test is n=r, the “‘ best ” estimate of mean-life, 0, is 
given by 0,, =T,,,/r where T,., is the accumulated 
life on test until the rth failure occurs. Furthermore 
100(1—«) per cent. confidence intervals for 6 can be 
found using the chi-squared distribution. Similarly, 
expressions are available where the life test is dis- 
continued after a fixed amount of total life T has elapsed. 

Two non-parametric procedures are outlined. In one 
case nm items are tested for time ¢* and the number of 
items that fail in [0, t*] are counted; the items that 
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fail are replaced. In the other case, items are drawn 
at random from some population and tested one at a 
time for a length of time ¢*. Testing continues until 
a preassigned number r have failed. In both cases, 
one-sided confidence interval estimates of both the 
reliability and the parameter 0 are given. 


(P. H. Randolph) 


GRAYBILL, F. A. & MORRISON, R. D. (Oklahoma State University, Stillwater) 4.4 (3.6) 
Sample size for a specified width confidence interval on the variance of a 


normal distribution—IJn English 


Biometrics, (1960) 16, 636-641 (2 references, 1 table) 


This paper is concerned with the determination of the 
sample size needed for the confidence interval estimation 
of a parameter when the length of the interval is specified. 
The authors cite the general problem that confronts the 
experimenter and then proceed to illustrate a method 
for solving this problem, the specific parameter being 
the variance of a normal distribution. 

The determination involves the following assumptions: 
(a) The probability that the confidence interval contains 
the parameter is pre-selected; (6) The probability that 
the width of the confidence interval will be less than a 
specified quantity is greater than or equal to a given 
value. The authors refer to these two probabilities as 
the confidence coefficient and width coefficient 
respectively. 

The problem of determining the sample size in order 
that (a) and (4) are satisfied is outlined in detail. The 
explanation is supplemented with a set of tables which 
is needed in the determination. The theory involved 
is not presented. Assuming an unbiased estimate of 
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the variance is known, the sample size is found by 
entering into the tables the following values: (1) (a) 
and (b) above; (2) An expression to be computed 
involving the specified confidence interval, the chi- 
square variable, the number of observations on which 
the variance estimate was based, and the estimate itself. 
The degrees-of-freedom can then be read, and hence 
the sample size obtained. The tables are presented for 
various values of the significance level and width level. 
An example illustrating their use is also given. 

The authors point out that these tables can also be 
used to solve a similar problem. If it is desired that 
the width of the confidence interval be a particular 
proportion of the variance itself, the required sample 
size can be determined. As usual, the confidence co- 
efficient and width coefficient are chosen. The definition 
of the latter must however be altered to fit the desired 
pre-selected conditions. The method for the deter- 
mination is outlined by the authors and an example is 
given. 


(R. H. Myers) 
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GITTELSOHN, A. M. (State Department of Health, Albany, N.Y.) 


=, ' 


4,3 (2.4) 


A model for the analysis of the distribution of qualitative characters in 


sibships—IJn English 


Biometrics (1960) 16, 534-546 (5 references, 1 table) 


The author considers the problem of estimating 
frequencies of recessive characters through sibship data. 
The basic treatment which assumed fixed sibship size 
and a binomial distribution for the number of recessives 
from the sibship is reviewed. A general model is 
introduced with potential sibship size as a random 
variable. The model is used to describe situations 
where family size is determined by the composition of 
sibs produced as well as natural size variability arising 
out of differences in fertility between couples. 

Each couple is thought of as possessing a potential 
or ideal family size, jointly mediated by fecundity and 
the desire for a particular number of children. The 
realised number of children produced will then depend 
upon the couple’s success in attaining their ideal family 
size. Failure may be through loss of fertility, faulty 
contraceptive technique or through a downward revision 
in their ideal family size as a result of having produced 
one or more children affected by the condition studied. 
In the latter circumstances, such selective termination 
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GITTELSOHN, A. M. 


sibships—ZIn English 


(State Department of Health, Albany, N.Y.) 
A model for the analysis of the distribution of qualitative characters in 


insures that realised family size will depend on family 
composition. 

Two types of ascertainment methods are considered. 
“Complete ascertainment through affected children ” 
whereby all sibships with at least one affected child 
are observable. The probability of ascertaining a 
sibship with at least one affected is therefore taken as 
certainty. The second method, “ incomplete multiple 
ascertainment through affected children ’’ whereby sib- 
ships come under observation with probability dependent 
upon the number of affected children and may be 
ascertained as many times as there are affected children. 

Representing sibship size as a geometric and as a 
Poisson variable, and the number affected and the 
number of independent ascertainments of a given sibship 
as binomials, maximum likelihood estimates of the 
parameters are obtained for the case of no selective 
termination of reproduction. Both methods of ascer- 
tainment are considered. Using the same model, 
maximum likelihood estimates are obtained for the case 


continued 


4.3 (2.4) 
continued 


Biometrics (1960) 16, 534-546 (5 references, 1 table) 


of termination of reproduction at the birth of the first 
affected r (r= 1, 2, ...) children. The model is also 
extended to populations where the termination rule may 
vary among couples. 

Estimates of the probability that a child is affected 
with the condition under study and the probability 
that an affected child is detected are given for both 
ascertainment method and termination rule. The effect 
of various assumptions in the model can be readily 
observed. 


(G. Krause) 
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"KIMBALL, A. W. (Oak Ridge Nat. Lab., Tennessee) Pere ‘ 
Estimation of mortality intensities in animal experiments—In English id Cage 2 Se hepe amie: 
Biometrics (1960) 16, 505-521 (17 Sie 2 tables, 1 figure) Nea ie DRS 
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This paper discusses nonparametric meinods for estimat- deaths are fied in advance fae me exp 


"ing mortality intensities from small samples in controlled _ time intervals are considered to be random vari j 
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OLIVEIRA, J. T. DE (University of Lisbon) 


Bivariate distributions with given marginal distributions and the estimations 


of their parameters—In French 
Second Hungarian Mathematical Congress (1961) 


In some problems concerning bivariate continuous 
random variables, the marginal distributions F(x, +00) 
and F(+ 00, y) are known; in other cases the marginal 
distribution is known with the exception of some para- 
meters, for example the location and dispersion para- 
meters. The estimation of these parameters from the 
samples of X’s or Y’s gives, in general, asymptotically 
the same results. The first case is considered. 

The uniform transformation U = F(X, +00), V = 
F(+ 0, Y) gives a random vector having the distribution 
function K(u, v) with 
m0) = KO.) = 0, K(1, 0) =p and K(u,-1) =z. 
Thus the density k(u, v) may take the form k(u, v) = 
1+h(u, v) where h(u, v) depends on two almost com- 
pletely arbitrary functions. Therefore the estimation 
of the parameters of K(u, v) is non-parametric for the 
marginal distributions. 

The parameter functionals 
P(K) 


= | oes [ates Dpeees Hye OAK (Uy 0, rc AK Uz, 0) 
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SCHMETTERER, L. (University of Hamburg) 
On the theory of unbiased estimators—Jn German 


Second Hungarian Mathematical Congress (1961) 


The problem of the existence and characterisation of 
local minimal unbiased estimates is dealt with; the 
question of their existence under certain suppositions 
is solved by means of the following theorem. 

The norm of a ratio space relative to a closed linear 
subspace of a Banach-space is always taken on if the 
subspace is reflexive. Local minimal estimators can be 
characterised under very general conditions by the 
vanishing of a Gateaux differential. With the aid of 
this characterisation the minimal estimators can be 
determined explicitly in some cases; in order to do 
this a known lemma concerning the representation of 
linear functionals in a ratio space is needed. 


(L. Schmetterer) 


4.8 (4.3) 


can be estimated by 
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method due to Hoeffding, the asymptotical distributions 
of which are multinormal; the same holds for the 
functionals, the first derivatives of which do not vanish 
for K. 

Suppose, for sake of simplicity, that K(u, v) has only 
two parameters 6, and 8,. One may find two functions 
g,(u, v) and g,(u, v) the mean values of which are 
functions of 0, and 6, independent of the parameters 
0, and @, respectively. Therefore, the asymptotic 
distributions of the estimators of 0, and 6, obtained by 
the inversion of the estimators of M(g,) and M(g,) are 
bi-normal. 


(J. T. de Oliveira) 
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SHUMWAY, R. & GURLAND, J. (Iowa State Univ., Ames, U.S.A.) 
Fitting the Poisson binomial distribution—In English 


4.3 (2.5) 


Biometrics (1960) 16, 522-533 (7 references, 4 tables) 


The authors have developed simplified methods for 
computing probabilities and for obtaining maximum 
likelihood estimates of the parameters, a and p, in the 
Poisson binomial distribution 


17.2) i 
Pik) °*)) is | Davee 
t=ot! k 

A table is provided which greatly simplifies computation. 

The density function can be written to contain explicitly 
the expression for the kth factorial moment, Urpy, Of 
the variable nX, where n is a known positive integer 
and X is Poisson with mean 4, 2 = ag", and gq = 1p. 
Exploiting this device, a recurrence relation for the 
computation of Pr(k) is 


Pr(k+1) =? rae Pry 


q (k+1 
where 


— Mik+1] 
Di hi em narae = 
Hex) 
Given values for pr,;, the computation is obviously 


simple. The likelihood equations are nép = k and 


L(p) = 
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SPECHT, W. (Universitat Erlangen) 


On the theory of elementary mean values—I/n German 


Math. Zeit. (1960) 74, 91-98 (1 reference) 


Let P, be the set of all n-dimensional, real vectors 
X = (X,, X2, --- X,) with positive components x, and 
n>1 and let p = (Pp, Po, ---, P,) be a distribution of 
positive weights p; with py+p.+...t+p,=1. The 
elementary mean value of the vector x is defined by 


M,(p, x) = (2 p,x,)'" 
for each real number s+0 and 
M,(p, x) = exp (2 p, log x,). 
The limits for s to —o and s to + © are 
M_,(p,x)=minx, and M,.(p, x) = max x,. 


Let O, , = Mt/Ms, then it is wellknown that Good Ghee 3 eal 
for all s and ¢ with s<t. 

QO40,-.1 =maxx,/minx, is called the relative 
range f(x) of the vector x. For all vectors x of P, 
with 1<f(x)<B and all distributions p of weights, an 
upper bound for Q,, ,(P, x) is given for —co0o<s<t<+0. 
This upper bound depends on f, s and B only and is 
proved to be optimal. 
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where the a, are observed frequencies and N is the sample 
size. 

The first estimate of p is made by some simple method; 
for example, sample moments or frequencies—and it is 
then calculated accurately by application of Newton’s 
method, using L(p) and a modification L’(p). Obtaining 
G once p is known is straightforward. 

As in computation of the probabilities, the procedure 
is not difficult if the pj,; and, in this case, the q,,; are 
known. The authors have computed the values for 
Ptxy and qy for different values of A, 2 = 0-10 (0-02) 
1:10 and fork = 0(1)9. Their methods for computing 
Px) and gj, are developed and exemplified in one 
section of the paper. 

Two examples are given showing the estimation of 
Poisson binomial probabilities and parameters. Some 
comments are made on expedient means of interpolation 
in the table pry, 943; there is also a brief discussion 
of the problem which arises when the zero class is 
missing, that is to say when the distribution is truncated 
on the left. 


(H. R. Baird) 


4.8 (4.0) 


Finally this research is extended to functions f(x), 
which are defined, positive and integrable in the interval 
[0, 1]. Let p(x)=0 be defined in the interval [0, 1] 
with the integral of p(x) between the limits 0 and 1 equal 
to 1. The integral average of f(x) is defined by 


M,(p, f) = | [Peon |" for s+0 
0 


1 
M(p, f) = exp | (low f(s) sa} 
0 


The definition of QO, ,(p, f) is analogous to the definition 
given above. For all functions f(x) integrable in [0, 1], 
for which 

O0<as<f(x)SA, forOSx<1 
the same upper bound of Q, ,(p, f) as above is valid 
for —o<s<t<+o if Bis replaced by A/a. 


(W. Uhlmann) 


BAHADUR, R. R. 
Stochastic comparison of tests—In English 


(Indian Statistical Institute, Calcutta) 


III IIE ETE 


5.1 (1.5) 


Ann. Math. Statist. (1960) 31, 276-295 (13 references) 


Let Ti” and T,?) be two sequences of real-valued test 
statistics for testing the hypothesis that @ belongs to a 
sub-set Q) of Q, where Q indexes a set of probability 
distributions. 

The method of stochastic comparison outlined below 
gives a measure of the relative efficiency of the two 
sequences at each parameter point in Q-Q). It is 
assumed that each sequence has a limiting distribution 
function independent of @ in Qo, and fer 0 in Q—Q,, 
T;?/,/n converges in probability (Ps) to a positive 
number 6,0). Let Ko = —2log[1—F(T,)], and let 
C,(@) be the limit in probability (P,) of K°?/n. Under 
certain assumptions made on F;, C,(@) = a,[b,(@)]? for 
8 in Q—Q), where a;>0 depends on the distribution 
function F;. The ratio $,,(0) = C,(0)/C,(8) is proposed 
as a measure of efficiency of the first sequence relative 
to the second sequence. The author gives some justifi- 
cation for this as a measure of efficiency on intuitive 
grounds and shows that in cases where the Pitman 
theory of efficiency is applicable, the Pitman efficiency 
is obtained as a limit of ¢,,(0) as 6 approaches some 
point in Qo. 
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BENNETT, B. M. & HSU, P. 


(University of Washington, Seattle) ‘ 
On the power function of the exact test for the 2x 2 contingency table—Jn English 


Furthermore, it is shown that, if for all n the power 
of the test 7’) never exceeds the power of T;?) at 0, 
$12(0)S1; and conversely, if ¢,,(0)<1, then the above 
inequality must hold for each n. In addition, it is 
shown that under certain conditions C(@) can be 
represented as the limit as n tends to infinity of 1/n 
times the expected ratio of the power of T,, to the size 
of T, where the size is chosen at random according to 
a certain fixed distribution. 

Several examples are given in which the relative 
efficiencies are computed for standard tests of the 
following hypotheses: the mean of a normal distribution 
is zero, the variance of a normal distribution is one, 
two distribution functions are equal, several distribution 
functions are equal, and that a bivariate normal distribu- 
tion has correlation coefficient zero. 


(D. R. Truax) 


5.2 (6.9) 


Biometrika (1960) 47, 393-397 (9 references, 10 figures) 


The exact tests of significance in a 2 by 2 table for the 
equality of two population proportions is discussed: 
this hypothesis was set out by Fisher in 1935. The 
conditional power function under the hypothesis that 
these proportions are not equal and also the overall 
power functions are discussed. 

Ten power diagrams for different sample sizes have, 
been calculated: (5, 5), (10, 10), (10, 5), (15, 15), (15, 10) 
(15, 5), (20, 20), (20, 15), (20, 10) and (20, 5). A com- 
parison with Sillitto’s approximation is made. 


(Florence N. David) 
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BHATE, D. H. (University Coll., London and Delhi School of Economics) 


000,00, CCOOECCSSSSS FF ee ee 


Approximation to the distribution of the sample size for sequential tests. 


Il. Tests of composite hypotheses—Jn English 
Biometrika (1960) 47, 190-193 (16 references) 


The author extends the procedure of his previous paper 
[“ Approximation to the distribution of the sample 
size for sequential tests. I. Tests of simple hypotheses ”, 
Biometrika (1959) 46, 130-139: abstracted in this journal 
No. 1/223, 5.7] to the case of composite hypotheses. 

Assume two normal populations with means p, and 
Mz and variances o, and a, respectively with A = 07/03. 
The hypothesis Hy is 4 = 1) and H, that A =A,. 
Given that at a certain stage there are samples, each 
of m measurements, available from each population, 
the author suggests the calculation of the quantity F,, 
which is the ratio of the two sample sums-of-squares 
and proposes the rules: 


(i) continue to the (m+1)st stage if F< F,, <F,, 
(ii) accept H, if F,, 2 F,,, 
(iii) accept H, if F< F,,. 


2/497 


BROSS, I. D. J. 
Statistical criticism—ZJn English 


Cancer (1960) 13, 394-400 (13 references, 2 tables) 


This is a discussion of the role and responsibility of the 
statistical critic. The author feels that an essential 
feature of statistical criticism is the implicit or explicit 
presentation of a hypothesis counter to the one being 
criticised. 

The following tentative criterion for criticism is 
suggested. ‘‘ The critic has the responsibility for showing 
that his counter hypothesis is tenable. In so doing, he 
operates under the same ground rules as a proponent.” 

Although this criterion would require proponent and 
critic to operate under the same rules in establishing 
their respective hypotheses, the burden of proof would 
rest on the proponent. The proponent would have 
to rule out all tenable hypotheses except his own, the 
critic need only show his hypothesis is tenable. The 
recommended minimum requirement for the establish- 
ment of tenability by the critic is that the effects predicted 
from his hypothesis be in line with the actual data both 
in direction and order of magnitude. The amount of 
additional argument required would depend on the 
hypothesis. 
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(Roswell Park Memorial Inst., Buffalo, N.Y.) 


F,, and F,* are the solutions, respectively, of the 
equations 


B/l—«, 
(1—P) fo = (Ag fay PY {1 + FnlAo) (A+ Flas} * 


where a and f have the meanings generally given in 
sequential analysis. At the end of this note some 
examples are given. 


(Florence N. David) 


5.0 (11.9) 


Discussion and examples of several types of criticism 
constitute the bulk of the paper. The topics treated are 


(1) ‘‘ Hit and run ”’ criticism (enumeration of real or 
fancied flaws in a study), 

(2) Dogmatic criticism (rejection, without further 
analysis, of data which do not satisfy certain 
statistical desiderata), 

(3) Speculative criticism (introduction of speculation 
into the critic’s counter hypotheses or conclusions 
without establishing tenability), 

(4) ‘ Tubular ”’ criticism (characterised by an inability 
to “see”? evidence unfavourable to the counter 
hypothesis). 


The examples are taken from the current controversy 
over the relationship between smoking and lung cancer. 


(Polly Feigl) 


5.7 (5.2) 
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EPSTEIN, B. (Wayne State Univ., Detroit and Stanford Univ., Cal.) 
Statistical life test acceptance procedures—In English 


5.2. (8.6) 


Technometrics (1960) 2, 435-446 (4 references, 3 tables, 2 figures) 


In this paper the author gives methods for testing the 
hypothesis Hy:0 = 0, against the alternative H POs 
with Type I and Type II probabilities of « and B, 
respectively, for a random variable which has an 
exponential distribution. Advantage is taken of the 
time-ordered nature of life test data to substantially 
shorten the time required to reach a decision. In the 
“ censored tests’, items are tested until r fail, where 
ris obtained from the operating characteristic or power 
curves. The average life of all items up to the time 
of the rth failure is computed and composed with a 
decision rule. 

In the “ time-truncated life-tests’’ a predetermined 
time Ty is given. If items are tested and rp fail, Hy 
is rejected. If less than r fail up to time To, Ho is 
accepted. 

Finally a sequential procedure is outlined: several 
examples are given. . 


(P. H. Randolph) 
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JOHNSON, N. L. (University College, London) 


On the choice of a sequential test procedure—In English 


Symp. Quantative Methods in Pharmacology (1960) 


In the construction of a standard sequential probability 
ratio test procedure we need to define two pivotal 
hypotheses Hj, H, and two approximate probabilities 
of error %, %,, such that 
Prireject H;| H,;}=a,, j = 1, 2. 
It is known that the sequential probability ratio test 
is then optimal in the sense that, among all procedures 
satisfying this expression, the test has nearly the smallest 
expected size of sample whether either Ho or H, is true. 
The conditions given above are the same as those 
often used ‘in determining the size of sample to be 
taken in a classical fixed sample procedure. In such 
cases the hypotheses Hy and H, together with the 
probabilities %) and «, are chosen in the light of the 
required sensitivity of the procedure. For example, 
H, may be a null hypothesis, and a a level of signifi- 
cance corresponding to a desirable level of quality; 
with H, representing a change in quality of such an 
amount that a high probability (1—«,) of detecting it 
is needed. If H, and H,, determined in this way, are 
used to define a sequential probability ratio test the 
full advantage of the sequential method may be lost 
unless it is especially desirable to average minimum 
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sample size when either Hy or H, is true. When a | 
hypothesis H, different from both Hy and H,, is true | 
there is no guarantee that this test procedure defined 
by Ho, H,, %, %,, will even have a lower average sample | 
size than a fixed sample procedure defined in the sameway. | 

A sequential probability ratio test satisfying the above | 
conditions need not necessarily be based on Hy and | 
H, as pivotal hypotheses. It is suggested that con- | 
sideration should be given to the introduction of 
additional hypotheses Hj, Hj, to be used as pivotal | 
hypotheses. Hg and H; would be chosen so that the 
saving in average sample size should be as effective as | 
possible. The associated approximate probabilities of i 
error &, a; (where Pr{reject Hj| HjXa;}; j=1, 2) | 
are then determined by the conditions given above. | 

In this paper detailed consideration is given to the | 
case when the sequential probability ratio test is regarded _ 
as a test of the null hypothesis Hj, and where we take | 
H4 = Ho. Incidentally to the main topics, some results — 
are obtained on the comparison of operating character- | 
istics or power functions of comparable sequential | 
probability ratio test and ancillary tables are presented. _ 


(N. L. Johnson) 


HEMELRIJK, J. (Mathematical Centre, Amsterdam) 


OO ee ee ee 


5.3 (5.0) 


Experimental comparison of Student’s and Wilcoxon’s two sample tests—In English 
Symp. Quantative Methods in Pharmacology (1960) 


The choice between the many available statistical 
methods is often difficult to make in practical situations. 
It is often, too often, based on personal preference 
without sufficient regard for the merits of the different 
methods available. A lack of knowledge about the 
small-sample properties of many  distribution-free 
methods tends to make their position uncertain in 
comparison with classical methods based on specified 
suppositions about the sampled populations; the 
properties of the latter are known, if these suppositions 
are true. 

In a small experiment, using Quenouille’s table 
(“Tables of random observations from standard 
distributions”’, Biometrika (1959) 46, 178-204; abstracted 
in this journal No. 1/325, 11.1] fifty independent pairs 
of samples of size ten were used to compare student’s 
and Wilcoxon’s two-sample tests. The samples were 
first taken from a standard normal distribution 
(Quenouille’s x,), then one of them was shifted to the 
right over distances of 0-81; 1-05 and 1-53. These 
shifts correspond to a power (for Student’s test) of 
0-40, 0-60 and 0-90 respectively. In all these positions 
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HEMELRIJK, J. (Mathematical Centre, Amsterdam) 


Experimental comparison of Student’s and Wilcoxon’s two sample tests—Jn English 
Symp. Quantative Methods in Pharmacology (1960) 


difference between the two tests was, in this case, at 
least as large as the difference in reaction to the shifts 
mentioned in the foregoing paragraph. 

Finally all observations, transformed to observations 
from a rather skew exponential distribution (Quenouille’s 
x4), were again tested in the same way. The shifts were, 
in this case, replaced by multiplication of the observations 
of one of the samples from the starting point of the 
distribution in order to keep this point invariant. The 
factors were: 1 (identical distributions); 2; 3 and 4:5, 
resulting in power function values of roughly 0-025, 
the level of significance; 0:25; 0-45 and 0-70 respectively. 
The two tests were now approximately equally powerful. 

All these results are dependent in a statistical sense, 
being based on the same 1000 observations. The 
investigation was too small, although laborious, to 
warrant strong conclusions: it is only the beginning 

’of more extensive research using more powerful means. 


(J. Hemelrijk) 
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both tests were applied one-sidedly with level of signifi- 
cance 0:025. As was to be expected, Student’s test 
being uniformly most powerful in this situation, the 
power of Wilcoxon’s test was the smaller one. 

The results indicated that in this case one would not 
be far off the mark in stating that, roughly, the loss 
of power was about ten per cent. of the power of | 
Student’s test as long as the latter is not too close to 
unity. This also holds under the hypothesis tested: 
the true level of significance for Student’s test is 0-025 
and 0-022 for Wilcoxon’s test in this example owing 
to its discrete character. In about 90 or more per cent. 
of the cases, however, both tests led to the same con- 
clusion about the hypothesis tested. 

A second point of investigation was the reaction of | 
the tests to an observational error in one of the original 
observations. For each pair of samples an observation 
was chosen at random from the sample with the larger 
median and this observation was shifted over a distance 
two to the right. In accordance with expectation 
Student’s test proved to be far more sensitive to this 
mild kind of slippage than Wilcoxon’s test. The 


] continued 


5.3 (5.0) 
continued 


ay 


JUVANCZ, I. (Math. Inst., Hungarian Academy of Sciences, Budapest) 


5.0 (11.0) 


Contra-indication of non-parametric tests in medical experimentation—Jn English 
Symp. Quantative Methods in Pharmacology (1960) 


According to the author’s personal experience from the 
design and/or analysis of over 700 experiments, most 
of them medical, the statisticians engaged in medical 
experiments show a reluctance to use non-parametric 
methods, though numerous theorists enthusiastically 
advocate them. A statistical survey is made of papers 
from some medical, biometrical and _ statistical 
periodicals, published in various countries. The results 
show the same tendency as the author has noted at home. 

The reluctance of those “‘ qualified in medicine or, 
if not so qualified, at least well soaked in it” is caused 
by some characteristic features of medical experiment- 
ation, well known to them. So, for example, in medical 
experiments the interest is mainly centred on the 
estimation of mean-values and their fiducial limits, and 
on material significance. There is great variation 
between and within individuals; the measurements are 
poor; the number of observations is small. Experi- 
ments are sometimes dangerous. Numerous factors 
disregarded in the statistical analysis influence the final 
decision. Experimental results are retested several times 
by others before acceptance and are often intended 
for individual use. 
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These characteristics compel the research workers to 
use statistical methods which yield results convertible 
to the original unit of measure. Order statistics and 
the ?-tests are not of this kind. To minimise the loss 
of information, non-parametric methods are too wasteful : 
consequently non-parametric tests are economically 
inefficient both qualitatively and quantitatively. Often 
non-parametric tests are insensitive to changes regarded 
as medically substantial and too sensitive to insubstantial 
ones. An example of this is the fact that y?-tests are 
insensitive to consistent though small increases and 
sensitive to great though inconsistent ones; this is not 
the case with the f-test. It is often impossible to decide 
what caused the statistically significant change indicated 
by a non-parametric test (mean, type of distribution, 
etc.). On the other hand, tests of the mean are very 
“robust” against non-normality. Only a small field of 
application is left for non-parametric methods, cases 
where the experimental techniques are too weak. 
Consequently, the use of non-parametric methods 
should be avoided, if possible by the invention of better 
experimental methods. 


(I. Juvanez) 


5.3 (6.7) 


On samples with two general populations and tests of hypotheses connected 


with them—ZJn Ukrainian 


Dopoy. Acad. Nauk Ukrain, SSR (1960) 150-155 (4 references) 


The samples X14, ..-; X1,, and X,2, -.., X2,, of two 
populations obeying the laws of distribution Pr(x: «,, 0,) 
Pr(x: «>, 0) respectively are considered by the author 
in this paper. 

The following hypotheses are examined: 

(i) H,—a, = a; 0, and 0, are known 

(ii) H,—a, = a,; if 0, = 0, for Pareto’s law. 
The tests deduced are analogous to the tests obtained 
by Epstein & Tsao [Ann. Math. Statist. (1953) 24, 458- 
466] for exponential distributions. 


(B. V. Gnedenko) 
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Approximation tests for m rankings—In English 
Biometrika (1960) 47, 476-480 (7 references) 


One or two sets of rankings of n objects are given. 
The sets of size m, and m, are regarded as a random 
sample from a population of rankings. Approximate 
tests for hypotheses regarding these rankings are derived. 

The first hypothesis, for a single set of markings, is 
that the mean ranks of a certain subset of the n are 
identical. The second and third hypotheses, regarding 
two sets of rankings are: 


(i) the n mean ranks are identical in the two popula- 
tions, 

(il) the coefficients of correlation are identical in the 
two populations. 


The author uses techniques which have been previously 
put forward by Stuart [Biometrika (1951) 38, 33-42]. 
He chooses statistical criteria on the basis of their 
reasonableness, obtains the first tvo moments of these, 
and then uses as approximations to their distributions 
the simplest tabulated statistical distributions: that is to 
say, the normal distribution for criteria within the 
range —« to +a, and the chi-squared distribution for 
criteria with range zero to +0. 


(Florence N. David) 
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MAURICE, Rita (University College, London) 


SESS S'S ee 
LINHART, H. (Nat. Inst. Personnel Res., S. African Coun. Sci. Industr. Res.) 


5.6 (6.5) 


5.3 (1.8) 


A different loss for the choice between two populations—Jn English 


J. R. Statist. Soc. B (1959) 21, 203-213 


A decision is to be made between two populations 
characterised by the values 0,, 0, of an unknown para- 
meter, on the basis of a sample of size n from each 
population. 

It is argued that, for example in an industrial process 
producing a series of N objects, the first 2” of which are 
to be given one of two treatments and the production 
of the rest based on the decision between these, but, 
where the test pieces may also be sold, the expected 
loss-function should be 

[((N—2n)p+n]@ 
where p is the chance of a wrong decision and 
0 = |0,-8, |. 

Reasons are given for this function to apply to a 
range of cases where WN has the wider interpretation of 
“extent of use”. The differences between this and the 
related loss-function of Healy et a/.and Box are discussed. 

It is pointed out that the minimax principle cannot 
be applied (to the problem of obtaining optimal n) 
without further specification since the most unfavourable 
value of @ gives “infinite” loss; an upper bound, t 
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to 0 is therefore supposed known. The set-up where 
the parameter is the mean of a normal population with 
known variance o” is then treated in detail and it is 
found that there is no simple expression for optimal n; 
even the form of the maximigatation equation depending 
on t, but numerical solutions are tabled against N and 
T = t/o,/2 and it is noted that unless T is as small 
as one or two, No, the minimax n does not vary much 
with 7. 

The author next considers sequential tests, &(n), in 
the minimax case, being compared with no. The 
sequential scheme is much superior on this count but 
it is pointed out that, owing to approximations, the 
comparison is only valid for large N and even then, it 
may not be experimentally possible to conduct sequential 
trials. 

Finally, the expected loss for both types of scheme is 
plotted against 6/o for N = 100, 2500 and it is seen 
that the sequential plan is very nearly uniformly better. 


(D. E. Barton) 


__ PEARSON, E. S. (University College, London) 


f 


Editorial note of Bennett and Hsu’s paper “On the power function of the 
exact test for the 2 x2 contingency table ”—Jn English 


0 IDEOEEOSSS S.-C a ee 


5.2 (6.9) 


Biometrika (1960) 47, 397-398; Original paper abstracted No. 2/469, 5.2 


Bennett & Hsu [Biometrika (1960) 47, 393-397] gave 
a comparison of Sillito’s approximation with their own 
computed values. The Editor of Biometrika here 
discusses Sillito’s and Patnaik’s approximations, pointing 
out that Patnaik’s power values are only “exact” for 
the particular significance points obtained by assuming 
that the hypergeometric distribution under the null 
hypothesis could be represented by a norma! distribution. 

Attention is drawn to the fact that the Bennett & Hsu 
power diagrams bring out clearly the difficulty of dis- 
tinguishing between two widely different probabilities 
when the samples are small. 


(Florence N. David) 
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5.6 (5.3) 
Bell Telephone Labs.) 


Contributions to the theory of rank order statistics. The two-sample 


censored case—IJn English 


Ann. Math. Statist (1960) 31, 415-426 (4 references) 


Rank order theory is developed for the two-sample 
problem in which censoring of the observations has 
Occurred: . More explicitly, let X,,:.<:, Xi. Yas +05 Yo 
be independent random variables where the X’s have 
distribution function F and the Y’s have distribution 
function G. Tests of the hypothesis F = G are con- 
sidered. A rank order z is defined as z = (z,, ..-, Zm4n) 
where z; is zero or one according as the ith smallest 
in the combined sample is an XY or a Y. Rank orders 
may be described as paths of horizontal and vertical 
steps of unit length on the integer lattice from (0, 0) 
to (m, n) where a zero is a horizontal step, and a one 
is a vertical step. 

In this paper a censoring scheme may be described 
as a set S of lattice points not containing (0, 0) such 
that each path from (0, 0) to (m, n) meets S. Experiment- 
ation is started at (0, 0) and continued until a point 
in S is reached. Under some censoring schemes the 
total sample size is a random variable. In the case 
where experimeritation is continued until m* of the 
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X’s have been observed the probability distribution of 
the total sample size is given for the case F#G and for 
the case F=G. For the general type of censoring 
scheme described above the conditional distribution of 
the rank orders given the total sample size is presented 
in terms of F and G. 

If F(x) = H(x, 0) and G(x) = H(x,. 0), where 0 is a 
real parameter, a locally most powerful rank order test 
can be constructed by choosing a test which maximises 
the derivative of the likelihood ratio at 0 =0. For 
the general censoring scheme the derivative of the likeli- 
hood ratio at 0 = 0 is given when H(x, 0) is the normal 
distribution function with mean 0, when H(x, 0) = 1 
—[1—J(x)]'*®, and when H(x, 0) = (1—6)J(x) + 0J?(x) 
where J is a fixed distribution function. In the un- 
censored case these formulas reduce to known results. 

In case F’ and G’ have a monotone likelihood ratio 
a partial ordering of the likelihood ratio of rank orders 


is given in the censored case. 


(D. R. Truax) 


REIMANN, J. & VINCZE, I. (Math. Inst., Hungarian Academy of Sciences, Budapest) 


5.3 (5.6) 


On the comparison of two samples with slightly different sizes—In English 
Publ. Math. Inst. Hung. Acad. Sci. (1960) 5, 293-309 (8 references) 


Let F,(x) and G,,(x) denote the empirical distribution 
functions of two samples of size n and m (m>n), taken 
independently from populations with continuous dis- 
tribution functions F(x) and G(x). The authors deter- 
mine the distributions of the following statistics 


Bim = Max [n F,(x)—mG,,(x)], 
(x) 
and 
Bum = nes | nF,(x)—mG,,(x) +(m—n)/2 | —(m—n)/2 


and the limiting distributions for the case when the 
sizes of the two samples only “ slightly differ’, that is 
if noo and (m—n)*/(m+n)>4c”, where c>0 is a 
constant. It is proved furthermore that in this case 
the tests based on the statistics B,,, and B,,, are asymp- 
totically consistent against all continuous alternatives 
F(x) G(x) and F(x)>G(x) respectively. 

The following test is suggested in such cases. (In 
the finite case to compensate for the lack of closer 
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(Inst. Theor. Biol., Leyden University) 


investigations c<1 is recommended.) Let the com- 
bination of the above-mentioned samples be denoted 
by (,<(,<...<¢,,,, arranged in order of magnitude. 
Let S; = nF,(C;+0)—mG,,(6,+0) and R,*, and T,7, 
respectively denote the first and last of the indices i 
for which the sum S; is maximal and R,,, the first index 
for which | S;+(m—n)/2 | —(m—n)/2 is maximal. The 
authors determine the joint distributions and limiting 
distributions of the pairs of statistics (By,, Ri); 
(Bea. ey and (Bains Rom): 

As special cases for n = m, that is, where c = 0, the 
distributions of Gnedenko & Korolyuk [Dokl. Akad. 
Nauk SSSR (1951) 80, 525-528] and Vincze [Pub/. 
Math. Inst. Hung. Acad. Sci. (1957) 2, 183-209] are 
obtained. The authors’ work is connected with the 
method used by Hodges in his paper [“‘ The significance 
of probability of the Smirnoff two-sample test”, Ark. 
Mat. (1958) 3, 469-486]. 


(K. Sarkadi) 


On the robustness of Wilcoxon’s two sample test—In English 
Symp. Quantative Methods in Pharmacology (1960) (4 tables, 4 graphs) 


The main object of this paper is to compare the asymp- 
totic (for large sample sizes) behaviour of the probability 
of type I error for Wilcoxon’s and Student’s two-sample 
tests in the following situation: a sample of size m is 
drawn from a N(u;, 07) distribution, and a sample of 
size n is drawn from a N(u>, 03) distribution. Then, at 
significance level «, the null hypothesis Hy; HM, = My is 
tested against H; .>j,, by means of the appropriate 
one-sided Wilcoxon’s and Student’s two-sample tests. 
The probability of rejecting Hy when in fact it is true, 
depends on o,/c, and is asymptotically equal to Q(1,/W) 
for Wilcoxon’s test, and to Q(/,/S) for Student’s test 
[for two-sided tests the corresponding probabilities are 
20(A,/W) and 20(/,/S)]. Here « = QO(A,) = Pr(v>2,), 
v being a (0, 1) variate. Put q = (oj—93). (o7+063)', 
then 
W? = —34+122-'(m+n)~1[m.arecos(3(1+4)*) 
+n.arccos(4(1—q)*)] = 1+2.3%.27*.q.(n—m) 
(ude oe on bg ...; 
S? =[m+n—q.(m—n)].[mt+n+q.(m—n)]"* =1 


+24 .(n—m)(n+m)~!+2q?.(n—m)(n+m) 7? +... 


A number of conclusions may now be drawn. 
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One result is for equa! sample sizes, m = n: whereas 
QO(/,/S) is then known to equal « for any value of q 
(that is, for any ratio o,/0,), Q(A,/W) is found to vary 
with qg and in fact to be larger than Q(A,/S) for each 
g #0. On the other hand, if the size of one sample 
is k(=2) times the size of the other sample, Wilcoxon’s 
O(A,/W) deviates from a less than does Student’s 
O(/,/S), again for each g £0. Finally, if the size of 
one sample is k(1<k<2) times the size of the other 
sample, Wilcoxon’s test behaves better only on part of 
the g-range. In the immediate neighbourhood of g = 0, 
that is, of ¢, =>, the ratio O(A,/W)/Q(A,/S)¥3* .n* 
~0-55. Four tables are given for numerical values of 
O(4,/W) and Q(,/S) for n =m, n= 3m, n= 2m, 
n = 3m, (« = 0-05) together with four graphs. 

Summarising, in the two-sample problem in the case 
of two normal distributions with unequal variances, the 
asymptotic behaviour of the probability of type I error 
for Wilcoxon’s test is not uniformly better than for 
Student’s test; but Wilcoxon’s test behaves rather worse 
only for rather extreme values of o,/¢,, combined with 
certain values of the ratio m/n between $ and 2. 


(H. R. van der Vaart) 


5.3 (2.3). 


VINCZE, I. (Math. Inst., Hungarian Academy of Sciences, Budapest) 


On a special type of randomised tests—In English 


Second Hungarian Mathematical Congress (1961) 


Let thin = tnm(Sis G25 +++> ns Nis +++ Nm) be a Statistic used 
for deciding whether the null-hypothesis F(x) = G(x) is 
true. Here ¢,,..., ¢, and m,, ..., 7,, are samples taken 
from populations with the continuous distributions 
F(x) and G(x) respectively. Let us consider the case 
when ¢,,, is not asymptotically consistent against all 
continuous alternatives. 

The author shows some examples in which f¢,,, may 
be imbedded in a parametric set of statistics 


Tam(D) Se cn (a se\rlg eh Ue ee, Nm a)d eQ, Tam(9o) = baa 


such that (i) in case the null-hypothesis holds the 
distribution of 1,,,,(9) does not depend on 9, (ii) choosing 
9 randomly, the randomised test obtained will be with 
probability one (relative to the measure defined on Q) 
asymptotically consistent against all continuous alterna- 
tives; that is to say, for any given G(x) = F(x) with 
probability zero will be chosen so that the test based on 
Tam() Will not be consistent. 


(I. Vineze) 
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5.2 (2.6) 


A test of fit based on the largest sample spacing—Jn English 


J. Soc. Indust. Appl. Math. (1960) 8, 295-299 (3 


A test-of-fit is presented based on V, the largest of the 
n+1 sample spacings obtained when the probability 
transformation is made on v independent, identically 
distributed random variables (x,, x2, ..., x,) from a 
continuous unknown distribution. The completely 
specified function F(x) used for the transformation is 
the hypothesised cumulative distribution function. 

The ordered values on the interval [0, 1] are denoted 
Y,5Y,S..-SY,. Yo denotes 0 and Y,,, denotes 1. 
T, denotes Y;,— Y;-, fori=1, ....n+1. 7), T), ---, 
T,,, are the n+1 sample spacings of which V is the 
largest. 

When the hypothesis is true the distribution of V is: 


prved=¥ @ : (11 -ajy" 


J 
Ais atl Lh 
O0<j<min (je n+?) 
and 


n+1 


n> oo 


1 oe 
lim Pr & eens | = exp(—e ) 
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references) 


The test rejects when V2d,(«). The critical value 
d,(«) depends on the sample size, and the desired level 
of significance a. 

Properties of the test were investigated by the author 
and it was shown to be admissible and consistent. 


(D. T. Searls) 


5.3 (5.6) 
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ARMITAGE, P. & BARTSCH, G. E. (London School of Hygiene and Tropical Medicine) 
The detection of host variability in a dilution series with single 


observations—IJn English 


6.8 (5.6) 


Biometrics (1960) 16, 582-592 (13 references, 4 tables, 1 figure) 


Current techniques are investigated for testing the hypo- 
thesis of equal host susceptibility to infection from a 
single particle of an inoculum. Attention is restricted 
in the paper to the case of single observations that 
form a two-fold dilution series. The Stevens range 
statistic R, the number of dilutions between (and 
including) the first at which not all observations are 
positive, that is to say the infection observed, and the 
last at which not all are negative is compared with a 
new Statistic J, the number of positive results at dilutions 
beyond that at which the first negative occurs. It is 
shown that the asymptotically efficient statistics ¢ and 
’ due to Armitage are more closely related to J than 
to R. The distribution of J under the null hypothesis 
of equal susceptibility is computed for a two-fold dilution 
series. The power function for the J test is calculated 
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(Central Planning Bureau, The Hague, Netherlands) 


with host variability a binomial (half of the hosts with 
susceptibility level po, half with level 1—p,) variable. 
The J-statistic appears to have a steeper power curve 
for small departures from the null hypothesis, while 
for large values of the departure parameter log[(1 —po)/po] 
this advantage seems to disappear. 


(D. F. Morrison) 


6.1 (6.0) 


Residual terms—a source of ordered information in regression analysis—Jn Dutch 
Statist. Neerlandica (1960) 14, 319-341 (24 references, 6 figures) 


Although in regression analysis residual terms seem to 
represent chaos, they can be used to complement the 
original hypothesis in a theoretically acceptable way. 
This is demonstrated with an example taken from the 
construction of a wage-rate equation for the Dutch 
economy 1923-38 and 1949-57. Extrapolation of the 
equation for the years 1958 and 1959 and its application 
to the United States economy 1931-58 provide a primitive 
test. 


(A. P. Barten) 
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BERKSON, J. (Mayo Clinic, Minnesota) 


Nomograms for fitting the logistic function by maximum likelihood—Jn English 


6.7 (11.1) 


Biometrika (1960) 47, 121-141 (8 references, 4 tables, 8 nomograms) 


In this paper the author uses the logistic function as a 
model for bio-assay and other experiments where a 
quantal response is observed at different values of a 
stimulus. In such experiments, as well as estimating 
the two parameters, « and f of the function, it is 
frequently desired to estimate the L.D.59,, given by 
their ratio. 

If r animals out of n die when they are given dose x 
then the maximum-likelihood estimates of «, B and y 
are functions of the two quantities; the sum of r over 
all values of x and the sum of the product rx over all 
the values of x. It was pointed out that, for any specified 
dosage arrangement, maximum-likelihood estimates 
could be tabled against an argument of these two 
quantities. The author points out that for the required 
accuracy the tabular presentation would be formidably 
elaborate and the computation would be quite tedious: 
he has therefore preferred nomographic presentation. 
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BLOEMENA, A. R. 


(Mathematical Centre, Amsterdam) 
Random associations of points on a graph—IJn Dutch 


The nomograms presented provide maximum-likeli- 
hood estimates of y and $f for equally spaced dosages 
and equal n at each dose, for 3, 4, 5 and 6 doses. A 
separate nomogram is given for y and for f. 

The variances of the estimates, obtained from large 
sample formulae, are determined by and the estimates 
of y and f. The author gives in eight pages of tables 
the standard error of each estimate, multiplied by Ait 
to two significant figures. 

Two examples of the use of the nomograms and 
tables are given. 


(N. W. Please) 


6.9 (2.9) 


Statist. Neerlandica (1960) 14, 267-274 (14 references, 1 figure) 


The author supposes a graph consisting of n points, 
numbered from 1, ..., n, and m,; joins between points 7 
and j. From this graph r points are chosen at random. 
This paper deals with the properties of a random 
variable based upon the number of joins between chosen 
points. Applications to a test for randomness, and to 
the order-disorder problem are considered. 

Analogous results for other statistics and also for the 
case in which the points are sampled independently, 
are given. These results are described in Bloemena’s 
report: ‘On probability distributions arising from 
points on a graph” [Statist. Dept., Math. Centre, 
Report S 266A (1960); abstracted in this journal 
No. 1/540, 1.3]. 


(A. R. Bloemena) 
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CSAKI, P. & FISCHER, J. (Math. Inst., Hungarian Academy of Sciences, Budapest) 


On bivariate stochastic connection—Jn English 


Second Hungarian Mathematical Congress (1961) 


In this paper the problems of bivariate connection are 
discussed with the aid of Hilbert space theory. Let 
L* = L°(Q, S, P) be the Hilbert space of random 
variables of finite variance defined on the probability 
space (Q, S, P). Let € and 7 be arbitrary random 
variables with joint distribution function H(x, y) and 
marginal distributions H,(x) and H,(y), respectively, 
which generate the measure P on the plane (x, y) and 
the measures P, and P, on the real line. Let Li, 
(i = 1, 2) denote the space of functions f(x) for which 


f?(x)dH(x) is finite; and L2 = L(Q, S., P) where 
S_ denotes the smallest o-algebra with respect to which 
€ is measurable; L*, is defined analogously. Li, and 
Lz are isomorphic. The sub-space of L? consisting of 
its elements with zero expected values will be denoted 
by Lz 9: Li, o is defined analogously. 

Let A, and A, be the operators transforming the 
elements of L* into their orthogonal projection on Lj 
and Li respectively. Thus the conditional expected 
value of C€L? on y is M(C/n) = A,¢. The maximal 
correlation and the mean-square contingency of € and 
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FISCHER, J. & CSAKI, P. 


(Math. Inst., Hungarian Academy of Sciences, Budapest) 


6.2 (6.9) 


y are S(C, 4) = || A, || = || Ag || and C(g, 1) = II A, Ill 
= ||| Ag ||], whenever the domains of A, and Ag are 
restricted to Lz» and Li, respectively. 

The operator A defined on Lj,, by the isomorphism 
of Lz and Lz, which corresponds to A,/A, is defined 
analogously on Dee 

The equivalency of the following statements is proved: 

(i) P(A | y)<P,(A) and P,(B | x) <P,(B) 

(GUD) Ia aS ee 

(iii) A, and A, are integral operators, where P,(A | y) 
and P,(B| x) are the conditional distributions, A and 
B are Borel sets on the real line: this is Theorem 1. 

It is proved that the above definition of the mean- 
square contingency is equivalent to that of Rényi [New 
version of the probabilistic generalisation of the large 
sieve, Acta Math. Acad. Sci. Hung. (1959) 10, 217-226; 
abstracted in this journal No. 1/252, 6.0; and On the 
measure of dependence, Acta Math. Acad. Sci. Hung. 
(1959) 10, 441-451; abstracted in this journal No. 1/658, 
6.0. 

Theorem two provides a method for replacing any 
given distribution on the plane with maximal correlation 
S by a symmetric one with maximal correlation S?. 

(P. Csaki) 


6.2 (6.0) 


Contributions to the problem of maximal correlation—IJn English 


Second Hungarian Mathematical Congress (1961) 


Let € and 7 be arbitrary random variables defined on 
the probability space (Q, S, P) and S; denote the smallest 
o-algebra with respect to which € is measurable, further 
hg L-(Q, Sz, P). The subspace of Lz consisting of 
its elements with zero expected values will be denoted 
by 12. The o-algebra S,, the spaces Li and Ly, 
are analogously defined. The symbol A, denotes the 
operator of the orthogonal projection of the elements 
of LeN6 on yp the conditional expected value, and 
Az is defined analogously. ; 

The correlation ratio of a standard random variable 
feLz, on n is 0,(f) = || A,f ||, the maximal correla- 
tion of € and is S(E; n) = || A, Il = I Ae ll. 

If A, f = 4g and Azg = Af holds then we call 4 an 
eigenvalue and f, g a pair of eigenfunctions of the pair 
of operators A,, Az. Necessary and sufficient conditions 
are given for feLz,o and géL,, 9 to form a pair of 
eigenfunctions; theorem one. 

In the following the case S(é, 7) = 1 is classified and 
a necessary and sufficient condition for 7 = f(¢) is given. 
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In the case of linearly correlated random variables 
and 4 a necessary and sufficient condition for homo- 
scedastical correlation is described in theorem two. 
Further, theorem three provides a method for calculating _ 
the eigenvalues and eigenfunctions of the pair of | 
operators A,, Az whenever linearly independent systems _ 
b,€Le 03 W,€L;, 9 satisfying certain conditions are 
known. The eigenfunctions are polynomials if and 
only if for each n the nth conditional moment is at most 
an n-degree polynomial of the conditioning variable. 

Finally, examples are presented for calculating the 
maximal correlation. 


(J. Fischer) 


FRECHET, M. (University of Paris) 
Tables which have given margins—In French 
Trab. Estadist. (1960) 9, 3-18 


The author considers a rectangular table of q rows 
and r columns; the element of the J row, i column will 
be designed by n;,=0. 

The margins of this table are the sums: 


q 
2 where XN; = oe 


r 
N= 2) gs Ny = 
veo i= 
The following problem is considered: what may be said 
about a table when the margins are known? More 
precisely, when given the non-negative numbers N ‘> 
N>, ..-, Ng; Ni, .... Ny we have to find 
(i) If there exists a table such that the given numbers 
are their margins. 

(ii) If there are one or more. 

(iii) How to calculate the elements of a table which is 

a solution of the problem; supposing that this 
solution exists. 

The author proves that there exists at least one table 
which has the given numbers as margins. He finds 
an expression which allows him to obtain all the possible 
solutions; special care is given to the extremal solutions. 
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FRECHET, M. (University of Paris) 
Tables which have given margins—IJn French 
Trab. Estadist. (1960) 9, 3-18 


(ii) another application of these properties is to give 
a definition to the distance between two statistical 
distributions F(x) and G(x). 


In the second part of the paper the author considers 
analogous problems, by regarding X and Y as simultane- 
ously determined random variables and examines the 
properties of the correlation between X and Y when the 
probability laws of these variables are known. It is 
proved that there are at least two distinct solutions, 
H°(x, y) and H'(x, y), which are respectively the greater 
and the smaller of the solutions. The author also 
studies the case of a unique solution and he proves that 
the necessary and sufficient condition for the existence 
of such solution is that one at least of the random 
numbers x, y is a constant. In this case the unique 
solution H(x, y) is identical to H°(x, y) and to H'(x, y). 


(J. Béjar) 
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6.9 (6.2) 


As a generalisation, the distribution functions are 
considered: if F(x) is the frequency in which Y¥<x and 
G(y) is the frequency in which Y<y we have to determine 
a distribution function H(x, y) such that F(x) and G(y) 
are the marginal distributions. It is proved that the 
necessary and sufficient condition for this is that H(x, y) 
be a function included between H°(x, y) and H'(x, y), 
where 


0 
FEC = MEX 5) EGG) —1 


eaehX) 
H' (x, y) = Min 
(x, y) = Mi G(y) 
The following two applications to linear programmes 
are considered: 


(i) when the frequency laws of XY and Y are given, 
the minimum and the maximum of the former 
coefficient of correlation of X and Y is the 
coefficient given by H® and H', that is, is given 
by the extremal solutions to the right and to the 
left 


continued 


6.9 (6.2) 
continued 


STOR SOS FOS TCE at TT SO de i 


ee 


-and any point on the estimated curve. 


HEY, E. N. & HEY, M. H. (Radcliffe Infirmary, Oxford and British Museum, London) 


6.7 (4.8) 


The statistical estimation of a rectangular hyperbola—In English 
Biometrics (1960) 16, 606-617 (3 references, 3 tables, 2 figures) 


A simple method is described for fitting a function of 
the type (Y—a)(Y—b) = c to experimental data in two 
variables x and y, both subject to errors of measurement. 
The errors are assumed to be normally and independently 
distributed but in general the error variances may be 
different, not only for the two variables but also for 
different observations of the same variable. Maximum 
likelihood estimates of a, b and c are obtained by 
minimising the sum of the squares of standardised 
residuals in both variables. A direct and rigorous 
solution of this minimal problem is possible only in a 
very few special cases, the most important of which is 
when a straight line is to be fitted to the data and the 
ratio of the error variances is constant for all points. 

Brown [J. Pharmacol. (1952) 105, 139-155] fits a 
rectangular hyperbola to experimental data by minimis- 
ing the sum of the squares of the differences between 
the area of the rectangle defined by the asymptotes of 
the estimated curve and an experimental point and the 
rectangle of constant area c defined by the asymptotes 
The solution 
is not entirely satisfactory; the method gives excessive 
weight to points with high values of x or y. 

Zisol 


HUDIMOTO, H. 


(Institute of Statistical Mathematics, Tokyo) 


The present authors estimate the parameters a, b 
and c by minimising the sum of the squares of standard- 
ised “normal” residuals, that is, distances from the 
experimental points at right angles to the fitted curve. 
Two approximations are introduced. First, the quantity 
actually minimised is an approximation to the sum of 
squared normal residuals. Secondly, estimates of the 
parameters are required for the minimisation. These 
can be obtained graphically from the data. Jmproved 
estimates can be obtained by iteration but the authors 
have not found this necessary in practice as the solutions 
are not very sensitive to variations in the first estimates 
of the parameters. 

The computation procedure is illustrated by a 
numerical example. The estimated curve fits the data 
better than that obtained using Brown’s method. The 
present procedure is well adapted for either a desk 
calculator or an automatic electronic machine. Formule 
for the approximate variances of the estimates are also 
given; the paper concludes with an appendix giving the 
expectation of the best estimate of the root-mean- 
square standardised normal residual, used for testing the 
goodness-of-fit of the estimated curve. 

(V. Chew) 


6.5 (3.8) 


On a coefficient of unidimensional ordering for the individuals’ attitudes—IJn English 
Ann. Inst. Statist. Math., Tokyo (1960) 12, 27-35 (4 references, 4 figures) 


Suppose that responses for m questions (II, ..., I,,) 
which characterise a sociological or a psychological 
object A in some sense are available for the purpose of 
grading of individuals in a certain group according to 
the intensity of attitudes for A. 

In this paper the author considers the dichotomous 
case that an individual’s response to every II; can be 
reckoned to take a score of unity or zero according to 
favour or disfavour; the grading is carried out by their 
total scores. The adequacy of such a procedure is 
discussed, and the comparison of the set of individuals’ 
response patterns which actually occurred with the ideal 
set of individuals’ response patterns (denoting that the 
procedure is completely adequate) is given by the ratio of 
variances calculated from the total scores. It isshown that 
the value of this ratio is non-negative, and does not exceed 
unity. For the cases where the responses to different 
items are statistically independent, it is shown that the 
ratio of covariances is approximately distributed accord- 
ing to a normal distribution for a large number of 
respondent individuals. 


(H. Hudimoto) 
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LANCASTER, H. O. (University of Sydney) 6.2 (2.1) 
On statistical independence and zero correlation in several dimensions—Jn English 
J. Aust. Math. Soc. (1960) 1, 492-496 (4 references) 


The author has extended his earlier theorem on the 
independence of two random variables [Aust. J. Statist. 
(1959) 1, 53-56; abstracted in this journal No. 1/435, 
6.2] to an arbitrary number of dimensions. The follow- 
ing results are obtained: 


Theorem 1. Let {x{'?} and {x$”! be complete ortho- 
normal sets on the marginal distributions of two random 
variables, which have a joint probability measure 
P=P(x,, x2), and let x{° = x$ =1. Then a necessary 
and sufficient condition for the independence of x, and 
xX, is that every coefficient of correlation 


Pisin = [ rasa, 
J 


should be zero for every i, >0, i,>0. 


Theorem 2. The random variables {x,} are mutually 
independent if and only if the generalised coefficients of 
correlation p*!, corresponding to complete sets of 
orthonormal functions, are all zero. 


(J. Gani) 
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NADDEO, A. (Faculty of Statistics, University of Rome) 6.9 (11.0) 
On the measure of statistical dependency between two disjoint qualitative variables—Jn Italian 
Riv. Ital. Econ. Demogr. Statist. (1960) 14, 117-130 (3 references, 6 tables) 


The author, making use of the law of conservation of 
formal properties, establishes first an index of dis- 
similarity between two disjoint qualitative variables. 
Using this result he extends them to the joint distribution 
of the two variables the indices of connection by Gini & 
Castellano; the index of similarity, also by Gini, and 
the Pearson’s correlation ratio, and shows in some cases 
a few simple relations connecting these indices to the 
chi-squared distribution. 


(V. Levis) 
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PERKAL, J. & SZCZOTKA, F. (Math. Inst., Polish Academy of Sciences, Wroclaw) 


6.4 (6.3) 


A new method for analysis of multiple attributes—In German 
Biom. Zeit. (1960) 2, 108-116 (3 references, 8 tables) 


The analysis of multiple characters in classification and 
allied multivariate problems of a mainly descriptive 
nature by factor-analysis methods is cumbersome and 
frequently prone to interpretational fallacies. This 
paper is concerned with a simpler and more straight- 
forward method of extracting a reduced set of “ char- 
acteristic indices ’’ from a set of p correlated attributes; 
sampled on each of n individuals and measured in a 
standardised scale. The np sample variates x, = 1, 
5 DP} J=1, ..., m constitute a pxn matrix x of p 
row vectors of elements x; with zero mean and unit 
mean-square and n column vectors x; with means m,. 
In case (i) of a positive correlation matrix Ry = (1/n) 
(xx’) a single “‘ first order index”’ qj of the jth individual 
is uniquely determined and equal to the individual mean 
Mjs qj =m,. The standardised first residuals d;; = 
(xi;-@))(1/s,), where ns? = [x;—q"[x;—q'}', form a 
p xn matrix D, of p unit row vectors d;. The attributes 
are renumbered and arranged in s groups of p, row 
vectors d‘ in such a way that the s sub-matrices 
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TALLIS, G. M. 


(Commonw. Sci. Industr. Res. Org., Glebe, New South Wales) 


D,, (k =1, ..., s) have positive correlation matrices 
Ry, = (/n)(Dy,Dii), where Xp, Sp. 

Certain rules are given to avoid ambiguities in group- 
ing of the d;. The column means d™ of the matrices 
D,, constitute s sets of second order indices m}*. 
Higher order indices can be obtained in the same way 
as long as the residual correlation matrices contain 
positive sub-matrices. If the correlation matrix Ro is 
non-positive the matrix x is rearranged by rows and 
partitioned in such a way that the sub-matrices have 
positive correlation matrices. 

In a first example (p = 7 human body measurements, 
n = 195 male individuals) the data are fitted by extract- 
ing one first and two second order indices; in a second 
example comparisons with factor methods show that 
this method has a low relative error. 


(R. Wette) 


6.4 (4.1) 


The sampling errors of estimated genetic regression coefficients and the errors of 


predicted genetic gains—IJn English 


Aust. J. Statist. (1960) 2, 66-77 (12 references, 3 tables) 


The purpose of this paper is to develop a theory, general 
yet easy to put into practice, for the evaluation of 
sampling variances of estimated genetic regression co- 
efficients, and genetic gains. 

Two main methods are known for estimating genetic 
regression coefficients. Method (a) uses relationships 
between parents and offspring, and is equivalent to an 
ordinary multiple regression analysis. Method ()) 
depends on relationships between full-sibs or half-sibs. 

Sampling errors are obtained for genetic regression 
coefficients estimated by methods (a) and (6). Sampling 
errors of predicted genetic gains are also found for 
method (a); exact results cannot be derived for method 
(b), but an approximation is given. 

Some special cases of particular interest are con- 
sidered, and an extension of the method to the estimation 
of the mean of the progeny from selected parents given. 
The paper concludes with an illustrative example based 
on some wool data. 


(J. Gani) 
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CALINSKY, T. (Inst. Plant Breeding, Polish Academy of Sciences, Poznan) 


752) 


On a certain statistical method for investigating interaction in several experiments 


with plant varieties—In English 


Second Hungarian Mathematical Congress (1961) 


Put 
ey — m+, +0;+w,,+8,(¢ = 1, 2,....0; j = 1,2,..., k) 


where m, v; are constants, p;, w;; and ¢,, are mutually 
independent normal random variables, each with zero 
mean and with variances 0%, o;,,, 0;/v respectively. 


The paper deals with the statistical investigation of 


the interaction variances o2,,, 04, ..+5 02,- 
Put 
(a 
(n—1)(k—-1) 
n 1 k tat 1 n k 2 
Xy— = > Mo — Yo xyt+ — x 


and let s$ denote the pooled experimental error. In 
this paper the author develops firstly the effective test 
for significance of particular components of interaction, 
based on the “‘ mean-square ratio ” s3/s$, and secondly, 
the approximate test for heterogeneity of interaction, 
2s on the ratio of one sj to another; for example 
Sasa. 
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DANFORD, M. B., HUGHES, H. M. & McNEE, R. C. (School of Avn. Med., 


In both the tests the F distribution is used. While 
the development of the former test is comparatively 
simple, the latter requires a more complicated mathe- 
matical treatment. Among others, the non-central y? 
distribution and the negative binomial distribution are 
applied. 

The practical application of the method is demon- 
strated on the analysis of the results of a series of experi- 
ments with varieties of garden pea. 


(K. Sarkadi) 


7.5 (5.8) 
Brooks Air Force Base, Texas) 


On the analysis of repeated-measurements experiments—In English 
Biometrics (1960) 16, 547-565 (12 references, 6 tables, 3 figures) 


Suppose there are J levels of some treatment, with n; 
subjects exposed to the ith level, and each subject 
measured on some characteristic periodically for K 
times after treatment. The problems considered are 
those of obtaining tests of hypotheses about treatment, 
time, and treatment by time interaction effects. Diffi- 
culties are encountered with performing a univariate 
analysis of variance because the assumption of equal 
variances and covariances is not fulfilled. 

Univariate procedures are discussed, with the applica- 
tion of the techniques to the data of a particular example. 
Multivariate techniques are considered, together with 
their underlying assumptions. The latter are applied 
to the data of the example, and a comparison of the 
univariate and multivariate tests is made. A serial 
correlation model is also considered. It is noted that 
asymptotically the univariate and multivariate tests are 
identical. The conclusions are that, for the example 
given, essentially the same inferences are made from the 
univariate and multivariate analyses. 


(W. A. Glenn) 
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DEDOLPH, R. R. (University of Maryland) 


TO(O.19 


A suggested method for handling data obtained with an exponential (variable dosage) 


sprayer—In English 


Proc. Amer. Soc. Hort. Sci. (1960) 75, 789-798 (6 references, 5 tables, 4 figures) 


With the development of the exponential or variable 
dosage sprayer, described in detail by Dedolph, Basham 
& Stark [Proc. Amer. Soc. Hort. Sci. (1960) 75, 785-788}, 
there has arisen a need for a method of handling data 
taken from experiments utilising this machine. The 
sprayer, designed to facilitate concentration screening 
of growth regulators, herbicides, fungicides, etc., operates 
on the principle of continual dilution, aad applies the 
affecting chemical at a rate which decreases exponentially 
with row distance covered. Randomisation of con- 
centration is therefore impossibie, although the direction 
of spray runs, that is, replications, can be randomised. 
In order to apply the analysis of variance to data of 
this type, one must be willing to assume that the cor- 
relation between adjacent plots, that is to say, between 
adjacent row increments, is negligible. Using yield per 
row increment, or per segment of a row increment, as 
the dependent variable and distance in increments as 
the independent variable, the standard curvilinear 
regression procedures for fitting the response curve and 
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HARTER, H. L. (Aero. Res. Labs., Wright-Patterson Air Force Base, Ohio) 


for estimating confidence limits are applicable. Since 
row distance covered and concentration of chemical are 
functionally related, the response curve for yield can be 
expressed in terms of the concentration. The significance 
of the confidence limits in determining the initial con- 
centration of chemical to use is discussed. 


(T. R. Konsler) 


7.7 (11.1) 


Critical values for Duncan’s new multiple range test—Jn English 
Biometrics (1960) 16, 671-685 (8 references, 1 table) 


The purpose of this paper is to give revised, extended 
and additional tables of critical values for D. B. Duncan’s 
well-known multiple-range test [see, for example, 
Biometrics (1955) 11, 1-42 and (1957) 13, 164-176], 
The table entries are upper 100 per cent. points of the 
studentised range for a normal parent population, with 
parameters p, the number of variates whose range is 
being taken and v, the degrees of freedom on which the 
estimate of error is based. In turn, the “ 
level’ P equals (1—«)?~*, where a is the significance 
level. Entries are to four significant figures for the 
following values of the parameters: p = 2 (1) 20 (2) 40 
(10) 100, v = 1 (1) 20, 24, 30, 40, 60, 120, co and « = 
0-10, 0:05, 0:01, 0-005, 0-001. 


(H. A. David) 
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LE ROY, H. L. & LANDIS, J. (Swiss Federal Institute of Technology, Zurich) 7.4 (—.-) | 
The estimation of missing values in a-b-c- and a-b-experiments—In German | 
Biom. Zeit. (1960) 2, 98-107 (3 references, 3 tables) 


Unbiased estimators Xi ik of nm missing values xX; in (ii) let « = 1—a, B = 1—b and y = 1~—c and deter- 
a-b-c-experiments are obtained by minimising the abc- mine the elements of K according to 
interaction mean-square. The estimation procedure 5 Same ata 

4 ; ; je = —(a t-B JJ! snp tck") 
reduces to the solution of the system of equations a : 
& = K's, where 2 is the vector of estimates, s a vector where 0,,, is the kronecker-delta, that is according 
with arbitrarily ordered elements to concordance (i = i’ etc.) or discordance of 
the marginal row and column indices ijk and 


S.. = @07;,.+ acT.,, + bcT ,—aT..— .— aps : ; 
Jk PCC lay pueda Glee ln Che til.., i'j’k’ respectively. The solution of the system of 


where the 7;;, etc. are the observed sub-totals (i = 1, 
Bee yi ie. bo Ok ula eee. ays and AK - is the 
inverse Of a nxn matrix K, the elements of which are 
easily determined in the following way: 


(i) put the n missing values in the same arbitrarily 
chosen (preferably lexicographical) order as the 
Sj, according to index ijk, let this order deter- 
mine TOWS (Fil, 2 mead columns (s:= 1, :.., 
n) of K, and write down a marginal row and 
column of this indices. 
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equations then follows along the usual lines. 


In case of an a-b-experiment simply put 


and 


_ A bit”, R655" 
kis = a” B a 


and proceed as before. 
given. 


MORLEY, F. H. W. (Commonw. Sci. Jndustr. Res. Organisation, Canberra) 
An effect of genotype x environment interactions on estimation of genotypic covariances—In English 
Aust. J. Statist. (1960) 2, 97-101 (4 references, 3 tables) 


Kempthorne [Introduction to Genetic Statistics (1957) 
New York: Wiley] has discussed the estimation of 
components of genotypic variance in a model in which 
s sires were each mated to a random sample of d dams, 
the progeny of each cross being tested in r replicates 
of a randomised block experiment. 

Certain assumptions are made in this model which 
the author has found inapplicable to his own recent 
studies in alfalfa. A modified model is therefore intro- 
duced, its main difference being that the number of 
replicates may be random if necessary. In this way, 


unbiased estimates can be more readily obtained even 


if interactions are present, providing these are included 
in the new model. Examples of the analysis when 
replicates are random and fixed are given. 


(J. Gani) 


Two numerical examples are 


(R. Wette) 


7.2 (9.1) 
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RASCH, D. (Inst. Tierzuchtf. Landwirt.-wissensch., Berlin) 


7.5 (7.0) 


On the analysis of variance for two-way classification with unequal numbers 


in subclasses—Jn German 


Biom. Zeit. (1960) 2, 194-203 (8 references, 4 tables) 


After an explanation of the orthogonal case of the 
analysis of variance technique for two-way classification 
and an illustration of the fact that the additive property of 


the sums-of-squares is still valid if the numbers in sub- 


Classes are proportionate, the author presents three 
approximate methods for the non-orthogonal case of 
analysis of variance. 


(i) Snedecor’s ‘“‘ method of unweighted means ”’ 
[Statistical Methods (1946) Ames: Iowa State Press, 
pp. 285-301] where the marginal means are 
estimated by x;.. = 1/q ))x;;.. or x!;. = 1/p ¥.x;;., 


J 
respectively, and the total mean-estimate is 
Kee NG Dari 
ty J 
and the error-sum of squares is given by the 


intra-subclass sum-of-squares divided by pq/ > ae 
ig Nj 
squares of 


(ii) Kendall’s “‘method of weighted 


London: Griffin, pp. 220-228] where the variance 
between the X;.. (and the x.,.) is calculated using 


weights g}.=q7/> te OF -g.y=—p 7] >: es and 
7 Ni; T Ni 

the error sum-of-squares is given by the bare 

intra-subclass sum of squares; 

(iii) Cramér’s ‘‘ method of weighted squares of 
means” [Mathematical Methods of Statistics 
(1946) Princeton, pp. 536-547] where the marginal 
means are estimated by 


oe I/n.. Y/N. ;X;j-5 OTe eye = 1/n.. 9) 0}-Xjj-5 
a l 


respectively, their squares are weighted by 
2 2 
ni, Nj. 
gi. = n2/> —, or g"; = n2/¥ —, 
J Nj; i Nj; 
and the error sum of squares is given by the bare 
intra-subclass sum of squares. 


Annotations: x;;, is the kth variable in subclass (i, j), 


means ” [Advanced Theory of Statistics (1946) 2, k = 1,2,....n;;,i= 1,2,..-,p,j = 1,2,...,q. Pointed 
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in subclasses—In German 


Biom. Zeit. (1960) 2, 194-203 (8 references, 4 tables) 


subscripts are written for mean operations in x’s or 
summations inn’s. Marked letters stand for estimates or 
weights. 

It is stated that Snedecor’s mean estimates are not 
as efficient as Cramér’s, and that Cramér’s method 
is more powerful than Kendall’s, if 


ee gn ae GIL g- 


and if an analogous inequality simultaneously holds 
for j. The power of the different methods is illustrated, 
by a numerical example with x the run of food of pigs 
(until to a living weight of 110 kg.), 7 the sex of pigs, 
and j several years. 


(R. K. Bauer) 


2/534 


’ ier , 


stall Te 


Liat ary 4 
nceiahchigg . 


4 Nee a 


RNR TT | Eye PT UCL e Oar HEL el Ow eS eer eae ee 


SSS SSS 


mined by iterative methods. 


ROBSON, D. S. & ATKINSON, G. F. (Cornell University, Ithaca, New York) 


7.6 (6.1) 


Individual degrees-of-freedom for testing homogeneity of regression coefficients in a 


one-way analysis of covariance—In English 


Biometrics (1960) 16, 593-605 (5 references, 4 tables, 1 figure) 


The test procedure for the one-way covariance design 
presented in this paper is constructed specifically to 
have power against some restricted classes of hetero- 
geneous but still linear alternatives to the hypothesis of 
homogeneous (parallel) linear regressions. If the 
regression lines do not intersect at a common point, 
which is an imposed regularity condition, then the 
relationship between slope and intercept will be non- 
linear. Whatever the configuration of ¢ regression lines, 
the slopes may be expressed as a polynomial function 
of the intercepts, of degree at most (t—1). If the 
within-treatment slopes are homogeneous, then the 
degree of polynomial will be zero; otherwise, the degree 
will be greater than zero and will depend on the pattern 
of homogeneity. This fact is exploited by constructing 
the orthogonal polynomial regression of estimated 
slopes on the adjusted intercepts, and then testing 
homogeneity of slope by testing the significance of the 
individual coefficients of the orthogonal polynomial. 
Such a test procedure amounts to partitioning the 
ordinary F-test ““among regression coefficients’? into 
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WILKINSON, G. N. (Commonw. Sci. Industr. Res. Organisation, Adelaide) 
Comparison of missing value procedures—In English 


at most (t—1) individual tests designed to have greater 
power against specific patterns of heterogeneity. Mention 
is made that this procedure compares with Tukey’s 
single degree of freedom tests for non-additivity, and 
the mathematical justification follows that outlined by 
Tukey [Biometrics (1949) 5, 232-242]. 

The theory for testing the polynomial regression of 
slope on adjusted treatment mean is developed, and a 
numerical example showing the computational procedure 
is illustrated. Finally, mention is made of the dual 
result for testing homogeneity of intercepts. Here the 
assumption of constant slopes and variable intercepts 
is replaced by the dual assumption of constant intercepts 
and variable slopes. 

It should be noted that under certain conditions, the 
test developed in this paper also has power against 
non-linear alternatives to the hypothesis of parallel 
straight lines. Some caution must therefore be exercised 
in interpreting a significant F-value in this (or any other) 
test of linear homogeneity. 


(W. H. Beyer) 


7.4 (7.6) 


Aust. J. Statist. (1960) 2, 53-65 (12 references, 5 tables) 


When loss or rejection of some of the observations 
introduces fortuitous lack of orthogonality into the 
data of a designed experiment, a familiar technique of 
analysis is to replace the missing values by estimates 
determined by the method of least-squares. When more 
than one value is missing, the estimates are often deter- 
Alternatively, an adapta- 
tion of covariance analysis, with dummy variates to 
represent the missing values, may be used. 

In this paper, the writing down and solving of the 
simultaneous equations for the estimates, called the 
“¢ direct procedure ” is advocated. This method is set 
out in matrix terms, and then compared with the co- 
variance procedure. The two procedures are shown to 
be equivalent, and the correspondences between them 
are set out. In general, the direct procedure requires 
fewer steps, and in particular gives the correct residual 
sum-of-squares directly. 
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The iterative procedure is also discussed and compared 
with the direct procedure. The advantages of deter- 
mining explicitly the equations for the missing values 
are that (i) they enable exact standard errors for com- 
parisons affected by the missing values to be calculated, 
and (ii) they enable the most efficient and concise method 
of solution to be found. 


(E. J. Williams) 
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GALANTINO, F. (School of Statistics, Palermo) 


8.1 (8.7) 


Statistical method in the survey of social phenomena. Recent contributions 


to some sampling techniques—Jn Italian 


Ann. Fac. Econ. Com. Palermo (1959) 13, 1-98 (91 references, 2 tables) 


In this paper the author analyses some theories of 
sampling and, principally, the sampling theory for 
estimates based on fewer individuals than the number 
selected (Durbin), the theory of stratified and casual 
samples (Deming, Yates), the estimation of an optimum 
sub-sampling number (Brooks), and so on. 

He proposes some sub-sampling methods; with 
reference to the Brooks theory, he considers in the 
general scheme 

Viv = Ytujtei; 
where y;, is the value of the jth element of the ith unit, 
Y the population mean, u; the deviation of the mean 
of the ith unit from the population mean, e;,; the devia- 
tion of the jth element of the ith unit from the mean 
of the ith unit, the hypothesis 


CiCr— ciig= Kk. 


where C is the cost of executing the two-stage sample, 
C,, the cost of selecting a single unit for sub-sampling, 
o2 the variance of u,, 6, the variance of Y, K being 
a constant. 
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JANOSSY, L. & RUPP, E. 


(Hungarian Academy of Sciences, Budapest) 


The number n of units selected for sub-sampling is 
(1+62/ma2)\(1+mC,/C,), 


m being the number of elements sampled within each 
selected unit or sub-sampling number, 62 the variance 
of e;; and C, the cost of sampling an element within 
a selected unit. 

The optimum sampling number (n, of n) may be 
calculated when the optimum sub-sampling number 
(m, of m) and the ratios between variances and between 
costs are known. 

The author, however, generalises the known scheme 
C=nC,+nmC, by proposing the relation C=nC, 
+nmC, where 


ey |S 
yv=1 
is the mean of sub-sampling numbers. 

Generalising the questionnaire scheme of Gallup, the 
author shows some new practical schedules for social 
survey, with reference to problems of sampling pro- 
cedures for mailed questionnaires (El Badry). 


(F. Galantino) 


8.1 (4.3) 


Determination of parameters in case of exponential decay—Jn Hungarian 
Magy. Tud. Acad. K6zponti Fiz. Kut. Intézet. K6zl. (1960) 8, 75-81 (3 references, 3 figures) 


In this paper the authors divide the measuring time into 
several intervals and estimate the parameters of the 
exponential decay by means of the numbers of particles 
arriving during each interval. The suitable choice of 
the dividing points is considered. This is shown to be 
essentially equivalent to the special case of the optimum 


_ one-dimensional stratification problem (proportionate 


allocation) of Dalenius [Sampling in Sweden (1957) 
Uppsala: Almqvist & Wiksell]. 

A graph is given, representing the suitable stratifica- 
tion points as function of the first stratification point. 
If the preliminary value of the life-time parameter is 
known, the stratification points can be determined, by 
use of this graph, for given measuring time and number 
of intervals. In a special case, namely for the division 
of the interval (0, 00) into two intervals, they obtain 
the same result as Dalenius. 

Investigations concerning the efficiency lead to the 
conclusion that a measure-time which is ten times the 
mean life-time and its division into four intervals are, 
in practice, sufficient. 


(K. Sarkadi) 
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SARKADI, K. (Math. Inst., Hungarian Academy of Sciences, Budapest) 


8.8 (2.4) 


A rule of dualism in mathematical statistics—Jn English 
Acta Math. Acad. Sci. Hung. (1960) 11, 83-92 (19 references) 


Let us denote by ¢ the fraction defective, by «, the 
number of defective items found in a sample of size n 
taken from an infinite population. The “rule of 
dualism” states a connection between the confidence 
limits and the inverse probability limits of the fraction 
defective. In the case of uniform a priori distribution 
this can be expressed in the form 


() Pr(€<p| x, =k) = Pr(kye2k4+1| C= p), 


k =0, 1, 2, ...,.n; O<p<1. 


If the lot is finite and contains N elements, then the 


corresponding relation is as follows: 
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The above relations are treated by the author who 
also deals with several generalisations of the binomial 
case (i) and the hyper-geometric case (ii). 

One generalisation concerns the poly-hypergeometric 
extension of (ii) for the case when more than two 
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THIONET, P. (Inst. Nat. Statist., Etudes Econ., Paris) 


alternatives are considered. As a limiting form the 
polynomial case is obtained. A further generalisation 
is the Poisson case when, in a continuous production, 
the number x of defects in a lot follows the Poisson law. 
In this case 
PrA=ei| K) Pri lel re) 

holds, where 2 is the expected value of the Poisson 
distribution having a priori uniform distribution over 
the half-line (0, oo), that is on the interval (0, L); Loo. 
The final generalisation concerns the mixture of Beta 
distributions as the a priori distribution. 

The paper concludes with some remarks.on the same 
acceptance method of Oderfeld which concern the 
determination of the parameters of the a priori Beta 
distribution on basis of the past sample results. Some 
of the results have been previously published in 
Hungarian [Magy. Tud. Akad. Alkalm. Mat. Int. Kozl. 
(1953) 2, 275-293]. The proofs are based on prob- 
abilistic arguments instead of on direct computations. 


(I. Vineze) 


8.3 (4.6) 


The adjustment of sampling results to results of counts—Jn French 
Rev. Int. Statist. Inst. (1959) 27, 8-25 (11 references) 


This article reverts to a problem already studied by 
numerous authors from Deming & Stephan [Ann. Math. 
Statist. (1940) 18, 427-444] through to El Badry & 
Stephan [J. Amer. Statist. Ass. (1955) 50, 738-762]: the 
problem considered is how to improve a sample estimate 
by the supplementary information, from two independent 
samples, contained in the marginal distributions of the 
population when sampled with respect to two different 
stratification variables. 

If the population was known in its entirety, this 
would be:a classical case of “‘ re-weighted ”’ estimation, 
also called stratification a posteriori. In attempting to 
reduce the problem to this case by first reconstructing 
the distribution from the marginal ones the author has 
retained two particular cases from all possible distribu- 
tions as studied by Fréchet. 

The first solution of a certain quadratic programming 
problem, is the one which minimises the variance of 
the estimator concerned, irrespective of possible bias; 
the second is that of Morgenstern-Gumbel, adjusted by 
maximum likelihood methods. Experimentally a close 
analogy is noted between the results of the second method 
and the two routine procedures already known: their 
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common fault would be that they do not produce a 
distribution which resembles the sample. The same 
could be said of any other method: and one is led to 
suppose that this absence of resemblance is due to 
deformations of the sample. This paper departs from 
the classical theory of “‘ good” samples and gives pro- 
cedures for treating “bad” samples: it is essential 
that hypotheses about the nature of the deformations 
should not be in contradiction with hypotheses on 
which any given method is based. 

The first method seems to be appropriate if the 
deformations of the rate of sampling by non-response 
vary much from one cell of the contingency table to 
another. On the other hand, the second, which supposes 
the sample-units to be independent, would be appropriate 
if the sample is too small without being systematically 
deformed with respect to the stratification variables. In 
addition a third method is studied which supposes that 
all characteristics are in a linear regression with the two 
stratification variables. The separate samples of these 
variables allow of an estimation by regression: the 
regression plane being estimated from the sample. 


(P. Thionet) 


care 


10M ahr 
Rai tee hue 
hee Ata 


4 he 


hae 


frat Vd 
part ten 


BRADLEY, R. A. & PENDERGRASS, R. N. (Virginia Poly. Inst., and Radford Coll., Va.) 


Ranking in triple comparisons—Jn English 


9.5 (9.0) 


Bull. Int. Statist. Inst. (1960) 37, III 229-241 (17 references, 5 tables) 


A mathematical model for triple comparisons has been 
developed as an extension of the earlier developed 
model for paired comparisons [see Bradley & Terry, 
“The rank analysis of incomplete block designs. I. 
The method of paired comparisons”’, Biometrika (1952) 
39, 324-345]. Parameters 7,, ..., 2, are associated with 
probabilities of rankings of ¢ treatments in n repetitions 


of the possible [ 


Based on the model postulated, tests of treatment 
equality are outlined and a test of agreement of treat- 
ment rankings over sets of repetitions of triple com- 
parisons is given with tests of the appropriateness of 
the model. 

Estimation is based on maximum likelihood pro- 
cedures and large-sample distributions of estimators, and 
test statistics are based on likelihood ratio theory, 
The suitability of these large-sample distributions as 
approximations for moderate size samples is considered. 

A brief example illustrates the techniques developed. 


triple comparisons. 


(R. A. Bradley) 
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DeBAUN, R. M. & CHEW, V. . (Amer. Cyanamid Co., Stamford, Connecticut) 


9.4 (7.3) 


Optimum allocation in regression experiments with two components 


of error—In English 


Biometrics (1960) 16, 451-463 (2 references, 6 tables, 1 figure) 


In the literature, experimental designs for regression 
analysis are usually optimised with respect to the total 
number of observations to be taken. In this paper cost 
functions are introduced and optimum designs are 
derived for both extrapolation and interpolation, in- 
cluding the case where “‘ split-plotting ” occurs. 

The authors discuss the situation in which replication 
in full plots are very expensive but replications within 
plots are relatively cheap, as they might be when 
aliquots of experimental material are the within-plot 
elements and whole lots of experimental material 
represent the whole plots. The number of sub-plots 
within the whole plots for greatest efficiency (in the 
sense of most information per unit cost) depends on 
the magnitude of between- and within-plot errors as 
well as on the incremental cost of additional within-plot 
replications. This general subject has been studied 
extensively in terms of the usual experimental designs 
but this paper is primarily concerned with the problem 
of regression experiments, particularly those in which 
extrapolation is the principal purpose. Attention is 
called, by the authors, to the common use of the method 
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of “extrapolation to infinite dilution” used to find 
the molecular weight of polymeric substances. 

The primary section presents the basic development 
which includes expressions and tables for optimum 
group size and efficiency for estimation of means in 
terms of the ratio of the between-plot variance to the 
total variance and the proportional cost of taking an 
additional sub-plot observation. It is also noted that 
the efficiency function is relatively insensitive to variations 
in plot size in the neighbourhood of the optimum. The 
authors then go on to develop expressions to specify 
the efficiency and the most efficient designs for both 
interpolation and extrapolation for the following cases: 

(i) Regression on one independent variable—In- 
dependent observations, 

(ii) Regression on one independent variable—Non- 

independent observations, 

(iii) Bilinear regression—Non-independent errors. 

A numerical example of last case is presented to 
illustrate various features of the technique. 


(R. J. Taylor) 
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RAO, C. R. (Indian Statistical Institute, Calcutta) 


9.0 (7.6) 


Experimental designs with restricted randomisation—In English 
Bull. Int. Statist. Inst. (1960) 37, III 397-404 (2 references). Discussion abstracted No. 2/544. 


In this paper some designs for experiments involving 
several groups of treatments, providing precisions of 
different order for within group and between group 
comparisons are discussed. 
The principle adopted is that of restricted randomisa- 
tion. The treatments whose differences are to be 
estimated with higher precision, occupy neighbouring 
(contiguous) plots. Thus instead of randomising all 
treatments over all the plots in a block, the block is 
split up into sub-blocks of contiguous plots, the groups 
of treatments are assigned to sub-blocks at random and 
within each sub-block the treatments within a group are 
randomised over the plots. Here the within-group 
treatment comparisons are estimated more precisely 
than the between-group treatment comparisons. Some 
other methods of restricting randomisation to attain 
similar results are also suggested. 

Under this type of randomisation, the analysis is 
derived for the following cases: 


(i) When two groups of v treatments each are 
arranged in randomised blocks of size 2v. 
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RAO, C. R. (Report/of discussion on a paper by) 


(ii) When two groups of v, and v, treatments 
respectively, are arranged in blocks of size 2k 
so that each block contains k treatments from 
each group and any two treatments from the ith 
and the jth group occur in 4;, blocks (i, j) = 1, 2. 


(K. R. Shah) 


9.0 (7.6) 


Experimental designs with restricted randomisation—Jn English 


Bull. Int. Statist. Inst. (1960) 37, I 56-57 


Those taking part in this discussion were Finney, 
Madow and Tippett: the author briefly replied. 

Finney asked for further information on the question 
of the method of comparing two experiments: for 
example, was there a single scale of money or some 
nutritional standard? The author’s reply was to the 


effect that the experiments were designed primarily 


for evaluation of the response surfaces as a function 
of the proportion of mixture. Once these were available 
other things could be studied such as the effect suggested 
by Finney. 

Madow suggested that the response surface type of 
design seemed to be most appropriate and Tippett 
commented that there would be some difference if the 
mixed crops were harvested separately or in their mixed 
form: the author confirmed that the harvesting was 
separate. 


(W. R. Buckland) 
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RUMEKE, C. L. (Dept. Pharmacol. and Neurobiol., Free Univ., Amsterdam) 


9°3:(5.2) 


An efficient design for comparing the effects of two treatments—In English 
Symp. Quantative Methods in Pharmacology (1960) 


The variability of the reactions of groups of animals 
to drugs often interferes with an efficient comparison 
of the effects of two treatments with all-or-none responses. 
An efficient design for an experiment in which the 


probabilities of an effect after two treatments may be 


compared has been derived from Dixon and Mood’s 
“up-and-down ”’ method. Two groups of animals 
which received the two different pre-treatments to be 
compared are treated with logarithmic equidistant doses 
of a convulsant agent. In each group the series of 
experiments is interrupted as soon as the sign of the 
effect has changed three times. The average of the 
two means of the last two doses administered is cal- 
culated. Thereafter n animals of group I and n animals 
of group II are given this average dose. A test for the 
difference between the probabilities of an effect after 
the two pre-treatments is performed in a fourfold 
contingency table at a dosage at which the proportions 
of animals showing an effect may be expected to be 
distributed approximately symmetrically round 0°5. 
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This “‘ semi-sequential ’? method leads to a test for 
the comparison of the effects of two treatments with 
all-or-none responses. However, an estimate of the 
difference between the ED5;,)’s cannot be obtained by 
this method. This opportunity arises by the following 
modification of this design—though with some loss of 
efficiency. 

The experiment is started in the way described. The 
average of the two means of the last two doses admin- 
istered is determined. Thereafter two doses are chosen, 
which have about the same difference from this average. 
In group I as well as in group II 4n animals are to be 
treated with the lower dose and a same number of 
animals is to be treated with the higher dose. A com- 
bined test is performed on the two fourfold contingency 
tables which are to be constructed from the results. 
From the proportions of successes found after the 
different dosages, the ED59’s can be obtained with the 
aid of a probability paper or by some calculation using 
the arc sin root transformation. 


(C. L. Riimke) 
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AKAIKE, H. 


(Institute of Statistical Mathematics, Tokyo) 


10.6 (-.-) 


On a limiting process which asymptotically produces f * spectral 


density—In English 


Amn. Inst. Statist. Math., Tokyo (1960) 12, 7-11 (4 references) 


In some of the reports about the spectral analyses of 
roughness of runways or roadways the spectral densities 
approximately of the form c/f? (f: frequency, c: 
constant) are observed. In this paper a simple model 
is tentatively suggested to explain how this 1/f? pheno- 
menon takes place. 

The model shows that by applying the roughening 
and smoothing operations alternatively and repeatedly 
to the stochastic process, which is taken to represent 
the profile of a runway or roadway, the spectral density 
function of the process tends to some limiting form. 
In this model the roughening operation is represented 
as an addition of a stationary stochastic process of 
fixed type to the former profile and the smoothing 
operation as an application of a convolution trans- 
formation to the profile. By taking into consideration 
the physical meanings of roughening and smoothing, 
some natural assumptions are imposed on the form of 
the spectral density of the process which is added to the 
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BELYAEV, Y. K. (Moscow) 


profile and on the form of the smoothing function to 
get the 1/f* form. This approximation of the 1/f 7 form 
is only valid in the range of f>f,>>0 (fo: arbitrary 
but fixed). 

During the discussion of the model the author indicates 
the necessity of the adoption of some trend elimination 
procedure in the spectral analyses of runway or roadway 
profiles. 


(H. Akaike) 


10.1 (—.-) 


An example of stochastic process with mixing—ZJn Russian 
Teor. Veroyat. Primen. (1961) 6, 101-103 (3 references) 


The example given in this paper is a stochastic process 
é(t) with a continuous parameter and mixing, whose 
range consists of two states. It is shown that the 


Pp 
integral ¢(p) = | €(t)dt does not have any increase in 
(0) 


variance. 
(Y. K. Belyaev) 
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BEN-ISRAEL, A. & NAOR, P. 
A problem of delayed service. I.—Jn English 


(Technion, Israel Inst. Tech., Haifa) 


10.4 (2.5) 


J. R. Statist. Soc. B (1960) 22, 245-269 (4 references, 2 tables) 


A system is considered in which an operator is in charge 
of a number of machines which are liable to break down. 
The operator walks along the machines in a prescribed 
order, attends to them, and whenever he finds a broken- 
down machine he restores it to a working condition. 

It is assumed that attendance time at a machine is 
independent of whether the machine is working or not. 
This means that repair time is negligible compared to 
time spent either walking from machine to machine or 
on inspection and maintenance of all machines. Two 
forms of breakdown distribution are considered: firstly, 
Poisson I where the time elapsing between successive 
breakdowns in the system is a random variable possessing 
an exponential distribution. Secondly, Poisson II where 
the uninterrupted working of a machine is a random 
variable with an exponential distribution. 

Two types of models are considered: 


(i) Simultaneous models where the operator arrives 
at the set of machines and those machines which 
are broken-down are repaired within a time very 
much shorter than the time elapsing between 
successive arrivals of the operator. 
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BEN-ISRAEL, A. & NAOR, P. « (Technion, Israel Inst. Tech., Haifa) 


A problem of delayed service. II—Jn English 


(ii) Successional models where the machines are 
inspected, and possibly serviced, in some pre- 
scribed order. 


Distribution functions are defined for the inter- 
arrival times of the operator, the “ round-time distribu- 
tion”’ and for the time that elapses from the actual 
breakdown of a machine up to the time of repair. The 
relation between these two is found for simultaneous 
models with both Poisson I and II forms of breakdown. 
For successional models only the Poisson II case is 
treated fully, whilst some of the assumptions are slightly 
modified for Poisson I breakdowns. Three particular 
forms of round-time distribution are considered, 
exponential, constant and Erlangian. 

Two further random variables are considered; the 
number of idle machines as observed by a controller 
who observes the system at random instants, and the 
number of broken machines encountered by a repairman 
during one round. The means and variances for these 
variables are given for most of the situations. 

See also abstract No. 2/550. 


(C. Burrows) 


10.4 (2.5) 


J. R. Statist. Soc. B (1960) 22, 270-276 (2 references) 


_ This is a continuation of the previous paper [J. R. 


Statist. Soc. B (1960) 22, 245-269; abstracted in this 
present journal No. 2/549, 10.4] and presents a model 
which is intermediate to those considered earlier. 

It is assumed that the machines are attended to in 

groups, repairs of machines in the same group being dealt 
with in succession. Both Poisson I and IL forms of 
breakdown are considered. 
' Expressions are developed for the mean and variance 
of the number of idle machines as observed by a random 
controller and as observed by the repairman during one 
round. 


¢ 


(C. Burrows) 
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BIRNBAUM, Z. W. (Statist. Res. Lab., Univ. Washington, Seattle) 10.1 (10.4) 
Life length and failure of materials interpreted as stochastic processes—Jn English 
Bull. Int. Statist. Inst. (1960) 37, II] 214-217 (1 reference). Discussion abstracted No. 2/552 


ee ee 


A stochastic model is presented for life-lengths of 
structures which are subject to failure attributed to 
fatigue, brittle failure and similar phenomena. 

The model in this paper assumes that the failure rate 
of the structure at any time t is determined by the 
state of the structure St and the instantaneous damage 
At affecting the structure at that time, and that either 
St and At may be a fixed function of time or a stochastic 
process. It is shown how, under certain conditions, 
the probability distribution of the life-length of a structure 
can be derived from information on St and At. 

The results of these derivations in specific cases are 
reported: see also Birnbaum & Saunders [J. Amer. 
Statist. Ass. (1958) 53, 151-160]. 


(Z. W. Birnbaum) 


2/551 


BIRNBAUM, Z. W. (Report of discussion on a paper by) 
Life length and failure of materials interpreted as stochastic processes—Jn English 


Bull. Int. Statist. Inst. (1960) 37, I 99 


_ Those taking part in this discussion were Gumbel and 


Hemelrijk: the author replied briefly. 

Gumbel stated that most physicists believed in a 
difference in principle between static and dynamic 
loading. The irreversible molecular changes in the 
latter gave rise to the question as to whether the existence 
of the minimum life and endurance limit was foreseen. 
The author replied that it was possible to obtain these 
in accordance with the model proposed in the paper but 
it was difficult to decide whether material counterparts 
existed. 

Hemelrijk asked for clarification of certain symbols 
in the development of the argument and whether some 
variables were random or not. The author referred to 
a paper [J. Amer. Statist. Ass. (1958) 53, 151-160] and 
clarified the detail of the points raised. 


(W. R. Buckland) 
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BLOEMENA, A. R. (Mathematical Centre, Amsterdam) 
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10.4 (10.1) 


On queueing processes with a certain type of bulk service—JIn English 
Bull. Int. Statist. Inst. (1960) 37, III 219-227 (13 references) 


Customers arrive at a counter to be served and line up 
in the order of their arrival. The time between two 
successive arrivals has an exponential distribution with 
unit mean. A batch of (at most) n customers is served 
at the same time. The service times, which are supposed 
not to depend on the size of the batch, are identically 
distributed with mean yu. All inter-arrival times and 
all service times are mutually independent. Service is 
stopped when no customers are present and resumed 
at the moment a new customer arrives, a situation 
which often occurs in practice; for example, in the 
case of lifts, shiplocks, etc. In this respect the situation 
differs from that studied by Bailey [J. R. Statist. Soc. 
B (1954) 16, 80-87] and Downton [J. R. Statist. Soc. B 
(1955) 17, 256-261; and (1956) 18, 265-274], where 
service is supposed to be continued even when no 
customers are present. 

The queue-lengths at the moments that the serving 
of a batch is just completed constitute a Markoff chain, 
which is seen to be irreducible and aperiodic. The 
chain is ergodic if <n. 
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BRODY, S. M. (Kiev) 


The stationary state distributions of the queue-length 
and the waiting-time are studied. The generating 
function of the queue-length distribution and the 
Laplace transform of the waiting-time distribution are 
obtained by interpreting these functions for values of 
the argument between zero and unity as probabilities. 


(F. W. Steutel) 


10.4 (2.5) 


On an integro-differential equation for systems with t-waiting—Jn Ukrainian 
Dopoy. Acad. Nauk. Ukrain, SSR (1959) 571-573 (3 references) 


The author considers the Poisson’s flow of demands 


on the service of one instrument. The time of service 
is a random variable having exponential distribution. 
The demand entering the system is served immediately 
if the instrument is not engaged and has to wait in line 
if the instrument is engaged. The demand waits for 
service no longer than the time t, after which it is lost. 
Integro-differential equations for the time of waiting 
for service are obtain in the paper. 


(S. M. Brody) 
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. of the random variables /,, /,, 


CHIANG, C. L. (University of California, Berkeley) 


A stochastic study of the life table and its applications. I. Probability distributions 


of the biometric functions—In English 
Biometrics (1960) 16, 618-635 (12 references) 


This is the first of a series of papers which will deal 
stochastically with life tables by treating the biometric 
functions in the tables as random variables. In this 
investigation, distributions, expected values and co- 
variances of these functions are derived. A second 
paper of the series will be devoted to the derivation of 
formulas for the sample covariances. A third paper will 
be concerned with applications to studies of patients 
afflicted with certain diseases. 

The author discusses briefly the two forms of life 
tables in common use: (i) the cohort (generation) life 
table and (ii) the current life table. In the derivation 
of the distribution of /,, the number of individuals 
surviving the age interval (0, x), the age x is treated as a 
continuous variable. Consequently, based on the 
derived distribution, a formula is obtained for the 
probability of /; individuals surviving any specified age 
interval in the life table, (0, x;). The joint distribution 
---9\ ly, given that there 
were /, individuals at age 0, is shown to be the product of 
w binomial distributions: and the author shows that the 
mumber of deaths do.=\05,.d; = 6;, ...3.d, = 0, in the 

2/555 


CISTYAKOFF, V. P. (Moscow) 


10.9 (2.4) 


life table have a multinomial distribution with prob- 
abilities of po;q; of an individual dying in the ith interval 
(X;, X;+1)3 Po; being the probability of surviving the age 
interval (0,.x;), q; the probability that a member of 
age x; will die in the ith interval, and /) the size of the 
original cohort. 

The expectation of the proportion of deaths in the 
age interval (x;, x;.,) is given in the fifth section. For 
large values of J), an approximate formula is derived 
for the variance of the observed proportion of survivors, 
p;. The author shows that the covariance between the 
proportion of survivors in two disjoint age intervals is 
zero but points out that 6; and ,; are not independent. 

The distribution of é,, the observed expectation of life 
at age x, is discussed: in addition to the general expres- 
sion for é,, the computational formula for é, is given . 
based on the assumption that the distribution of deaths 
in each interval is uniform. A_ theorem states 
that as /, increases the distribution of the observed 
expectation of life at age x, is asymptotically normal: 
formule for the mean and variance of this asymptotic 
distribution are given. 

(B. J. Trawinsik) 


Transient phenomena in branching stochastic processes—Jn Russian 
Teor. Veroyat. Primen. (196!) 6, 31-46 (7 references) 


Let u,(t) = {uy,(t), +--+) Myn(t)} be a branching process 
with n types of particles and let 
Pr {u,,@Q) = @;; j= 1,..., 0} = OP + petto(t) 

GAM n) 


when 7-0. Here @ = {a,, ..., @,}, 0, =1 for o, = 1, 


wo, =0 (j#k) and 6; =0 in other cases. Having 
defined the generating functions by 
Ce Pad, Gay aN yi. geerae ume Cas Sacre) 


the author denotes the factorial moments by 


= i b&® = Of 
Ox; |x = ee bets Ox;OX ; |x = A 
(heer | 
ee Ox;Ox OX, |. a 
Let 2 be the compact set of an undecomposable 
matrix a = || a, || (k, j= 1, -.., m); 4 = max (Re d)) 


1Sisn 
where the numbers /; satisfy the equality | a—A,;E| = 0 
(E being the unity matrix) and v = {v;}7- 1, u = {uj}i=1 


aij 
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10402) 
satisfy the equalities aw = Au, va = Av, 
SS pe = 1, y U0, = 1. 
k=1 k=1 
Let K(%, 6, B, c) be a class of {f,(x)} with a e Y, 
ead a rs bie B 005 el es e100, 


esa 
The following asymptotic formula for too, 4-0 holds 
true uniformly for all {f,} eK 

1—Pr{u;,(t) a 0, j a i vey M | L;>O}~uk(t, A, 0), 


n 


hy oh, Li; ,(t). 


j= 


The probability distributions 


SEO; OhOF) Yn) 
la ‘| A) n>} 
M(t; | Hx > 9} 


converges to an exponential distribution as f—>oo, 
4-0 uniformly for all { f,} € K. 


<y;j=l,....n 


(V. P. Cistyakoff) 
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CONOLLY, B. W. (Saclant ASW Res. Center, La Spezia, Italy) 
Queueing at a single serving point with group arrival—In English 


OO 


10.4 (2.5) 


J. R. Statist. Soc. B (1960) 22, 285-298 (8 references, 4 tables) 


This paper studies the effect on a single server queueing 
system of arrivals taking place in groups of a fixed size 
instead of singly. The time intervals separating the 
arrivals of successive groups are independent samples 
from a fixed distribution. The service time of each 
individual is assumed to be negative exponentially 
distributed. The groups are dealt with on a first come, 
first served basis, though the order of service of 
individuals within a group is indifferent. A steady state 
is also assumed. This system is formally equivalent 
to the queue given in Kendall’s notation as G//E,/1. 

By considering the instants of arrivals of groups, a 
set of integro-difference equations for the state prob- 
abilities is set up. These are solved using Laplace 
transforms and the class of the inter-arrival distribution 
is restricted to those which give a particular form to 
this solution. Formulae are then derived giving the 
distribution of queue length, the mean idle time of a 
server, the mean waiting time for a customer and also a 
generating function for the joint probability and 
probability density p,(t) that a busy period, that is 
when the server is fully occupied, lasts for time ¢ and 
then # groups are served. 

DiOST 


CSISZAR, I. (Eétvés University, Budapest) 


A numerical example is given with group sizes of 1, 
2 and 4 but with the proviso that the arrival intensity 
for individuals remains constant. The inter-arrival 
intervals are taken as constant. This shows that when 
the group size varies but the mean number of arrivals in 
unit time and the average length of a service period 
remain constant then: 


(i) average queue size increases with increasing 
group size; 

(ii) the average non-zero idle period of the server 
increases with increasing group size, but if the 
average is calculated over all occasions when a 
customer completes service it remains constant; 

(iii) the average waiting time of the first member to 
be served of a newly arrived group increases 
when the average only takes into account non- 
zero waiting times, but decreases when it includes 
all occasions when a group arrives; 

(iv) the average number of groups served during a 
busy period decreases with increasing group size, 
though the number of customers served increases. 


(C. Burrows) 


10.5 (1.5) 


Some remarks on the dimension and entropy of random variables—In English 
Second Hungarian Mathematical Congress (1961) (2 references) 


The dimension of a random variable € is defined as 
the limit d(é) = lim H,(€,)/log, provided this exists; 


where H,(é,) denotes the Shannon entropy of the 
discrete random variable €, = 1/n[né]. The d-dimen- 
sional entropy of € is defined by 

H,() = lim {Ao(¢,)—4 log n} 


always provided that the limit exists and it is finite; 
[x] denotes the integral part of x. 

These notions were introduced by Balatoni & Rényi 
(“Remarks on entropy”, Publ. Math. Inst. Hung. 
Acad. Sci. (1956) 1, 9-40; and also Rényi, “ On the 
dimension and entropy of probability distributions”’, 
Acta Math. Acad. Sci. Hung. (1959) 10, 193-215; 
abstracted in this journal No. 1/313, 10.5]. 

According to theorem 1, if Ho(¢,) = Ho([€]) is finite, 
there exists the limit 

lim {Ho(é,) log n} 
which is either finite or equal to —oo. In the proof 
the authors use only elementary properties of the 
entropy. If the dimension of € is absolutely continuous 


with the density function f(x) then the above limit is 
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— 00 
where the integral surely exists: for the case of finite (, the 
equality of these two limit expressions was proved by 
Rényi. 

Theorem two states that if the distribution of € is 
not absolutely continuous then the first limit is equal 
to —oo, that is to say the one-dimensional entropy of 
€ does not exist. For the proof it is shown that if the 
distribution of € is not absolutely continuous then there 
is an e>0 such that for every 6>0 and for every 
n>n,(o, €) there exist such integers k,, j=1, ..., 5, 
that s<no and 2 Pr(é, = k/n)>«. 

The sum representing H,(€,) is composed of two parts, 
the first of which consists of the members corresponding 
to the k;—j; j=1,.-..,s. Using the Jensen inequality 
we obtain an estimation which proves the statement. 

The second theorem cannot be sharpened: 

H(¢,)—log n 
may tend arbitrarily slowly to —oo for appropriate 
singularly distributed random variables. The results 
obtained can also be generalised for random vectors. 


(I. Csiszar) 
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10.2 (1.6) 


Statistical analysis of a stochastic difference equation—Jn Ukrainian 
Dopov. Acad. Nauk. Ukrain, SSR (1959) 120-124 (1 reference) 


The author considers a random process x, where f is 
an integer satisfying the difference equation 


XT O4X,-4+...+4,xX,_, +a = Hgcet Dagan gt: -Dycr sy 


where a; and b; are unknown parameters and €, are 
independent variables with a mean equal to zero and a 
variance equal to unity. 

The parameters a;, b; are estimated by a sequence of 
observations x,, x2, ..., Xy of the process. The case of 
biased estimates is discussed and estimates of the a;,- 
parameters are given, these estimates being consistent 
and their joint distribution tending toward the normal 
law when Noo. 


(A. I. Dorogovtsev) 
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DREMICK, R. F. 


(Bell Telephone Labs., Murray Hill, N.J.) 
Mathematical aspects of the reliability problem—ZJ/n English 


10.1 (1.2) 


J. Soc. Indust. Appl. Math. (1960) 8, 125-149 (18 references) 


This paper deals with certain theoretical problems 
encountered in the study of equipment reliability. The 
primary purpose of the paper is to report or attempt to 
establish some uniformity of approach to the reliability 
problem. A new concept is introduced which is more 
general than those in common use but which specialises 
to the usual uses under appropriate conditions. 

In the author’s formulation it is visualised that an 
enterprise is to be undertaken which consists of operating 
a piece of equipment or a group of equipments of the 
same kind, for a certain period of time. The entre- 
preneur expects a certain economic return from it, 
but he also expects his return to be adversely affected 
by the incidence of malfunction. He will have established 
in his mind an exact idea of the undesirability to him 
of every kind of poor performance of his equipment. 
He can then define the reliability of his equipment in 
this operation as the mean gain, or loss, which he can 
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expect of it. In this paper this mean gain is adopted 
as the definition of reliability. The paper brings together 
many of the concepts and problems of reliability and 
devotes considerable space to the relationship with 
renewal theory. 


(W. T. Wells) 
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FINCH, P. D. (London School of Economics) 10.4 (2.5) 
On the transient behaviour of a simple queue—In English 
J. R. Statist. Soc. B (1960) 22, 277-284 (5 references) 


A single server queueing process is considered with 
Poisson input and a general distribution of service time. 
The assumption is made initially about the queue 
discipline although it is assumed that the server is never 
idle when there are customers to be served. 

The generating function for the probability distribu- 
tion of queue size after the nth departure from the system 
is found. An explicit solution is given for the case 
where the service-time distribution is a negative 
exponential. The relation between queue size on an 
arrival and a departure is also considered for first come, 
first served queue discipline. Limiting distributions are 
given for both negative exponential and general service- 
time distributions. 


(C. Burrows) 
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GHOSAL, A. (Central Fuel Research Institute, Bihar, India) 10.4 (2.5) 
Problem of emptiness in Holdaway’s finite dam—lIn English 
Bull. Calcutta Statist. Ass. (1960) 9, 111-116 (6 references) 


_ A probability theory of dams has been developed by 
Moran, Prabhu, D. G. Kendall, Gani and others: see 
Moran [Theory of Storage (1959) London: Methuen]. 
In this paper, the probability V(u) of the dam drying 
before it fills up, for a given initial level u, has been 
derived for Holdaway’s finite dam, assuming the input 
to be exponential. 

If x, be the input during the interval (¢, +1), and z, 
the content of the dam at time f¢, the release scheme of 
Holdaway’s dam is as follows: given /<k,d<m 


(i) If z,+x,2k+m, release z,+x,—k 
(ii) If /+mSz,+x,Sk+m, release m 
(iii) If /+d<z,+x,S4+m, release Z,+x,—/ 
(iv) If dSz,+x,S/+4+d, release d 
(v) If z,+x,Sd, release x,+ 2, so that z,,, = 0. 


It is shown that V(u) satisfies an integral equation 
and a method of solving the equation is discussed. The 
results, however, are not in a form which can be 


immediately applied in practice. 


(C. $8. Ramakrishnan) 
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FURSTENBERG, H. & KESTEN, H. (Princeton University) 


Products of random matrices—Jn English 


10.1 (1.5) 


Ann. Math. Statist. (1960) 31, 457-469 (4 references) 


The authors consider a stochastic process X', X?, ... 
with values in the set of square matrices with real or 
complex entries of a fixed dimension and they investigate 
the asymptotic behaviour of the product "Y! = 
X"X""*...X". They generalise a result and prove a 
conjecture of Bellman [Duke Math. J. (1954) 21, 491- 
500] that these products obey a central limit theorem. 
By taking the norm | X || = max )’| X;,|, where the 
U J 


subscripts refer to the entry of the matrix, it is proved 
for a stationary process that the limit as no of 
tq = I/n{é log ||" Y* ||} exists; call it E. If, additionally 
the process is metrically transitive and 


*@ max {log || X? 


|, 0}<0o 


then also £8, = 1/n || log” Y* || +E. The limiting result 
for «, is proved straightforwardly. The result for £, 
is obtained in two stages. First, showing that the 
lim sup §, SE and the second is accomplished by con- 
sidering an auxiliary process (X,, Z,), (X2, Z2), +--+, 
where X, are matrices in the range of the X” and Z,, 
are matrices satisfying a recursion condition involving 


- the norm and such that the projection of the first 
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FURSTENBERG, H. & KESTEN, H. 
Products of random matrices—Jn English 


(Princeton University) 


coordinate is a member of the sample space of the 
X-process. The idea is to obtain uniform integrability of 
max {log || X7Z* |, 0} with respect to a mixture of 
measures induced by translation on the (X, Z)-process. 
Thus a suitable subsequence {v;'! of these mixtures 
converges weakly to a stationary measure yw for the 
(X, Z)-process and its projection agrees with the given 
probability measure on the X-process. The main step 
of the proof being to prove the relation 


E = lim inf &(v,) | XY" || S&() log || X? ¥' || 


where 6(.) represents expection with respect to the 
measure in parenthesis. This yields lim inf £,>E. 

Next the authors prove a result on the asymptotic 
behaviour of the entries. Assume that the possible 
matrix values M for X' satisfy the two conditions 
M;;>0 and 1<max M;;/min M,,<C<oo, where the 
extrema are both over all i, 7, and that if Q is an event 
for the X-process not depending upon the first (m+n) 
coordinates, then 

| Pr{Q | X*, ..., X"}—Pr{Q} | < DA"Pr{Q}; 
where D and / are fixed positive constants and 4<1. 


continued 


10.1 (1.5) 
continued 


Ann. Math. Statist. (1960) 31, 457-469 (4 references) 


The theorem is: If &| log Xy; |?*’<oo for some 
6>0 then for a, b (their defining equations are given) 
we have, if b>0, that {log ("Y’);,;—na}/,/(nb) converges 
in law to I1(0, 1) and if b = 0 then {log("Y'),;—na}/,/n 
converges in probability to 0. The proof largely follows 
the treatment in Doob’s book, [Stochastic Processes 
(1953) New York: Wiley]. The paper concludes with 
two examples indicating directions in which generalisa- 
tions of these results are not possible. 


(S. C. Saunders) 
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GNEDENKO, B. V. (Moscow University) 


On some aspects of the theory of queues—In German 


10.4 (-.-) 


Math. Tech. Wirtschaft (1960) 4, 162-166 (14 references) 


This article presents a review of some results of the 
theory of queueing, with special regard to some more 
recent models which have been dealt with by Russian 
mathematicians. 

After a general introduction in which the author 
emphasises the inter-dependence between theory and 
practice in the field of mathematics, with quotations by 
Aristoteles, Tchebycheff, Gauss, etc. he deals with the 
theory of queueing. The models as considered by 
Gnedenko refer, above all, to the following generalisation 
of the simple problem: the units entering the system 
are not willing to accept unlimited waiting time, but 
may possibly leave the system without waiting for 
completion of the service or without even entering the 
service station. Some more simple modelsas developed by 
Barrer, Takacs, Jaroschenko and Farinitsch are presented, 
in which various aspects of deserting are considered. 

The author deals with a comprehensive model by 
Kowalenko, which contains all the cases already con- 
sidered as special cases. The essential novelty suggested 
by Kowalenko results from the fact that the time which 
ihe individual units are prepared to spend on the service 
proper is a random variable, the distribution of which 
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GNEDENKO, B. V. (Moscow University) 


depends on the actual waiting time beginning with the 
time of joining the queue and ending with the entrance 
to the service point. Further, the time which the units 
are willing to wait in the queue is a time-independent 
random variable. 

Assuming exponentially distributed service-time and 
Poisson-input into the queue, an integro-differential 
equation for F(x, t) is established. F(x, f) indicates the 
probability that all units arriving before the moment t 
have left the system in the interval (t, x+1). It can be 
shown that in the model of Kowalenko F(x)= lim F(x, 1) 

t7o 
exists. With the aid of the function F(x) the most 
important characteristics with respect to the utilisation 
of the service-procedure can be expressed; for example 
the distribution of waiting time, of the total time spent 
in the system, the probability of loss (that is the case of 
the unit departing without completion of service). 

Kowalenko’s solution is illustrated by an example 
applying the simple model of Barrer in which the waiting 
time spent in the queue has a fixed limit. All the models 
presented are restricted to service systems with a single 
channel. 


(F. Ferschl) 


10.4 (2.5) 


On a generalisation of Erlang’s formulae—In Ukrainian 
Dopoy. Acad. Nauk. Ukrain, SSR (1959) 347-350 (4 references) 


The paper gives a generalisation of Erlang’s formulae 
for the case when devices serving demands may become 
out of order and require expenditure of the operator’s 
time for their repair or replenishment. The case is 
considered when the demands are served by n devices, 
and these devices by r operators. The inflow of demands, 
as well as the flow of impairments of the devices are 
assumed to be simple Poisson flows. The probability 
distribution of the time of service of both the demands 
and the instruments is assumed to be exponential. 


(B. V. Gnedenko) 
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GRIFFING, G. 


(Commonw. Sci. Industr. Res. Organisation, Canberra) 


10.9 (-.-) 


Accommodation of linkage in mass selection theory—In English 
Aust. J. Biol. Sci. (1960) 13, 501-526 (10 references, 1 table) 


The purpose of this paper is to develop mass selection 
theory which will accommodate not only linkage but 


_will provide for different recombination frequencies in 


the two sexes. 
The theoretical aspects of the linkage problem are 
developed in three stages: 


(i) The mass selection theory for two loci is extended 
to accommodate different recombination values 
for the two sexes. 

(ii) A method is developed by which the generalised 
two-locus model may be used to cope with genetic 
situations which are considerably more complex. 
This method requires the estimation of the 
recombination value averaged over all possible 
pairs of loci. 

(iii) The expectations of the half-sib and full-sib 
covariances for a random-mating population are 
generalised to permit different recombination 
values for the two sexes. This allows unbiased 


estimates of genotypic variance components to be. 


obtained. 
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HIDA, T. (University of Kyoto, Japan) 


Finally, application of the more general mass selection 
theory to the problem of detecting the influence of 
natural selection in modifying the effectiveness of 
artificial selection, is discussed. 


(G. Griffing) 


10.1 (1.0) 


Canonical representations of Gaussian processes and their applications—Jn English 
Mem. Coll. Sci., Kyoto (1960) 33, 109-155 (22 references) 


This paper gives systematic establishment of the Lévy’s 
theory, which is concerned with the canonical represent- 
ation of Gaussian processes, and a development of 
some new facts. In section one, the general theory of 
representation is introduced. After showing, as Lévy 


_ proved, that a canonical representation of any Gaussian 


process is uniquely determined if it exists, the author 
gives a necessary and sufficient condition for the process 
to have a canonical representation; using Hellinger- 
Hahn’s Theorem in the theory of Hilbert Space. 

The author gives a criterion to determine whether a 
given representation is canonical or not, which is a 
generalisation of Karhunen’s theorem on a stationary 
process. This criterion is different from Lévy’s, which 
is given with the help of Hellinger integral. 

In section two, multiple Markoff Gaussian processes 
are discussed by using the general theory of representa- 
tion. In this paper the author generalises the notion 
of N-ple Markoff processes, which were defined by 
Doob or Lévy, in order to treat the processes which 
are not always differentiable. It is shown that the 
canonical kernel of the N-ple Markoff process is a 


Goursat kernel of order N, which generalises the fact 
obtained by Lévy, and special N-ple Markoff Gaussian 
processes which are differentiable up to N—1 times are 
discussed. As the generalisation of Doob’s Theorem, 
it is proved that the stationary N-ple Markoff process 
is the sum of general simple Markoff processes. Further, 
exact forms of predictor are given for certain classes of 
multiple Markoff processes. 

Section three is concerned with Lévy’s M(t) process. 
Let M,(t) be the spatial average of an ordinary Brownian 
motion with N-dimensional parameter over the sphere 
with centre O (origin). Lévy obtained the canonical 
representation of this process only when WN is odd. 
Simplifying his proof by transforming M),(t) into a 
stationary process, the author gives the canonical 
representation of M(t), which can be applicable even 
when JN is even. It is also shown that the number of 
different representations of M,p,,(¢) is just 2°73 
including canonical one. Further, it is proved that 
M,,(t) is not a multiple Markoff process, though 
M,p.,(t) is a (p+1)-ple Markoff process. 


(M. Huzii) 
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KEMENY, J. G. & SNELL, L. (Dartmouth College, Hanover, U.S.A.) 10.1 (0.6) 


Finite continuous time Markoff chains—In English 
Teor. Veroyat. Primen. (1961) 6, 110-115 


In a recent book [Finite Markoff Chains (1959) Princeton: 
Van Nostrand] the authors developed a method for 
computing the basic parameters of a finite Markoff 
chain. Matrix expressions were obtained for these 
quantities in terms of certain basic matrices easily 
obtainable from the transition matrix. In this note 
corresponding expressions are given for the basic 
quantities for a finite continuous time chain. 


(B. V. Gnedenko) 
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KHASEN, E. M. (Moscow) 


10.1 (10.2) 


Evaluation of the single-dimensional probability density of random process 


in the output of a nonlinear system—Jn Russian 


Teor. Veroyat. Primen. (1961) 6, 130-138 (10 references) 


In this paper the author considers nonlinear systems 
with feed-back. These systems contain nonlinear 
functional reformers, which are separated by linear 
filters with rational spectra. 

In order to evaluate the single-dimensional probability 
density and the moments of a random process in the 
output of the system, the author constructs a multi- 
dimensional random Markoff process; and employs 
the Kolmogoroff diffusional equations. 

The exact results obtained by using this method are 
compared with the results of an approximate “ statistical 
linearisation ” approach for one system. The two sets 
of results are found to agree well. 


(E. M. Khasen) 
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KRZYZANSKI, M. (Math. Inst., University of Krakéw) 10.1 (10.8) 


The diffusion under the action of an exterior force and the Kolmogoroff 
equation—IJn French 


Second Hungarian Mathematical Congress (1961) (2 references) 


A sequence of stochastic processes relative to a diffusion the differential equation of the diffusion under the 
is considered where an exterior force expressed by a_ action of an exterior force. This equation was obtained 
function f(x) of the abscissa X is in action. Systems of by Smoluchowski independently from the above formula. 
registration, moments of the abscissa of a moving A theorem of Khintchine, given in his work [Asymp- 
particle 


correspond to these processes. Let W(x, y, A), 
where Aj” = tf?,—t{”, be the probability density of 
The function W,(x, y, AM) is assumed 
to be given by the formula of Smoluchowski [Annalen 


the passage. 


BO 4) ef 


hy 


der Physik (1915) 48, 1103-1112] 
(i) W(x, y, Af”) 


—x— BAS/C? 
= [4nDA®1-# exp 2 LY =x = PAR AOD]’ 
L yal P| ADA 


totische Gesetze der Wahrscheinlichkeitsrechnung (1933) 
Berlin], is applied in the proof. This required the 
verification of some theorems concerning the properties 
of the solutions of the parabolic equations. 


(M. Krzyzanski) 


Under certain hypotheses it is proved that the limiting 
density for noo, max A{”)—=0 satisfies the equations 
of Kolmogoroff; where the second one is identical to 
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LEONOFF, V. P. (Moscow) 


On the dispersion of time-dependent means of a stationary stochastic 
process—In Russian 


Teor. Veroyat. Primen. (1961) 6, 93-101 (9 references) 


10.1 (-.-) 


Let €(t) be a stationary process in the wide sense with 
discrete (continuous) time Mé(t) = 0 


op = yy &(1) (: = le cod), 


b,=M|C,|’. 


In this paper the author studies the behaviour of b, for 
poo. 


(V. P. Leonoff) 


MORAN, P. A. P. (Australian National University, Canberra) 


10.9 (-.-) 


The survival of a mutant gene under selection. Il.—Jn English 
J. Aust. Math. Soc. (1960) 1, 485-491 (2 references) 


Suppose M haploid individuals are of genotype a or 
A, with a having a small selective advantage over A 
so that the expectations of offspring are in the proportions 
1+s:1 respectively. The paper shows that if the 
population consists initially of kg individuals a, the 
absorption probability P,, that the total population 
eventually becomes of type a is approximately 
Pag = (1— 27 2%)/(1 = 0-29), 

where s/(1+s)S0<s. For k, = 1, this reduces (when 
s*M is small) to Fisher’s result 2s/(1—e~?*™), [see The 
Genetical Theory of Natural Selection (1930), Oxford: 
Clarendon Press ]. 

The argument used for this approximation is extended 


to the case of diploid individuals with distinct sexes, 
and a similar absorption probability is approximated. 


(J. Gani) 
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MOTOO, M. (Institute of Statistical Mathematics, Tokyo) 
+ b'(x) aos —In English 
Ox' 


a 


Bos 0 
Diffusion process corresponding to 52 ; 


x/2 


10.0 (1.0) 


Ann. Inst. Statist. Math., Tokyo (1960) 12, 37-61 (5 references) 


The author first constructs the diffusion process corres- 
ponding to the operator 
y) N 


ahaa an 
pm Oa, 


by using the stochastic integral. The method of con- 
struction is the multi-dimensional generalisation of that 
which is treated in Maruyama’s paper. The author 
shows that the semi-group associated with the process 
can be restricted to certain sub-spaces of continuous 
functions and it is strongly continuous on these spaces. 

Assuming the Lipschitz condition for b'‘(x)’s, the 
Hille-Yoshida’s generator of the semi-group of the 
process is investigated. Dynkin’s generator of the 
process is given in the form independent of the process. 
This gives the exact meaning of the correspondence 
discussed in constructing the process, and will serve to 
solve the problems of the differential operator A by 
the probabilistic method. . 
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Finally, using the fact that the process constructed 
above and the Brownian motion are absolutely con- 
tinuous with respect to each other, the author gives 
some probabilistic properties which are easily derived 
from those of the Brownian motion. 


(M. Motoo) 
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MUTH, J. F. 


(Carnegie Institute of Technology, Pittsburg) 
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10.2 (6.1) 


Optimal properties of exponentially weighted forecasts—In English 
J. Amer. Statist. Ass. (1960) 55, 299-306 (12 references) 


For time series of the form 


MV ifeaa Ep » opis 


or of the form 
Vi=Vitn, 


t 
Y, SF Ds gj 
heme 
the ¢, are independently distributed with mean zero 
and variance a? and the n, are independently distributed 
with mean zero and variance o%, it is demonstrated that 
predictions for the period t+T of the general form 


where 


vier = BY A=, 


are appropriate. 
Assuming that « is given by 
C, = ay,+0, 
6, distributed independently of jy, and y{, an estimate of 


2575) 


ONICESCU, O. (University of Bucharest) 


a may be obtained by replacing j, by the prediction for 
period ¢+1, yfi1, and by oe the ordinary least 
squares formula, that is 
a = UCVr+ - 

X(yi+1) 
The bias in this estimate is investigated and shown to 
be small in large samples. 


(W. A. Fuller) 


10.5 (6.2) 


The entropy-function, correlation and quantity of information—Jn French 


Second Hungarian Mathematical Congress (1961) 


Let (Q, 7, P) be a probability space; Q is an arbitrary 
space, & a o-algebra of subsets of Q(Q € #) and P(A) 
a probability measure defined for Ade .#. Let © be 
a partition of Q into a finite number of disjoint sets 
01, 05, ---, 0, (0, € % Pr(O,)>0). The entropy function 
H,(A) of the partition © is defined as follows: H(A) 


is a o-additive set function defined for A € by 


n 


y) Pr(6,A) log 
k=1 


ee Pr (0,). 

Here the product of two sets denotes their intersection. 
Clearly H,(Q) is equal to the usual entropy of the 
partition ©. It is shown that the most important 
properties of the entropy remain valid for the entropy 
function too, for example, if 6’ =(0{, 6, ..., @,,) and 
0” = (01, 03, ..., 0”) are two independent partitions and 
© = 0’. 0” denotes the product of the two partitions, 
that is the partition formed by the sets 05.0; then we 
have identically for A € A 


H(A) = Ho(A)+Ho(A). 
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The quantity of information is generalised similarly by 
defining the information-function Jy, (A) of the partition 


Oa (075605, 0.) es With Hespect to. the. partition 
d = (1, bo, ---» bm) Which is also a set function, and 
is expressed by the formula 
M m Pr (0 Dp ) 
I, 4(A)= PIO ,O Ay Oe 
Fie 2 Beit A) los Pr (0,)Pr(,) 


Of course J, 4(Q) is equal to the usual quantity of 
information. The information function for the general, 
not discrete, case is defined by passing to the limit. 


(A. Rényi) 


EpitoRIAL Note: The original paper used abstract 
Boolean algebra for the purpose of defining the entropy 
and information functions. 
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RENYI, A. (Math. Inst., Hungarian Academy of Sciences, Budapest) 


10.5 (4.0 


Some fundamental problems of information theory—In H ungarian 
Magy. Tud. Akad. III Oszt. Kézl. (1960) 10, 251-282 (24 references) 


In this paper the different possibilities of measuring the 
amount of information are considered. The author 
investigates whether it is possible to find quantities other 
than Shannon’s measure of information. It is shown 
that the quantities satisfying certain reasonable postulates 
concerning the amount of information are only the 
so-called measures of information of order « contained 
in the value of the random variable having the discrete 
distribution P = (p,, Pp, ..., p,), which are defined by 
I{P) = 1/1—«) log, Zp. The limiting case for «71 
gives Shannon’s measure of information, which is 
characterised within the family of quantities of inform- 
ation /,(P) by certain of its properties. 

The author investigates the possibilities of extending 
some known theorems concerning Shannon’s measure 
of information to the measure of information of order 
a #1. He deals with the measures of information 
concerning the so-called “‘ generalised incomplete- 
probability distributions ’’ introduced by him; that is, 
distributions of random variables which are not certainly 
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SADAO, I. 


(College of Science & Engineering, Nihon Univ., Japan) ~ 


observable as well as with the relative information of 
order «. The different definitions presenting themselves 
concerning the measure of information of order « 
contained in the observation of a random variable ¢€ 
relative to a random variable 7 are investigated; and 
taking into account the advantages and disadvantages 
of these different definitions, it is shown which of them 
is the most suitable measure of the relative information 
of order a. 

The quantities of information of order «, discussed in 
this paper have, in several respects, similar properties 
to the Shannon’s measures and they have certain 
advantages too; for example they are often easier dealt 
with. The application of the notion of information in 
statistics is also discussed. It is shown that when 
dealing with statistical applications the measures of 
information of order ~ + 1 lead to the same conclusions 
as Shannon’s measure of information; for example 
they are in the same relation with Fisher’s “‘ quantity 
of information ”’. 


(K. Bognar) 


10.5 (-.-) 


A remark on the convergence of Kullback-Leibler’s mean information—ZJn English 
Amn. Inst. Statist. Math., Tokyo (1960) 12, 81-88 (2 references) 


Under fairly general assumptions, the author gives two 
sufficient conditions for the convergence of the mean- 
information for finite measure. These are generalisa- 
tions of the conditions given in Lemma 2.1 and Theorem 
2.2 of Chapter 4 in Kullback’s book [Information Theory 
and Statistics (1959) New York: Wiley]. 

In this paper the author considers the conditions 


ess.sup | 1— t | —0, instead of the condition 


lim (4) S41 
I 
uniformly which was given by Kullback. The author 
says that, in general, the conditions in Lemma 2.1 and 


the finiteness of I(f:g) in Theorem 2.2 of Kullback’s 
book are not necessary. 


(M. Huzii) 
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RENYI, A. (Math. Inst., Hungarian Academy of Sciences, Budapest) 


On different measures of informaticn—In English 
Second Mathematical Congress (1961) 


Let [Q, /, P] be a probability space. A generalised 
random variable is defined as an ./-measurable function 
on a sub-set Q; of Q (Q, € ¥, Pr(Q,)>0). This can be 
interpreted as the result of an experiment depending on 
chance which is not always observable, only for w €Q,. 

Let ¢ be a generalised random variable taking on 
only a finite number of values x,, k = 1, 2, ..., n and 
B an arbitrary event (Be &) with positive probability. 
Put p, = Pr(c, = x), d = Pr(é&, = x, |B) (p,>09), 
P= (P1, Fils Pn)s Q a (41, Gimp In): 

The amount of information concerning €, gained by 
observing the event B is denoted by I(Q | P); it is 
supposed that this information depends only on Y and 
Q. Six postulates are given, such that any reasonable 
measure of information-gain should fulfil all of them. 
It is shown that if 1(Q | Ff) satisfies these postulates, 
we have either 


1(Q |P) = 1[(a—-1)log, » («tle 1a.) = 1,(0 |) 
(a + 1) 
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RENYI, A. (Math. Inst., Hungarian Academy of Sciences, Budapest) 


On different measures of information—In English 
Second Mathematical Congress (1961) 


An additional postulate is given, which characterises 
I,(Q || P) among all 1,(Q || Y). 


The “‘ entropy of order « ” of the generalised distribu- 
tion FP = (p;, ---, p,) is defined in accordance with the 
above formulae. If € is an ordinary random variable, 
then the expression reduces to the well-known formula, 
of Shannon for the entropy of a probability distribution. 
The case of continuous random variables is also dis- 
cussed together with the connection of entropy of order 
a with kinetic theory. 


(I. Csiszar) 


ee 


10.5 (4.0) 


1(Q |) = {Sactoetalnot | Ya =1,(0|P) 
= lim 1,(Q ||P) 


1,(Q || P) is called the “ information-gain of order «” 
concerning the generalised random variable € from 
the observation of the event B; or, from replacing 7 
by Q. If € is an ordinary random variable 


(Sn.= I, Fa, =1) 


then J,(Q | FY) reduces to the well-known formula of 
information gain used in the theory of information. 

It is shown that in many respects 1,(Q || P) behaves 
similarly for all positive values of «. Fisher’s measure 
of information can be deduced from 1,(Q || A) in the 
same way as from J,(Q || A). The same remark applies 
to the distance measure of Mahalanobis, the 77-distance, 
etc. 


continued 


10.5 (4.0) 
continued 
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SEREGIN, L.. V. (Moscow) 


10.1 (-.-) 


Continuity conditions of stochastic processes—In Russian 
Teor. Veroyat. Primen. (1961) 6, 3-30 (8 references) 


Let X, OSt<c<oo be a separable stochastic process 
in metric space X. The author states that the main 
purpose of this paper is to derive conditions under 


which almost all sample functions of process x, are 


continuous. We designate p(x, y) the distance between 
the points x, J eX. 
Let P(...) be a Markoff transition function, satisfying 

for each e>0 

Sup Pr (5) x dx) =-0(1), et 0, 

x, S, t 
where xe X; s, te[0, cl], O<t—s<h and V,(x) = 
{y : p(x, y)2e}. Then almost all sample functions of 
Markoff process x, are continuous if and only if for 
each e>0 


c—h 
| Pr P(%uria) > ehat = o(h), hd 0; 


0) 
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SEGUCHI, T. 


Some statistical process controls. 
II. Controls of finite Markoff chains—Jn English 


(Kyushu University, Japan) 


Almost all sample functions of a Martingale (semi- 
martingale) x, are continuous if and only if for h | 0 


c—h 
| Pin, ax, ye at = O(n), 


(0) 


ch 
| Pr (>), X,44<a}dt = oh) 
0 
for each a and b, (a<b). 


(L. V. Seregin) 


10.1 (8.9) 


I. Controls of additive processes, 


Mem. Fac. Sci. A, Kyushu (1960) 14, 139-178 (11 references) 


In part I of this paper which occurs between pages 139 
and 153, the author considers the controlling system 
of a discrete time parameter additive process. The 
system works to bring the process back to the reference 
value when the process shows a deviation larger than 


-a preassigned constant from the reference value. The 


action of the system may contain a random error which 
is independent of all the other variables. The author 
gives the stationary distribution of the Markoff process 
representing the process under the control by this 
system: this section includes six references. In part II 
which follows (on pages 154 to 178) is given a controlling 
system for a Markoff chain. The system is supposed 
to work in such a way as to make the chain return to 
some preassigned state when it reaches some unfavour- 
able state. After the detailed descriptions of the necessary 
quantities, the author gives a proof of the fact that the 
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best controlling system or the system which gives the 
least stationary probability of the process being at some 
unfavourable state, is the one that makes the process 
return to the state which is the farthest, in the sense of 
the mean-length of the first passage time, from the 
boundary or the set of unfavourable states. 


(H. Akaike) 
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dition is very often not fulfilled. 
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SOEST, J. L. van (Technical University, Delft, Netherlands) 


Order and disorder—In English 


10.0 (10.5) 


Statist. Neerlandica (1960) 14, 249-258 (15 references) 


The degree of order and of disorder are defined such 
that they are complementary and, both, vary between 
0 and 1. Examples are given in problems of human 
life, and of human community with anarchistic and 
totalitarian extremes; also in thermo-dynamic entropy 
and information theory. In the latter problems the 


statistical quantity )° p(i) log p(i) appears; it has 
ih 


quite the same meaning in both cases: the quantity of 
statistically weighted alternatives. 

The relation to entropy and to quantity of selective 
information is discussed. The degree of order 1, degree 
of disorder 0, appears if this statistical quantity is 0, 
one of the p(i) being equal to 1: “frozen” or ‘‘ pure ” 
state. The degree of order 0, degree of disorder 1, 
appears if this quantity is log n, all p(i) being equal to 


1/n. The degree of disorder is equal to }) p(i)” log p(i). 
ta 
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STORMER, H. (Siemens u. Halske A.G.,-Miinchen) 


In information theory the degree of order is identical 
with redundancy. Order/disorder problems in patterns 
are related to Markoff chains of higher order and are 
not discussed here. 


(J. L. van Soest) 


10.4 (-.-) 


A queueing problem in telephone exchanges—/n German 
Zeit. angew. Math. Mech. (1960) 40, 236-246 (3 references) 


The well-known distribution, first given by Erlang, of 
waiting times occurring when the central steering 
mechanism is working with a constant holding time 
Ty, applies to the case where the calls are handled in 
the order in which they arrive. In practice, this con- 
Instead, the waiting 
calls are dealt with according to chance. 

In the present investigation the distribution of waiting 
calls is first treated under the assumption that the calls 
are being handled according to a system of priorities 
possessing m different steps. It is then assumed that 
there is a chance correlation between any one waiting 
call and one of these priority steps. In the limit m—oo 
a simple closed expression is obtained for the distribution 
of the waiting times. > 


(H. Stormer) 
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STAM, A. J. (Physical Lab., Gov. Def. Org., The Hague, Netherlands) 
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10.0 (10.5) 


Some mathematical remarks on information theory—In Dutch 
Statist. Neerlandica (1960) 14, 259-265 (18 references) 


This paper is the elaboration of a lecture given recently 
by the author. It contains definitions of entropy and 
the amount of information as given by Shannon; some 
remarks on their interpretation and their generalisations 
to abstract probability spaces. Some remarks are also 
made on the mathematical problems arising from these 
definitions. Shannon’s definition is compared with the 
concept of “intrinsic accuracy’ or information used 
in statistics that was defined by Fisher. Some properties 
which these different kinds of “ information” have in 
common are mentioned. 


(A. J. Stam) 
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THEODORESCU, R. (Institute of Mathematics, Bucharest) 


Multiple Markoff processes—Jn Russian 


10.1 (10.8) 


Rey. Math. Pures Appl. (1960) 5, 341-362 (14 references) 


. In this paper a direct study of multiple Markoff processes 


is suggested by means of integro-differential equations. 
Starting from the functional equations of a multiple 
Markoff process [see Onicescu, Probability theory (1956) 
Bucuresti: Ed. Tehnica] which are replacing the 
Chapman-Kolmogoroff equation, in section one the 
corresponding differential equations are deduced for 
the case of a finite set of states. These equations are 
analogous to those given by Kolmogoroff for simple 
Markoff processes. 

In section two multiple Markoff processes on the 
real axis are considered; the corresponding integro- 
differential equations are deduced and certain ergodicity 
properties are investigated. Further, conditions for 
existence and unicity for the corresponding solutions 
are provided. 

Finally, the case leading to systems with partial 
derivatives similar to the Fokker-Planck equation, are 
also investigated. 
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The present paper represents a shortened Russian 
version of a sequel of papers published previously in 
Rumanian [Bul. stiintific, sect. st. mat. fiz. (1955) 7, 
763-774 and 776-794; Analele Univ. C. I. Parhon, ser. 
gt. naturii (1956) 10, 23-24]. 


(M. losifescu) 
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THEODORESCU, R. (Institute of Mathematics, Bucharest) 10.1 (10.8) 


| Certain stochastic models for learning and their interpretation. 
Linear chains—ZJn French 


Second Hungarian Mathematical Congress (1961) (6 references) 


A systematic treatment of stochastic models occurring 
in learning theory was given by Bush & Mosteller in 
their monograph [Stochastic models for learning (1955)]. 
In the present note it is shown that certain stochastic 
models for learning, considered by these authors, may 
| be obtained by means of schemes involving chains with 
| complete connections, introduced by Onicescu & Mihoc 
| [see Bull. Sci. Math. (1935) 59, 174-192]. A recent 
: monograph deals with the detailed treatment of these 
processes [see Ciucu & Theodorescu Processes with 
complete connections (1960) Bucharest: Ed. Acad. 
R.P.R., in Rumanian]. 


(R. Theodorescu) 
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VOLKONSKY, V. A. (Moscow) 10.1 (-.-) 


Construction of inhomogeneous Markoff processes by means of a random 
time substitution—IJn Russian 


Teor. Veroyat. Primen. (1961) 6, 47-65 (6 references) 


It is proved by the author in this paper that a continuous 
single dimensional Markoff process y(t) with wide 
restrictions can be obtained from the Wiener process 
x(t) in the following form: y(t) = W[x(t,), t], where 
W(x, tf) is a continuous function, monotonic in x for 
a given ft, and t, is a nondecreasing random function 
of ¢ (theorem 1). 

The author gives conditions which should be met by 
the Markoff process x(t) in abstract space and the 
random function t, so that the process y(t) = x(t,) will 
also be a Markoff process (theorem 2). 


(V. A. Volkonsky) 
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WILKINS, C. A. 


WIGGINS, A. D. 


(Hanford Laboratories Operation, G.E.C., Richland, Washington) 


10.2 (10.9) 


On a multicompartment migration model with chronic feeding—Jn English 


Biometrics (1960) 16, 642-658 (6 references) 


A biological system with K compartments is considered, 
and the probability of a particle migrating from com- 
partment 7 to compartment j in a small time, h, is 
4; ;h+0(A), where the A;,; are constant migration rates. 

A set of simultaneous differential-difference equations 
is derived for the transition probability of a particle 
starting in compartment i = 1 having migrated to com- 
partment j at time ¢. When particles are introduced 
into the first compartment according to a known 
* feeding ’ function, the state of the system is described 
by the number of particles in the various compartments. 
By using generating functions, the stochastic behaviour 
is found in terms of the single-particle process. In 
particular, it is shown that the expected numbers of 
particles in the compartment satisfy the same set of 
equations as the original probabilities, except for the 
addition of non-homogeneous terms. 

If experimental determinations of the numbers of 
particles in each compartment are made at discrete 
times, these may be used to estimate the migration rates 
d jj. 
replaced by the observations and their time differences 
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(University of New South Wales) 
On two queues in parallel—In English 
Biometrika (1960) 47, 198-199 (1 reference) 


In a recent paper Haight, F. A. [Biometrika (1958) 45, 
401-410; abstracted in this journal No. 1/128, 10.4] 
investigated a two-queue system in which an arrival 
is assigned to the shorter queue, or if they are of equal 
length, to a specified one, the “‘ nearer ”’ one. 

The results given by Haight are here briefly extended 


to the general case where, if X is the length of the near 


queue, Y the length of the other, and W(XY) the 
probability of an arrival joining the near queue, 


G) W(x, y) = 1, x<y 
Gi) W(x, y) = w(x), x = y 
(ili) W(x, vy) = 0, x>y. 
The modifications of Haight’s results are set out by the 
present author. 


(Florence N. David) 
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The expectations and their time derivatives are 


in the set of equations, and a method similar to least- 
squares autoregression fitting is used to obtain the 
estimates. To facilitate the solution, an isomorphism 
is constructed between certain matrices of a block- 
pattern type and matrices of a lower order; multiplica- 
tion of a matrix of the first type by one of the second is 
defined. The inversion of a block matrix is required 
for the solution, and this is achieved by proving that 
each block is a cyclic matrix whose inverse is therefore 
a power of itself. 

Three unsolved problems are posed. In view of the 
randomness of the observations, do the linear estimating 
equations have a solution with probability one? What 
are the sampling properties of the estimators? Would 
their accuracy be greater if the observations were 
spaced at shorter rather than longer time intervals? 
On the basis of an artificial example with known 
migration rates, it is suggested that the answer to the 
third problem is in the affirmative. 


(G. A. Watterson) 
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YEO, G. F. (University of Western Australia) 10.4 (-.-) 
The time-dependent solution for a dam with geometric inputs—In English 
Aust. J. Appl. Sci. (1960) 11, 434-442 (6 references) 


Let the content Z, = 0, 1, 2, ..., of an infinite dam at The stationary probabilities P; = lim P,t) of the 


the discrete times ¢ = 0, 1, 2, ..., be defined by . : ee Oe 
process give results previously obtained by Prabhu [ Ann. 
Zi44 = Z,+ X,—min (Z,+ X,, 1), Math. Statist. (1958) 29, 1234-1243; abstracted in this 


where X,=0, 1, 2, ..., is a random additive input journal No. 1/310, 10.4}. 

during the interval (¢, t+1); the process Z, is known 

to be a Markoff chain. This paper is concerned with (J. Gani). 
the explicit time-dependent solution for the probability 


Pine a). 1 = Oy. 1) 2... 
the initial content u being a non-negative integer when 
X, follows a geometric distribution. 

Using occupancy methods, the probability P,(f) of 
emptiness of the dam is first found. The generating 
function P(s,t)= }) P((t)s' (|s|<1) is shown to 

i =.0 


depend on P,(t), and the author makes use of residual 
methods to obtain explicit expressions for the P,(t) 
from it. 
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ZITEK, F. (Math. Inst., Czechoslovak Academy of Sciences, Prague) 10.0 (0.5) 
The theorem of Saks concerning the functions of random intervals—Jn French 
Second Hungarian Mathematical Congress (1961) (1 reference) 


In the theory of the Burkill integral a theorem of Saks 


is known stating approximately the following: if f 


is an integrable, in the sense of Burkill, and derivable 
interval function; then its integral is also derivable 
and the two derivatives are almost everywhere equal. 

The aim of the present note is to show that an 
analogous theorem holds for the functions of random 
intervals. The theory of these functions and of the 
(BB) integral, which is a generalisation of the Burkill 
integral, was developed in detail in the author’s papers 
[Czech. Math. J. (1958) 8, 583-609; abstracted in this 
journal No. 1/150, 10.0] and a further paper in the press. 
See also Cas. pést. math. (1959) 9, 83-89; abstracted in 
this journal No. 1/321, 10.0]. 


(F. Zitek) 
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Sa Ar. aw. 


HOSTELET, G. (41 Avenue W. Churchill, Bruxelles) 


11.0 (0.0) 


_ The co-operation of mathematics in the measurement and interpretation of 


characteristics of statistical phenomenae—In French 
Bull. Int. Statist. Inst. (1960) 37, III 375-384 


Mathematical statistics mean the co-operation of mathe- 
matics to investigate the measurement and the experi- 
mental analysis of specific statistical groups. An 
explanation is given of four out of the six main tasks of 
mathematical statistics with special emphasis on the 
relations between mathematical and experimental analysis. 

In the first part the empirical character of the algebraic 
fitting of frequency polygons or of time series is 
emphasised. The second step consists of obtaining, 
with the help of the algebraic smoothing formula, the 
algebraic expressions of significant characteristics of a 
statistical feature, where these characteristics had 
beforehand been defined by experimental analysis. The 
quantitative conditions of the validity of the experi- 
mental signification of these algebraic expressions are 
explained. 

The third step is the discovery of new characteristics 
by the experimentally explanatory analysis of the 
equation of the frequency curve or time series. Here, 
once more, the quantitative conditions of the experi- 
mental validity of new characteristics are explained. 
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HOSTELET, G. (41 Avenue W. Churchill, Brussels) 


The co-operation of mathematics in the measurement and interpretation of 


characteristics of statistical phenomenae—ZJn French 
Bull. Int. Statist. Inst. (1960) 37, III 375-384 


the n statistical observations are compared with n 
drawings in a lottery system of which the composition 
must be adapted to the studied data. Attention is 
drawn to the quantitative and qualitative conditions 
which ensure the experimental validity of the lottery 


system used and of the results derived from it. An 


example of the determination of significant character- 
istics in statistical observations is analysed. 

The fifth task of mathematical statistics: the evalua- 
tion of probable degrees of approximation of calculated 
statistical characteristics; and the sixth: the working 
out of sampling methods have not yet been dealt with. 
The paper aims to explain the quantitative conditions 
of the experimental validity of formulae and equations 
relating to the measurement of socio-statistical data, 
and states that these conditions should be more often 
taken into consideration. 


(G. Hostelet) 
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The other tasks imply the application of lottery systems: 
the author states that he has discriminated between 
two kinds of probabilities: 


(i) the probability of the existence or the production 
of a statistical fact; 

(ii) the probability of the exactitude of the measure- 
ment of a statistical fact. 


The fourth step is the search for measurements of 
observations which do not conform with the formulae 
of the lottery system, to which the total of observed 
statistical characteristics has been referred. It will 
thus be a matter of discovering the experimental signifi- 
cance of this disagreement. 

It is known that a group of observations relating to a 
statistical feature, for which the number of units is less 
than that of another group, could provide an estimate 
closer to the characteristic sought than the one provided 
by the other group. And without it being known. In 
statistical observations only probable degrees of accuracy 
can, therefore, be established. In order to obtain this 


continued 
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FREUDENTHAL, H. (Schubertstraat 44, Utrecht) 
Order and disorder—In Dutch 
Statist. Neerlandica (1960) 14, 235-247 


Order is a quite definite notion as long as it means a 
special kind of order. Disorder in the sense of absence 
of any order is a vague notion. Statistical notions do 
not depend on notions like order and disorder. Instead 
of disorder, independency is the fundamental statistical 
notion. Statistical independency can be obtained by 
creating disorder. But pseudo-random-material shows 
that it can also be obtained by suitable strictly causal 
methods. It is not possible to test whether a given 
material is produced by stochastic devices. It is danger- 
ous to try disorder as a null-hypothesis in order to 
refute it by showing some regularities in the underlying 
material. This is shown by the analysis of statistical 
sham proofs for the impossibility of generatio spontanea, 
such as tried by some biologists. 


(H. Freudenthal) 
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SMIRNOFF, S. V. & POTAPOFF, M. K. (Moscow) 
Nomogram for probability functions y7—Jn Russian 
Teor. Veroyat. Primen. (1961) 6, 138-140 (3 references) 


In this paper a nomogram is constructed for the function 
1 Sa ae oes 
D ss n—-1 ,—22/2 
Pr(x*, n) = paeeal. were dz 
or the variables P, x’, 1 lying within the following limits: 
1<n<110, 1<y*<150, 0:001<P<0-999. 


The relative error in the middle part of the answer 
scale of P does not exceed 3 per cent., for 0:15 P<S0°9, 
and 10 per cent. at the ends of this scale. 


(S. V. Smirnoff) 
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NEYMAN, J. (University of California, Berkeley) 


11.9 (10.0) 


Indeterminism in science and new demands on Statisticians—In English 
J. Amer. Statist. Ass. (1960) 55, 625-639 (14 references, 3 figures) 


The words “‘ indeterministic study ”’ are used to designate 
research aiming to determine how frequently a quantity 
X characterising the phenomena considered assumes its 
various particular values. If the purpose of research 
is to establish the exact value of X as a function of 
other variables, then this research is ‘‘ deterministic ”’. 
In the history of indeterminism in science four over- 
lapping periods are discernible. 

(i) Period of “‘ marginal indeterminism.” This was 
the period, symbolised by the names of Laplace and 
Gauss, in which research in science was all deterministic 
with just one domain, that of errors of measurement, 
treated indeterministically. 

(ii) Period of “ static indeterminism,”’ is symbolised by 
the names of Bruns, Charlier, Edgeworth, Galton and Karl 
Pearson. Here, the main subject of study was a “‘ popula- 
tion” and efforts were made to develop systems of frequency 
curves to describe analytically the empirical distributions. 

(iii) The third discernible period, may be termed the 
period of “static indeterministic experimentation.” It is 
marked by the name of R. A. Fisher and by his book 
The design of experiments. The typical problems con- 
sidered were: do these two populations have the same 
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distributions of X? This and similar questions led to the 
development of basicideas of tests of statistical hypotheses 
and of estimation, and also of the appropriate techniques. 
All of these are currently at the disposal and in constant 
use by an applied statistician. 

(iv) The fourth period in the history of indeterminism, 
currently in full swing, the period of “ dynamic in- 
determinism,” is characterised by the search for evolu- 
tionary ‘chance mechanisms capable of explaining the 
various frequencies observed in the development of the 
phenomena studied. The chance mechanism of carcino- 
genesis and the chance mechanism behind the varying 
properties of the comets in the Solar System exemplify 
the subjects of dynamic indeterministic studies. One 
might hazard the assertion that every serious con- 
temporary study is a study of the chance mechanism 
behind some phenomena. The statistical and prob- 
abilistic tool in such studies is the theory of stochastic 
processes, now involving many unsolved problems. In 
order that the applied statistician be in a position to 
cooperate effectively with the modern experimental 
scientist, the theoretical equipment of the statistician 
must include familiarity and capability of dealing with 
stochastic processes. (J. Neyman) 
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